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The Decay of Shock Waves 


B. CassEN AND J. STANTON 
Pasadena Physics Section, Naval Ordnance Test Station, Pasadena, California 


(Received July 2, 1947) 


The instantaneous rate of decay or intensification of a shock front in any fluid for which 
dissipative effects are localized in a shock zone is shown to be directly related to the coexisting 
pressure gradient in the wake of the wave adjacent to the front. The higher time derivatives of 
the shock velocity are also explicitly related, although in a more complicated fashion, to the 
higher spatial derivatives of wake pressure. Thus the experimental observation of the velocity 
or some related variable of a blast or ballistic shock wave can be used to ascertain properties 
of the fluid motion in the wake of the shock. Conversely, a knowledge of the wake profile at any 
instant can be employed to determine the corresponding rate of decay of the shock front. The 
general theory is applicable to shocks of all strengths where an equation of state can be applied 


graphically or analytically. 





INTRODUCTION 


T is commonly known that blast waves decay 
with time, not only on account of spatials 
divergence, but also because of the dissipative 
processes occurring within the shock front. Free 
energy is flowing toward the shock discontiniuty 
at a rate determined by the local pressure, par- 
ticle velocity, and the gradients of these quan- 
tities in the immediate wake. Ordinarily the rate 
of free energy degradation at the front, plus the 
rate of weakening caused by divergence, exceed 
the forward influx of free energy and the wave 
decays. But intensification may occur if the 
influx predominates. This latter event is rarely 
observed in blast waves. However, photographic 
evidence has been obtained in two instances of 
an intense blast which, emerging with strong 
spherical divergence from a small opening in a 
container, maintained nearly constant shock 
velocity for several hundred microseconds. 
During this interval, reinforcement from the rear 


compensated the divergence weakening although 


it was insufficient to cause intensification of the 
shock. 


THEORY OF FIRST DERIVATIVE 
RELATIONSHIPS 


The theory of the relationship of the rate of 
change of quantities pertaining to a shock front 
to the fluid motion in the wake is general and not 
dificult to develop. Relationships between 
second and higher derivatives are more involved, 
although still obtainable by straightforward 
algebra. The general theory of the first deriva- 
tives can be developed simultaneously for plane, 
cylindrical, and spherical shock waves by use of 
an index n, which is equal to 0 for the plane case, 
1 for the cylindrical, and 2 for the spherical. The 
theory also applies to plane oblique shocks; such 
shocks, as is well known, can be treated as 
normal shocks subject to a tangential velocity 
transformation. 
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Fic. 1. Wake pressure distribution to produce decay or 
intensification of plane waves. 


The equations of motion and continuity are 
written in Lagrange form as follows: 


dv/dt = —(x/z)"(0P/dz) ; (1) 
V = (x/z)"(dx/dz). 


Here v=particle velocity, x=spatial position 
or radius, P=p/po=pressure/initial density, 
V=po/p=initial density/density, z=Lagrange 
label of particles, and ¢=time. The adiabatic 
equation of state which obtains in the wake of 
the shock is expressed quite generally. If S is the 
entropy, let this be written as: 


V=V(S, P). (2) 


Time differentiation of the equation of con- 
tinuity gives 


(x/z)"(dv/dz)+ (nVv/x) =dV/dt 
=(0V/dP)(dP/dt)+(AV/dS)(AS/dt). (3) 


It is assumed that in the wake of the shock 
0S/dt=0, but not necessarily 0S/dz. It follows 
that the equations determining fluid motion in 
the wake are reduced to 


dv/dt= —(x/z)"(0P/dz); 


(x/2)"(dv/dz)+(nVv/x) =(0V/dP) (dP/dt). (4) 


These equations are most conveniently con- 
sidered in a coordinate system at rest with 
respect to the shock front in Lagrange labeled 
space, i.e., one for which the particle labeled z is 
replaced by Z=z—z, and T=¢. Here 2;(t) is the 
Lagrange position of the shock front, and 
U=dz,/dt is the shock velocity. 

Then 0/dz2=0/0Z (represented by dash super- 
script in what follows) and 0/dt+ U(d/dz) =8/dT 
(represented by dot derivative in what follows). 
After transformation, Eqs. (4) become 


: b— Uv’ = —(x/z)"P’; 5 
(x/z)"v’ +n Vo/x=dV/aP(P— UP’). (5) 


These two equations are linear in v’ and P’ and 
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Fic. 2. Wake pressure distribution to produce decay or 
intensification of diverging waves. 


can be solved explicitly for these quantities as 


z\"/0OV. nvVv 
ra|-+0f2) GP )I/ 
x oP x 
sy" OV s/2\" 
(=) +a(-) | 
z OP\x 
OV /2\" 
vA-) #|/ 
OP\x 
x\* aV/2\" 
(=) +752(-) | 
Zz OP\x 


When specialized to particular values at the 
shock front, where x =z=2), one obtains 


nV w,U 
8) /asRv 
21 
(7) 


(6) 


OV. nVv 
v-|—P-"—+ 
oP x 


P,’ =| —U(@y— URPy) — 


3 





* nV Ww 
v= U(RPy- URw)-"—"| /{1+R0%] 
21 


5 


where R=(0V/0P)z~» and U used as a subscript 
indicates differentiating with respect to U at the 
shock front. As will be shown, all quantities 
appearing on the right-hand side of these ex- 
pressions can be determined in terms of U and U, 
when the equation of state of the fluid is given. 

For a steady plane wave n=0, and U=0, so 
that P’=0 and the tangent to the pressure dis- 
tribution in Lagrange space is horizontal at the 
shock front. If the pressure falls off behind the 
front with a tangent of positive slope, the shock 
will decay, and conversely a negative slope 
indicates a shock that will build up (see Fig. 1). 

For a diverging wave, at a given value of 2; 
from the origin or axis, there will be a critical 
value of gradient, depending on the strength of 
the shock and the nature of the fluid, at which 
the free energy influx will balance the divergence 
loss. For gradients smaller than this critical 
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value the wave will decay; for steeper gradients 
it will be intensified (see Fig. 2). 

The nearer the shock front is to its source, the 
steeper will be its critical gradient. The critical 
gradient is obtained by setting U =0 in (7), giving 


nVw,U 
(1+RU?)(21)_ 





, 
e P qh =|- 


These considerations apply equally well to the 
formation of a tail or trailing shock wave such 
as those observed in the case of ballistic N-waves.* 
A pressure gradient in excess of the critical must 
be produced in the wake of the projectile as a 
necessary preliminary to the development and 
growth of the shock. As the critical gradients are 
still steeper in the case of spherically diverging 
waves it appears very unlikely, but not impos- 
sible under unusual conditions, that a blast wave 
will ever develop a tail shock. Also, from this 





point of view, a projectile with a long tapering 
tail will probably not give a ballistic wave with 
a tail shock. 

The relationships between wake gradients and 
time derivatives expressed by (7) can easily be 
specialized for an ideal gas. The Rankine- 
Hugoniot relations relating the shock variables 
can be expressed analytically in this case for 
which the adiabatic equation of state is 
P=A(S)V~, in the well-known form 


P,=C*[2y(U/C)?+(1—y)]/(v¥i+y7)), 

v= 2, U—C?/ U)/(y+1), (8) 
Vi=(C(y—1)+2(C?/U*)]/(y+1), 
R=—V/yP 


(where C is the normal velocity of sound) and 
corresponding time derivatives. This relation 
between the shock variables makes possible the 
explicit expression of P’ and v’ in terms of U 
and U: 


Py . {(1+C2/U2][2yU2+(1—y) C2] +2[ U2(-y—1) +2C2]} 


+1) (U= C4) 


and 


v= 


2i(y+1)*U? 


2U[3+C?/U?][(y—1) U2+2C?7] 2n 


{K(y—1) U?+-2C?)[27U?+(1—y)C*J}, (9) 





(y+1)?U(U?—C?) 


Some applications of these expressions will be 
considered after a brief discussion of the general 
relationships connecting higher derivatives. 


THEORY OF HIGHER DERIVATIVES 


The relationship for second and_ higher 
derivatives may be obtained directly from (6) by 
differentiation and substitution but unfortu- 
nately involves a large number of terms, although 
only in simple algebraic combinations. The dot 
derivatives should be taken first, i.e., P’ and #’, 
since the evaluation of these is required for the 
later calculation of P” and v’’. Recalling that 


2=U, «z=dx/dt+ U(dx/dz) =v+ UV(2/x)", 
R=(0R/dP)P+(aR/dS)S, S’'=S8/U, 
R' =(0R/dP)P’+(S/U)(dR/aS), 
V’ =RP’+(S/U)(aV/aS), 


a DuMond, Cohen, Panofsky, and Deeds, J. Acous. Soc. 
Am. 18, 97 (1946). 
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(L(y—1) U?+2C*][2yU?+(1—7)C*]}. (10) 


2,U*(y+1) 





the process is straightforward and is not shown 
explicitly. The result P,’’ for P” at the shock 
front can be shown to be 


~ Py’ =—(WR+1)?{ (nP’/2,:)(3W’R+1)(V—1) 


—2UR(—(nv/2;)P’+ Uv! — (8°%v/dT?) 

+2U*R[n/2,(2v’ V—v'+vRP’ 

+0($/U)(8V/aS) —(V%/z;)) 

+(UP’— P)((aR/aP)P’ + (S/U)(aR/AS)) ] 

+(U?R—1)[RP'U —RP—(aR/aP)(P— UP")? 

+ (n/2,)(vv' + Vi— UVo' — Vo? /2, 
+vRP—vRUP’)]}, (11) 


where all quantities on the right are evaluated at 
the shock front. In principle then, a knowledge of 
U, U, U,—along with the equation of state and 
derived Hugoniot equation of state—would 
enable a Taylor's series to be constructed for P 
or v as a series in Z and T with coefficients func- 
tions of T. 
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An interesting special case involving higher 
order gradients is that of the constant velocity 
spherically divergent wave found experimentally 
as previously mentioned. In this case U, U, 
vanish and in the wake of the shock S is inde- 
pendent of z as well as ¢. Then 


P,’ =n Vw,U/(2,(1 +Rl ?)) =f U) 21, 

Py" = —(U?R+1){ (nP’/2:)[(Vi-—1)(3U?R+1) 
+v,UR(U?R+ 3) ]+ (no1'/2,:)[U(3ViL?2R 
+ Vi—2U°R)+0,(U?R—1) ] 
+ (nv V/2,?)[ — U?R(2U(*/8Vi) +0:) +01] 
+ U°P,”(dR/dP)(U?R+1)} 

= g(U)/z,?. 

This is of form g(U)/z,?, since P;’ and 2’ are of 

the form f( U)/z;. Similarly, P)’” =h( U) /z;*. Thus 

P=P\+P/(Z4+P,"27/2!+ ere 


= P\+(f(U)/U)(Z/T) + (h(U)/2U*)(Z/T)? 
+:---=P(Z/T). (13) 


(12) 


At the origin, Z=—UT, so that the pressure 
there is always constant and determined only by 
U. From dimensional considerations an explosion 
with constant pressure at the origin must produce 
a distribution which is a function of Z/T.** For 
such an explosion the pressure distribution might 
look like that indicated in Fig. 3 and be expected 
to be found where a high pressure is suddenly 
produced in a large container with a small 
opening so that the source of pressure is main- 
tained for a considerable length of time. This 
situation is in contrast to the opposite extreme 
in which the energy release occurs in a neg- 


shock 
pressure 


ome Py 











Wa Distance 








Fic. 3. Propagation, at constant velocity, of a shock origi- 
3 nating from an explosion at constant pressure. 


** This follows as pressure of dimensions (ML~!7~-*) and 
the undisturbed air density of dimensions (ML~*) are the 
only relevant variables. Their ratio of dimensions (L?7~*) 
is independent of mass and is the only combination of L 
and T determined by the data. 
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ligibly small time interval. In most explosions 
energy is released over a small but finite time 
interval, and this generation of energy is spread 
over a finite region of space. For such a finite 
source of given explosive energy release, the 
more the concentration of the explosion in space 
and time the greater will be the strength of the 
initial shock but the faster will be its initial rate 
of decay. If the pressure is maintained for a 
longer time by slower release of energy, the 
initial shock will be less intense but also will 
decay less rapidly since it is augmented at a 
later time by the influx of free energy. However, 
at long enough times compared with the temporal 
spread of the explosion and at great enough dis- 
tances compared with the region of energy 
release, the shock strength should approach a 
value determinable from the total energy of the 
source and the distance from the origin. 


APPLICATION TO WEAK SHOCKS 


It is of considerable interest to observe the 
forms that the wake pressure derivatives take 
when the Mach number approaches unity. 
Equation (9) reduces to the form 


Py’ = —2(C*/y+1)[(4U/U?-C)+n/21] 
= —2(C/y+1)[ (261/21) + (nC/21) J. 


Similarly, Eq. (12) reduces to 
Py" =(4C/y+1)[3(0*v/dT?),/0°+462/0;]. (15) 


(14) 


In the plane case, n=0, the differential equa- 
tions (4) can be solved for weak waves where the 
entropy is practically constant. For this case 


dv/dt= —OP/dz; dv/dz=R(AP/dt). (16) 
Eliminating v through differentiation gives 
— (#P/d2*) = R(e’P/dt*) + (dR/dP)(dP/dt)?. (17) 
The equation of state reduces to 
V=[Py/C}", 


where A(S) is considered constant. Also, 


R=0V/dP=—V/yP; (18) 
whence, 
| aR/aP =(1+y/y*)(V/P2). 
Then (17) becomes 
—yA-"NyP + 7192P /dz* 
=0P/d?+(1+7/y7)-(1/P)(aP/dt)?. (19) 
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A time-dependent solution of this equation is 


P = p/po= (po/po)(2/ct)?/*7, (20) 
giving 
v= (2¢/y—1)[(b/po) 7-1], = (21) 
and 
Op /dz=(2ypo/1+y7)(1/ct)(2/ct)O-P 7, (22) 


if p=po(1+6), where 6, the over-pressure in 
atmospheres, is small compared with unity, then 


v/e=(2/(y—1) 01+ ((y—1)/2y)6+----1] 


=6/y. (23) 


But this is also the approximate Hugoniot con- 
dition between particle velocity and over- 
pressure for very weak shocks. Therefore, the 
abolve solution (20) for pressure can be cut off 
by a shock boundary at any value of z where 6 
is small enough. As a case in point, Fig. 4 depicts 
how this time-dependent solution to the hydro- 
dynamic equations can be terminated by head 
and tail shocks, H and 7, to represent a plane 
N-wave. The pressure distribution is not a 
straight line between the terminating shocks. 

When z=ct, p= pf» so that the plane of pressure 
Po travels with constant velocity c in Lagrange 
space and therefore also in actual space for weak 
disturbances. 

The pressure gradient at the terminating 
shocks will differ only very slightly from the 
slope of the pressure at p=» so that the rate of 
decay of the shock terminations can be closely 
determined by putting this value of gradient in 
(14), giving 

(2y/1+-)(popo/Ct) =(—4C/y+1)(b1/01), (24) 
whence 

b,/v3= —1/2t or 


61 ‘piC=v,=kt, (25) 
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Fic. 4. Progress in time for a plane N-wave constructed 
from a time-dependent solution of the hydrodynamic 
equations terminated by weak shocks. 


where & is a constant. Therefore, to this degree 
of approximation both the head and tail shocks 
should decay inversely as the square root of the 
time. To do this they must move away from the 
plane z=Ct, at which p= po, as time progresses. 
Assuming the length \ of the N-wave is \=2; 
— Ct, and recalling that 


(1+6)(1+y)/2y=1+(1+7/y7) 61 


=2,/Ct=(A+Ct)/Ct=1+A/Ct, (26) 
it follows that 
A= (1+y/y)Ct6. 
Therefore, from (25) 
A=(1+y/y)Ct-kpoCt, ~t. (27) 


The length of the plane N-wave increases as the 
square root of the time. It should be noted that 
these laws of decay of weak plane waves are 
independent of the y of the fluid. 

Because of the explicit occurrence of x in Eqs. 
(1) for cylindrically and spherically divergent 
waves, the method developed here for the study 
of the decay of plane waves is not directly 
applicable. Further studies are being made to 
determine whether as rigorous a method can be 
developed for diverging waves. 
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Surface Currents on a Conducting Sphere Excited by a Dipole* 


C. H. Papas AND RONOLD KING 
Cruft Laboratory, Harvard University, Cambridge, Massachusetts 


(Received January 23, 1948) 


This paper treats the problem of determining the current distribution on the surface of a 
perfectly conducting sphere when driven by a dipole antenna erected on its surface. Curves 
of the real and imaginary parts of the surface currents are given for the case of a half-wave 


dipole and various radii of the sphere. 


I. INTRODUCTION 


HE distribution of currents on the surface 

of a metal spheroid excited by an adjacent 
antenna or by a cavity is of great practical 
interest (see Fig. 1). The analytical problem 
involved is closely related to the determination 
of the diffracted field of a plane wave incident 
upon a sphere or spheroid,' which has been 
studied by many investigators. The rigorous 
solution of the general case in which the wave- 
length of the incident plane wave is of the same 
order of magnitude as the dimension of the 
diffracting sphere is due to L. Brillouin. In 
practice a plane wave is approximated by the 
field of a radiator situated a great many wave- 
lengths from the diffracting sphere or prolate 
spheroid. This field induces surface currents 
which, in turn, set up an electromagnetic field. 
This induced field is so disposed that on the 
surface of the sphere or spheroid (assumed 
perfectly conducting) the sum of the tangential 
components of the incident and induced E-fields 
vanishes and the sum of the tangential compo- 
nents of the H-fields equals the surface current 
density. Since the scalar wave equation is 
separable in spherical (and spheroidal) coordi- 
nates, these boundary conditions permit the 
determination of the unknown constants of the 
infinite series solution for the induced E- and 


’ H-fields. 


If the distant radiator is brought close to the 
diffracting body, the incident field is no longer a 
plane wave and the analytical problem becomes 


*The research reported in this document was made 
possible through support extended Cruft Laboratory, Har- 
vard University, jointly by the Navy Department (Office 
of Naval Research) and the Signal Corps, U. S. Army 
under Contract N5 ori-76, T. O. I. 

1J. A. Stratton, Electromagnetic Theory (McGraw-Hill 
Book Company, Inc., New York, 1941). 
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so complicated that the general procedure de- 
scribed above has to be abandoned. An alter- 
native procedure has been formulated by Feld? 
who has applied Lorentz’ lemma to a linear 
antenna over a perfectly conducting sphere. The 
current in the antenna is assumed to be sinus- 
oidal; the coefficients of the modal currents on 
the sphere are determined by the well-known 
Fourier-Lamé method. 

By applying the general reciprocity theorem 
(Appendix A) to a closed metallic surface with a 
linear antenna erected upon it, an integral equa- 
tion is obtained. In order to solve this integral 
equation by the Fourier-Lamé method it is 
necessary to restrict the closed metallic surface 
to a shape that may be described simply using 
coordinates in which the homogeneous scalar 
wave equation, 


Vy +k*y =0, 


is separable.‘ Another limitation (which is of 
purely practical importance) is the availability 
of numerical tables of the functions used in the 
solutions. For example, if the closed metallic 
surface is a sphere, the ‘“‘natural’’ coordinates are 
spherical which permit separation of variables 
and give as angular solutions the Legendre 
polynomials.and as radial solutions the spherical 
Hankel functions. Both the Legendre poly- 
nomials® and spherical Hankel functions* have 
been tabulated. Similarly, if the metallic surface 
is a prolate spheroid, the prolate spheroidal 
coordinates! £, », @ are used. These also permit 


2 J. N. Feld, Doklady U.S.S.R. 51, No. 3, 203-206 (1946). 
eae) Grunberg, J. Phys. U.S.S.R. 10, No. 4, 301-320 
1946). 

4 L. P. Eisenhart, Annals Math. 35, No. 2, 284-305 (1934). 

5 Tables of Associated Legendre Functions (Columbia 
University Press, New York, 1945). 

6 Tables of Spherical Bessel *Functions (Columbia Uni- 
versity Press, New York, 1945). 
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separation of variables and yield the functions 
Se(n) as “angular” solutions and Re;() as 
“radial’’ solutions. These functions have been 
tabulated.’ It is evidently possible to solve for 
the surface currents on spheres, prolate sphe- 
roids, and oblate spheroids, and from these to 
determine the transmitting and receiving char- 
acteristics of the array consisting of linear 
antenna and closed surface. 

Although the spheroidal functions have been 
tabulated, the additional computation necessary 
to determine the surface currents is considerable. 
Therefore, the simpler problem of a linear 


antenna over a sphere has been carried through 
first. 


Il. LINEAR ANTENNA OVER SPHERE 


For an array consisting of a linear antenna of 
height h over a perfectly conducting sphere of 
radius a, the total current crossing a parallel of 
latitude @ on the sphere is (Appendix C) 


I,s(0) = 3 B,.P,,(cos6), 


(1) 
where 
B,= (Imax/2)[Lon—1(kd)/pn—1' (ka) ] 
—coskh{_p,»—1(ka) /pn—r' (ka) ] 
ay [pn+i(kd) /pn+1'(ka) | 
-+coskh[_pn+i(ka)/pny1'(Rka) j], (2) 


and k=2z/x, a is the radius of sphere, h is the 
height of antenna, d=a+h. It was assumed in 
the derivation of (1) that the current along the 
antenna has the form 


T,a(r) =Imax sink(d—r), 


and that this distribution is independent of 
currents on the sphere. For very thin antennas 
driven at the junction with the sphere a sinus- 
oidal distribution is a good approximation for 
lengths differing considerably from \/2. In gen- 
eral, a sinusoidal current cannot be maintained 
by a single generator but requires a distribution 
of generators along the entire antenna. 

When the length of the antenna is \/4, the 
current is maximum, Imax, at the base (r=a) 
and zero at the top (r=a+h). Since kh=7x/2, 


7Stratton, Morse, Chu, and Hutner, Elliptic Cylinder 
and Spheroidal Wave Functions (Massachusetts Institute of 
Technology Press, Cambridge, 1941). 
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which may be solved for surface currents by the method 
of this paper. 


the expression for B, reduces to 


B= (Imax/2)[Lon—1(kd) /pn—s' (ka) } 
i (en+1(kd) /Pn+1' (ka) J}. (3) 
Calculations have been made for various 
values of the radius a with kh=2/2 using (3). 
In order to compute terms such as the ratio 
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Pn—1(kd)/pn_i'(ka) it was found convenient to 
follow a computational procedure used by L. 
Infeld.* By factoring this ratio it is seen that 
one of the factors is a function which satisfies 
Riccati’s equation and for which there is a 
recurrence relation. That is, 


pn—1(Rd)/pn_1' (ka) = pn—1(ka)/pn—s' (ka) 
‘pn—1(ka+kh) /pn—1(ka) - €n-1(Ra) 
pn—i(ka+kh)/pni(ka), (4) 


where 
Fn—1(R@) = pn_i(ka)/pn_i'(ka). 
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This zeta satisfies Riccati’s equation and equals 
the zeta used by Infeld with the equal to his 
2m+1. Letting x =ka, the recurrence formula is 
On(x) = [nxt n_1(x) — x? ] 

[(x*@—m?)fn-1(x)+nx], (5) 
and ¢_,;(x)=-++7. Starting with ¢_;(x)=-++17, the 
entire set of functions fo(x), €:(x), --- was 
obtained. 

The other factor which is of the form 
Pn—i(Ra+kh)/p,_1(ka) was computed from the 
fundamental definition of the functions: 


Pn-i(Ra+kh) (x+kh)[(9/2(x+kh)) Tins (x +kh) + (i/cosnm) (9 /2(x+kh))! Toray (x +h) 





pn—1(ka) 


The functions 
(w/2x)*Tngy(x) and = (#/2x)!J_¢n44)(x) 


are the spherical Bessel functions, tables of which 
are available. 


§L. Infeld, Q. App. Math. 5, No. 2, 113-132 (1947) 
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xf (w/2x)*Tna3(x) +i(1 /cosne) (e/2x)8IT—cnay(x)] 





Using (6) and (5), the ratios p»_:(ka+kh) / 
pn-1'(ka) were computed for n=0, 1, 2, 3, 4---11 
when kh=x/2 and ka=7/8, 4/4, 2/2, 32/4, z=, 
52/4, 3x/2, 77/4, and 2x. These ratios were 
substituted into (3). The values of B, thus 
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obtained were substituted into (2) and J;s(@) 
was determined. 


Ill. RESULTS 


The real and imaginary parts of the current 
I,s(@) crossing the @ parallel of latitude have 
been computed for several values of the radius a 
of the perfectly conducting sphere in the special 
case of an antenna of length h=)/4 erected on 
the sphere. Plots of these currents are shown in 
Figs. 2-10. The ordinate of these plots is twice 
the ratio of J(@) to Imax. Since the positive 
direction of current along the antenna is towards 
its top and the positive direction of flow on the 
sphere is from its North pole towards its South, 
it is necessary for J,s(0) to be opposite in sign 
to Ia(a)=Imax. In all the plots for 6=0 the 


imaginary part is zero and the real part equals 
—2. That is, 


211s(0)/Imax = —2 or I,s(0) seat — Imax. 


At the bottom of the sphere, @= 180°, the current 
is zero. 
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In the theoretical derivation it was assumed 
that there was a generator or a distribution of 
generators somewhere in the antenna. The only 
restriction upon this generator is that it be 
independent of the angle ¢. For example, the 
generator may be a single idealized slice generator 
or a distribution of these along the antenna. In 
the computed case the current distribution is 
assumed to be sinusoidal with a node at the top 
of the antenna and a loop Imax at its base. This 
requires a single slice generator to be placed at 
the base of a very thin antenna as indicated in 
Fig. 11. The current distributions shown in the 
plots are for just such a model. 

Since an actual physical antenna of finite 
cross section driven from a coaxial line over an 
imperfectly conducting sphere necessarily differs 
from the ideal model assumed in the theory, 
actual currents must be expected to differ some- 
what from the computed currents represented in 
the plots. However, it is safe to assume that the 
currents in the idealized mathematical model are 
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Fic, 6. Total current on sphere versus polar angle for 
ka=z. 
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Fic. 7. Total current on sphere versus polar angle for 
ka = 5/4. 


a fair approximation of actual currents in a 
suitably arranged physical set-up. The determi- 
nation of such actual distributions using tech- 
niques developed in this laboratory for the 
measurement of magnitude and phase of surface 
currents is contemplated. 
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APPENDIX A 
The General Reciprocity Theorem 


A given source or distribution of sources 
maintains a current density 7; with an associated 
electromagnetic field E,, H;..A second source 
maintains a current density iz with an associated 
electromagnetic field E2, Hz. Then, according to 
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the theorem of Lorentz, 


V-(E,X He) —V-(E2X Mi) =i: E2—te: Ey. (A1) 


Up to this point there is no relation either 
between the two sets of fields and currents or the 
regions in which they are defined. Let it be 
assumed that the two sets of fields are defined in 
tree*space and are bounded by the same geo- 
metrical surfaces. The boundary conditions on 
these surfaces need not be related in any other 
manner than that they are applied at surfaces 
that satisfy the same mathematical equations. 
Let (A1) be multiplied on each side by an ele- 
ment of volume dr and integrated over a volume 
V in empty space that is bounded by a closed 
surface S and extends to infinity where all fields 
vanish. Since all volume densities of current 
vanish in empty space, the right side of (A1) 
when so integrated vanishes. The left side may 
be transformed using the divergence theorem, 
since the fields are continuous in the volume V. 
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The result is 
[tv Gxt -E.xHy ide 
§ 
Ss 


where 7 is an external normal to the volume V 
on the bounding surface S. The only restrictions 
on the two sets of fields is that they satisfy the 
field equations, that the volume of integration 
V be empty space, and that the bounding surface 
S have the same geometrical shape for both sets 
but not necessarily the same physical properties. 

For. periodic time dependence (A2) is valid 
for complex amplitudes provided both fields 
vary with the same frequency. 


APPENDIX B 
Derivation of the Integral Equation 


The problem under consideration is to deter- 
mine the distribution of current on the surface 
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of a perfectly conducting sphere of radius a when 
excited by a given distribution of current along 
a linear antenna of very small cross-sectional 
area. The antenna is placed along the line 6=0 
in a system of spherical coordinates r, 0, @ with 
origin at the center of the sphere as shown in 
Fig. 11. Let the volume V be the entire volume 
in empty space outside the sphere and the 
antenna. Similarly, let the inner bounding surface 
S consist of the surface Ss of the sphere and the 
surface S, of the antenna. The outer boundary 
is at infinity where all fields vanish. 

The ‘electromagnetic field vectors EZ; and H, 
are defined to be caused by the actual currents 
in the antenna and on the surface of the sphere. 
On the perfectly conducting sphere they must 
satisfy the boundary conditions: 


nXE,=00n Ss, (B1) 
naXH,= Kis on Ss, (B2) 


where # is an external normal to the sphere and 
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hence an internal normal to the volume V 
bounded by S, and where Kis is the surface 
density of current on the sphere. The current 
i14 in the antenna and Kis on the sphere are 
maintained by an appropriate generator or 
distribution of generators in the antenna so dis- 
posed that no currents exist in the ¢-direction 
(i.e., around an axis through the antenna), in 
the antenna or on the sphere. Hence, 


Kis= 6K 1s(8); 
(B3) 


f K,s(@)-a sinéd¢ = I1s(8), 


where J,s(@) is the total current crossing a 
parallel of latitude on the sphere. Also the total 
current in the antenna is axial, i.e., radial in the 
system of spherical coordinates. Evidently, com- 
plete rotational symmetry is obtained about the 
axis @=0 through the antenna and the sphere. 
The line integral of the tangential component of 
the magnetic field H, around the surface of the 
antenna is equal to the total axial current J,4(r) 
in the antenna. Thus, with Rp the radius of the 
antenna is 


ptvas= f H,Rodo 


=f isalnas= tar), (B4) 


cross 
section 


The electromagnetic vectors E, and Hz are 
not actual fields related to the real currents of 
density 7; in the antenna or the sphere. They 
constitute an auxiliary, fictitious field that is 
introduced to facilitate the derivation and solu- 
tion of the integral equation. It has been pointed 
out that E, and H, must satisfy the field equa- 
tions and that they must be defined on the same 


. geometrical surfaces, i.e., on the surfaces of the 


sphere and the antenna, as E, and A). Similarly, 
iz and E» are defined in the same volume as i; 
and £,. Let it be assumed that the volume V, 
(occupied by the antenna for E,, H;) is empty 
space for Es, Hz. Also, let the volume Vz 
(occupied by the perfectly conducting sphere for 
E,, H,) be constituted to contain appropriate 
generators or distributions of generators to 
maintain the following prescribed field on the 
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surface Ss of the sphere 


Ex%=0; Ex=U sinéP,(cosé), 
Hou=0; H.,=0. 


(BS) 


Thus the entire field in empty space and its 
boundary Sy on the sphere has the components 
E>,, Exe, Ho, which must satisfy the field equa- 
tions and the prescribed boundary conditions 
(BS). 

After rearrangement of the vector products in 
(A2) they may be written as follows: 


Ss 


+f (Hs: (XE1) —Hy-(@XEs)\do=0. (BO) 
SA 


Let the integrals in (B6) be considered in turn. 
As a result of (B1) and (B2) the first integral 
becomes 


SS 


-f E.-Kisdo= E»K,s(0)do 
SS “SS 


-{ Eel, s(0) ade. (B7) 
0 


The second and third steps follow from (B3) 
and (BS). 

On the surface S, of the cylindrical antenna 
both (7X E,) and (7X2) are constant around 
each cross section because of rotational sym- 
metry. Furthermore, since E,=0, 


nXE=—-O0XE=—(6X#)E, — ; 
—(@XO0)E,=E,. (B8) 
Hence, since H = $H, for both fields, 


f {| H2-(nXE:) —H,- ("i XE») }do 


SA 
r=a-+h o=2 
” J E,dr f Hx,Redo 


r=a o=0 
r=ath o=2r 

- f End f HRudg, (B9) 
r=a o=0 


where Ro is the small radius of the antenna. 
However, from (B4) the contour integral of the 
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magnetic vector H,; around the periphery of the 
antenna equals the total current in the antenna. 
Accordingly, 


a+h 2” 
[Badr [tigate 
ra o=0 


ath 
-{ Es,I,a(r)dr, (B10) 


where J,4(r) is the total current in the antenna. 
Since the surface S4 (which is the surface of the 
antenna for E,H;) encloses only empty space 
for E2f2, the same contour integral with H, 
instead of H,; equals only the time rate of change 
of e9E2, integrated across the small cross-sectional 
area rR,?. Since Ry can be made as small as 
required, this quantity can be made negligible 
so that 


r=at+h o=2r 
f Eds f H2Roddr =0. (B11) 


r=a o=0 


Substitution of (B7), (B10) and (B11) in (B6) 


gives 


7 r=at+h 
f Exlss(0)ad0= f Eo,I14(r)dr. (B12) 


=a 


This is the desired integral equation for deter- 
mining J,s5(@). I:a(r) is assumed given, Eo¢ is 
the prescribed «field given in (B5), and Eo, is 
determined from the vector Ee, which satisfies 
the field equations and the boundary conditions 


(BS). 
APPENDIX C 
Linear Antenna and Sphere 


The solution of the integral equation (B12) 
for the current J,s(@) on a metal sphere on 
which is erected a linear radiator with current 
Iia(r) is carried out as follows. Due to the 
orthogonality of the Legendre polynomials, that 
is, 


1 
J PaW)Paw)du=[2/2n+1) Ban (C1) 
—1 
whee iawt 2" tite ah 
ere Omn = 0 if mxn' e oneckKe ’ 


is possible to expand a function in terms of 
Legendre polynomials. Without limiting the gen- 
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Fic. 11. The geometry of sphere excited by dipole. 


erality the current J,s(@) on the surface of the 
sphere can be represented as follows 


I,s(0) => B,P (cos). (C2) 
n=() 
B, are independent of @ and later it will be 
shown that they are functions of the radius a of 
the sphere and the height A# of the antenna. 
Physically (C1) means that the total current 
crossing the parallel of latitude @ is the sum of 
all the so-called modal currents crossing this 
parallel of latitude; m is the order of the mode. 
As n increases, the amplitude of the current in 
this mode decreases so that it is possible to 
obtain a rather precise evaluation of J,s(@) by 
summing over the first twelve modes. Mathe- 
matically this means that (C1) converges for 
0=6=r rapidly enough so that terms of higher 
order than twelve contribute very little to the 
sum. 
Expressed in appropriate coordinate form 
(B12) becomes the following for r=a: 


at+h 
f Fo,Iya(r)dr 


Ls 


= Ew > B,»P,(cosé)-ad@. (C3) 


0 n=0 


Furthermore, as was previously stated, E2,, Exe, 
and H2, are the only non-zero components of the 
field due to currents on the sphere; this means 
that the field is transverse magnetic TM. 
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It has been shown! that a TM field is uniquely 
derivable from a function u. Moreover, u satisfies 
the scalar wave equation : 


V2(u/r)+k(u/r) =0, (C4) 
where k=w/c. The general solution of (C4) is 


u= > AynP»(cosd) pm(kr), 


m=) 


(CS) 


where P,,(cos@) is the Legendre polynomial of 
order m, and 


pu(x) = (4x /2)'Hm4° (x) 


is the weighted spherical Hankel function of the 
second kind. The components of the field are 
obtained from 


E2,=([k?+(d2/dr?) Ju, | 
Ex9=(1/r)(0?/drd8)u, | 
Hog = —1(k/r)(0u/08). | 


(C6), 


(C6) 


From (C5) and 


follows: 


Eo, is evaluated as 


Ew = —(k/a) > AmPm'(COSO) pm’ (Ra), (C7) 
m=( 
where P,,,'(cos@) is the associated Legendre func- 
tion. In order to evaluate the A,, in (C7), (B5) 
may be written in the form 


Eo= U sinéP,(cosé) 
= U[[ Pn+:'(cos@) — P,_1'(cos0) ]/(2n+-1)]. 
(C8) 
Equating the coefficients of P,,'(cos@) in (C7) 
and (C8) leads to 
A nui= —aU/(2n+1)-k-pns1' (ka), (C9) 
An-1= +aU/(2n+1)-k-pn_i'(ka). (C10) 


Substitution of these values of A,, in the first of 


_ Eqs. (C6) results in the following expression for 
Eo: 





0? WaP »-1(cos@) pr_i(kr) 
motes 2] 


art Ik (2n-+1)-k- pa1/(ka) 
OP n+1(COS8) pnyi(Rr) 
* (2n+1)-R- pass’ (Ra) 


*P. Frank and R. v. Mises, Diff. Gleichungen d. Physik, 
(Mary S. Rosenberg, New York) p. 872. 





| (C11) 
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Now let it be assumed that the antenna current is 
Tia(r) a sink(d—r). (C12) 


With (C12) the left side of (C3) may be com- 
puted using (C11): 


oth Ual max 
f Esla(r)dr = 
e (2n+1)-kR- pn_1' (ka) 


° 


ath Ps be 
xf sink(d—1)-( t+ —) pasar) 
or* 


i) oth 
—_—— ———; f sink(d—r) 
(2n+1)-R-+ pri’ (ka) Y, 


oe? 
. (#+—)on, i(Rr)dr. 
or* 











This integration is carried out along the antenna, 
i.e., 6=0, and since P,(1)=1, the Legendre 
polynomials appearing in (C11) equal unity. 
Explicitly, the integration yields 














on | 
J Eola (r)dr =——_——_- 
a 2n+1 
Pn—1(kd) Pn—1(ka) 
x ——coskh —— 
Pn—1 (ka) Pn—1 (ka) 
pnss(kd) puii(ka)] 
_ ——-+-coskh ——j. (C13) 
Pn+i (ka) Pn+i (ka 


With (C8) the right-hand side of (C3) be- 


comes : 
Ew > B,P»(cosé) -adé 
0 n=0 
=B,-[ U2a/(2n+1)]. (C14) 


Equating (C13) and (C14), the following expres- 
sion for B, is obtained: 


B _Fmasf Pn—i(kd) a 
2 Lpn1’(ka) 

J Pn+i(kd) 

Pn+1 (Ra) 


By substituting B, from (C15) into (C2), Ji s(6) 
is determined. 


Pn—i(ka) 
Pn—1 (ka) 


Pn+i(ka) 


coskh 





+cos 


| (C15) 


Pn+ 1‘ (ka) 
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Spurious Signals Caused by Noise in Triggered Circuits* 


Davip MIDDLETON 
Cruft Laboratory, Harvard University, Cambridge, Massachusetts 


(Received January 23, 1948) 


The expected number of peaks/sec. (m¢4n), caused by 
noise or by a signal and noise, in excess of a given amplitude, 
or triggering level, is determined for broad-band noise (a) 
and for the envelope of a narrow-band disturbance (b). 
The signals considered here are pulses of various shapes 
(rectangular, triangular, trapezoidal), or a carrier modu- 
lated by these pulses. The expected number of peaks per 
second is directly proportional to the band width of the 
noise, and depends in some degree (20-30 percent) on filter 
shape, the gaussian response yielding a larger value of 
Non than the rectangular, for the same total input power. 
Further, 174" is greatest when the triggering level is the 
most probable amplitude. Pulse shape is not significant 
when the pulses are on only a small fraction of the total 
period of observation, but becomes so if (1) the time of 
observation is comparable to the pulse length, and (2) the 
pulses have sharp leading edges. The latter gives rise to a 
transient effect which is not ignorable. This transient 
contribution does not appear in the envelope of narrow- 
band waves, as long as the maximum time-rate of change 
of the wave form of the envelope is much less than the 
central or resonant frequency of the JF, i.e., as long as the 


concept of “envelope’”’ is meaningful. Rectification does 
not change ,,n, provided that there is no appreciable 
distortion of the envelope by video- or audio-filters and 
that the triggering level of the input is passed by the 
detector. 

Improvement, defined in ways appropriate to the sys- 
tems in question, of one circuit which can tolerate A; 
spurious signals/sec. in excess of a given triggering level 
over another circuit, which can tolerate only Az sueh 
signals (A,> Az), is calculated for broad-band signal and 
noise, and for narrow-band waves following half-wave 
linear or quadratic rectification. In general, the improve- 
ment is slight (a few db) if the trigger levels of the two 
circuits are not close to the most probable value of the 
amplitude and if the ratio of the number of spurious signals 
allowed/sec. to the band width of the disturbance is small. 
Significant improvement follows only when the maximum 
tolerable number of triggering peaks/sec. is the same order 
of magnitude as the band width of the noise. Tables and 
figures illustrating the above are included in Part I; the 
mathematical derivations are contained in Part IT. 





I. DISCUSSION 
1. Introduction 


N all instances in which signals of a certain 

magnitude are used to trigger a circuit, noise 
accompanies the triggering signal and may give 
rise to spurious responses as a result of the 
presence of noise peaks in excess of the excitation 
level, particularly if the signal voltage is com- 
parable to the r.m.s. noise amplitude. Thus it is 
important to determine how many noise peaks 
per second, or how many on the average in a 
given time interval, will exceed the threshold 
level. Further, we may ask to know the improve- 
ment in using a circuit which can tolerate A, 
spurious pulses per second over another circuit 
which can tolerate only Az (<Aj,) spurious 
pulses per second, the various filter widths and 
detector characteristics, etc. being essentially 
the same. In other words: how much greater is 


* The research reported in this document was made 
possible through support extended Cruft Laboratory, 
Harvard University, jointly by the Navy Department 
(Office of Naval Research) and the Signal Corps, U. S. 
Army, under Contract N5ori-76, T.O. I. 


VOLUME 19, SEPTEMBER, 1948 


the sensitivity of circuit No. 1 over No. 2, given 
A, and A:? This type of consideration and 
question clearly applies to all systems that de- 
pend for their performance on the knowledge of 
the number or average number of pulses received 
in a given time interval. Of course, the signals 
need not be in the form of pulses; for example, 
they may be sinusoidal, but the former are most 
widely used for trigger operations. 

The analytical treatment of our problem is 
based on the observation that the expected 
number of peaks in excess of a certain value Vo 
in the time interval (to, fo+¢) is the same as the 
expected number of crossings of the level Vo with 
positive (or with negative) slope in the range 
(to, to+z¢), since once this level is exceeded, we 
have a triggering pulse. It is also clear that this 
number of crossings is not the same as the 
number of maxima, since the minimum between 
two successive peaks may lie above the level Vo. 
However, when Vp is large, it may be shown! 


1S. O. Rice, ‘Mathematical analysis of random noise,” 
Bell Syst. Tech. J. 24, 46 (1945). See Section 3.6, Eq. 
(3.6-11) in particular. 
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that we do approach the condition that the 
number of peaks or crossings is, on the average, 
equal to the number of maxima. 

Rice? and Kac* have obtained, independently, 
expressions for the expected number of zeros 
(Ve=0) of a random wave. Our results (see 
Part Il) are a generalization to the cases where 
Vo#0 and where there is a signal along with the 
random disturbance (Sec. 1, Part Il). The ex- 
pected number of peaks in excess of Vo, in the 
interval to, t9+¢ follows from (II-1.5), and is 


tort 


Non Vo, t) -{ arf V W, ( Vo, V: t)d V 
to 0 


-{ p(Vo; edt, (I-1.1) 


ut) 
while the expected number of peaks per second is 


No+n( Vo, t) = No+n( Vo, t) t 


1 


to+t 
t Ht 


The subscript o refers to a possible signal and N 
to the noise. Here W,(Vo, V; £) is the probability 
density for the variables V(=V») and V, and 
b( Vo;t) is the probability density for the number 
of peaks exceeding Vo, both of which have been 
evaluated in Part I1 for a number of important 
cases. When the random process is stationary, 
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Fic. 1. Typical spectra of broad-band (a) and narrow-band 
(b) signal and noise waves. 


2S. O. Rice, Am. J. Math. 61, 409 (1939), and see also 
Eq. (3.3-5), reference 1. 


3M. Kac, Bull. Am. Math. Soc. 49, 314 (1943). See refer- 
ence 8. 
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i.e., the mechanism underlying the process does 
not change with time, we see at once that 


No+n(Vo, t) =t- p( Vo) 
and 


Ne+Nn( Vo) = p(Vo), (I-1.3) 


from (1I-1.1) and (I-1.2). 

Before turning to a consideration of specific 
results we note that there are two types of noise 
or noise and signal that are of interest here. 
These are (a) broad-band noise, the width Af of 
whose spectrum between half-power points is 
equal to or greater than fhe resonant frequency, 
and (b) narrow-band noise, whose spectral width 
Af is much less than the central or resonant 
frequency fo.‘ In either case, an accompanying 
signal is assumed to be spectrally of the same 
order of magnitude as the noise. These remarks 
are illustrated in Fig. la, b; here Wo» denotes a 
spectral intensity (units: power/frequency). 
Type (a) spectra occur more frequently in 
acoustical problems, while type (b) is of greater 
interest in the reception of radio or radar signals. 

When the disturbance consists of narrow-band 
noise and signal, for example a pulse-modulated 
carrier, it is usually rectified first, and the low 
frequency part® is then applied to trigger the 
circuit in question. Rectification does not destroy 
the validity of our results as long as the low- 
frequency filter following the detector is suffici- 
ently wide to prevent any appreciable distortion 
of the detected envelope because of frequency 
selection, and provided, of course, that the recti- 
fier passes voltages in excess of the triggering 
level. This is not surprising, since detection under 
these conditions merely distorts the wave with- 
out changing the number of peaks greater than 
the trigger level; relative wave shape is pre- 
served, although the detector characteristic may 
be quite non-linear. Equations (1-1.1)—(1-1.3) may 
still be used, except that Vo is now replaced by 
Ro, the triggering level in terms of the undetected 
envelope of the incoming signal and noise, and 
V by R, the time rate of change of this envelope. 
The formal justification of these modifications 


* We consider here symmetrical spectra only. 

5See D. Middleton, “Some general results in the 
theory of noise through nonlinear devices,’’ Quart. 
App. Math. 5, 445, Section 3, 1948, for the case of a 
signal as well as noise. 
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TABLE | 
Type of noise spectrum Analytic form bo be p =(b2/bo)t Af vs wF 
1. Rectangular 1 
‘ _9,f>fri Wwr Wwe "ay ori 
(a) broad band: w(f)= Wo, 0<f<friJ an a wri/(3) fi=>- 
, _ 0, f>fotfra; f et Wwre W wre , wre 
(b) narrow band: w(f)= We, fo-fra<f <Sotfrs —— ~— wre /(3)4 it 
2. Gaussian 
’ Rist , 
(a) broad band: w(f)= Woes?!" Wer/(e ln (QQ af 
8(7)} 4a 
aii , 
(b) narrow band: w(f) = Woe~@-#0)?/#*F2 W wr2/2(r)! ee wr2/(2)8 Afz=s2? 
TT 
3. Optical 
(a) broad band: w(f)= Wo/(1+2?/w*r:) Wor 2% x fat 
, w—wo\? Wwr wre 
(b) narrow band: w(f)= UW / (1+(—)) = x x Afp=—> 
. wre 2 2 
| is given in Part II, Section 2, and consequently by 
(I-1.1)-(I-1.3) remain unchanged.*® : 
. . P f (eo — wo)*2v(w) deo 
l - 2 0 
: 2. Results and Conclusions p=—= _ : 
. by - 
’ Let us consider the broad-band case (a) first, f w(w)dw 
1 and let us restrict the discussion to include, at 0 
1 most, rectangular pulses of amplitude V,=Ao w=2rf, wo=2mfo. (1-2.2) 
wy vaagaseed a — the cen — When the signal is on or off 
of the vertical leading edge can be ignored. . ‘ 
; . : : ; Vp= _ ; =V,/ } -2. 
y Expressions for a general signal are given in ~° (Vo Ve)/(2b0)', or vo=Vo/(2bo)!, (l 4-5) 
4 Part II, cf. (II-3.5). For the intervals that the espectively. The presence of the d.c. V. simply 
4 signal is on or off the process is stationary, and Serves to boost’ the noise, so that the effective 
. from (I-1.3) and (II-3.14) we have at once for ‘Triggering level is now Vo—V.. The quantity bo 
y the expected number of peaks per second in '§ the mean noise power passed by the filter, and 
i- cmenat Fe: w(f) is the mean noise power spectrum; wo is 
0: : 
g zero for the broad-band disturbances treated 
r : _ ; here. 
1- mo)+n (Vo) = p ty * p=P(Vo), (1-2-1) The general result for the expected number of 
n peaks in excess of Ro in the case of the envelope 
a where p is the band-width factor, proportional to narrow-band noise and/or a signal has been 
y the width Af of the noise spectrum and is given obtained soll Part II; see Sections 4 or 5, in 
Ly eT es particular Eqs. (II-4.9) and/or (Eq. II-5.11). 
ry * Care must be exercised in dealing with signals that Where there is no signal we have 
od have a large time rate of change, of which the idealized, ‘ 
rectangular pulse used here is an extreme example. Signals a 
id with steep leading edges contain transients whose contri- nn(Ro) =p{ ——e-™™ }, ro=Ro/(2bo)!, (1-2.4) 
e. butions to the expected number of peaks exceeding Vo (x)! 
often cannot be ignored. When we deal only with the : 
ns envelope these transients do not appear, to a first approxi- and for an unmodulated carrier 
mation, since they consist chiefly of frequencies outside 
the spectral limits of the envelope. A rigorous treatment ro 
ae would include the effects of finite filter widths, but this is Nea. (Ro) = p| ——e~"®?—#°"Io(2aoro) |, 
rt a refinement unnecessary for the purposes of the present (3)! 
“ work, since pulse shape in most instances does not have a 
critical effect on the number of peaks exceeding Ro. y= Ao/(2bo)', (1-2.5) 
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Fic. 2. Expected number of noise peaks per second in 
excess Of Vo for broad-band noise. 


where Ay is the peak amplitude of the carrier. 
Equation (1-2.5) applies equally well for a 
carrier modulated by a rectangular pulse® for the 
interval the pulse is on, while (I-2.4) holds 
during the time the pulse is off. 

Although the rectangular pulse is an idealiza- 
tion, other pulse shapes, such as the gaussian, 
triangular, trapezoidal, etc., give values of m.i~ 
that are not much different from those obtained 
in the rectangular case. Only for strong signals 
(relative to the noise) is the number of peaks in 
the rectangular case greater than for the other 
pulses (of the same maximum amplitude and 
width between half-power points), since in the 
former the noise is “‘boosted’’ during the entire 
duration of the pulse, while for the others this is 
not true. As mentioned in Section I, pulse shape 
is not critical in its effect on m,,», and is of 
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Fic. 3. The expected number of noise peaks per second 
in excess of ro (=triggering level/rms noise level). The 
d.c. component is included. 
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negligible importance when the duty-cycle 7 is a 
small fraction of the repetition period T, the 
most common case in practice. Thus, the assump- 
tion of a rectangular pulse, while not strictly 
valid, is quite good enough for most uses. 

The spectral width of the noise significantly 
affects the number of peaks exceeding our trigger 
level Vo or Ro. In fact, for noise alone and for the 
types of signal mentioned above, m(.)4.7 is directly 
proportional to the band-width Af of the noise mean 
power spectrum, whether broad- or narrow-band. 
From the physical point of view we would expect 
this, since the number of crossings at any level 
depends not only on the probability that the 
wave will reach or exceed such a value but also 
on how often the wave fluctuates, which, on the 
average, depends on the number of times the 
disturbance passes through zero. This number, 
in turn, is of the same order of magnitude as the 
mean frequency of the noise and hence is pro- 
portional to its band width. The constant of 
proportionality depends on spectral shape, al- 
though the variation from one distribution to 
another is not very pronounced. Table | lists 
data on several types of spectra considered in 
this paper. 

For purposes of comparison we normalize our 
spectra so that they all have the same maximum 
intensity W» and equal power; wr is accordingly 
different in the various cases. From the table 
we see immediately that 


(epg haus = ——(wr) rect. = ()*(wr)opt., 
aw) 


where 
(wr) (a) =2(wr)- (1-2.6) 


Then the band-width factors p are related by 


: 2\! 6\! 
(p) gauss = (-) (wr) rect. = (-) Prect.; 
T T 


Prect. = WF, (3)3. (1-2.7) 


It is interesting to note that for the ‘‘optical’”’ 
spectrum (single-tuned circuits), (3) Table I, p 
becomes infinite. The spectrum does not tail off 
to zero intensity with sufficient rapidity. In 
actual practice, however, this catastrophe does 
not occur, as damping or parasitic effects insure 
convergence; the ideal ‘“‘optical’’ spectrum is 
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modified so that the contributions to the band- 
width factor p for frequencies well away from 
resonance can be completely ignored (see refer- 
ence 1, Section 3.3 for further details). 

Curves showing the variation of the expected 
number of peaks per second with trigger level 
are included in Figs. 2 and 3. In general, this 
number is greatest at triggering levels for which 
the probability density of the amplitude or 
envelope is a maximum, or, approximately, for 
levels which are equal to the mean value of the 
disturbance. In the case of broad-band noise 
these points are V)>=0, or Vo=V, if there is a 
steady component (Fig. 2). For the envelope R 
of narrow-band waves the mean value (R) is 
always greater than zero, as Fig. 3 indicates. 
Note, moreover, that for triggering levels Ro 
below R there are actually fewer peaks/sec. than 
at Ro=R; maxima and minima above this level, 
and in the time interval between an upward 
crossing of Ry and the next downward one, 
merely act as a single peak exceeding Rp as far 
as triggering effects are concerned. Gaussian 
noise spectra, with the same power and maximum 
intensity as rectangular ones, always give a 
greater number of noise peaks. 


Improvement 


We are now able to determine the improve- 
ment (db) in the minimum allowable trigger 
level of one circuit which can tolerate A, spurious 
pulses per second over a simpler one which can 
tolerate only A» spurious pulses (A2<A;). Our 
circuits are assumed to have otherwise identical 
characteristics: (b9);=(bo)2 and pi= po. 

There are various ways in which ‘improve- 
ment’’ may be defined, depending on the mode 
of operation of our circuits. For example, we 
may use the following definitions: 

Let 


I,= {output level (after detection, if any) for 
which the expected number of spurious 
pulses/sec. for circuit No. 1 is A;}. 


Iz={output level, as above, for which the 


expected number of spurious pulses/sec. is 
A: (A2<A))}. 


Then the improvement is 


TY =20 logio(I2/I;) db; (I2>J]1; A2<A31). (I-2.8) 
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If the d.c. is excluded, as is frequently the case in 
audio or video filtering following a detector, our 
definition (1-2.8) is modified to 
T2—(La.c.)2 
r'(—d.c.) =20 logo —— lao > 1). (1-2.9) 
"- Tae.) 

The relation between the expected number of 
spurious pulses/sec. (A:, A2) and the triggering 
levels (11, I2), (Roi, Roz), etc., in our definitions 
of improvement may be obtained from (I-1.1), 
(I-1.2), and (II-3.16). For pulses (of any shape) 
whose duration 7 is short compared to the 
repetition period T (T>r), so that the effect of 
the signal is negligible, we have from (I-2.8) and 
(11-3.19) on solving for Vo; and Voo, 


Vor 
r=20 logu( ) 
Vor 


=10 logu( 





logio27 ye 








)ao (1-2.10) 


logi27y1 
for unrectified, broad-band waves. Here 
y=[A+a(vo,a0;T)]/p (<1/2m), (I-2.11) 


where a@ is a correction term caused by the 
transient effects of the leading edge of the pulse 
(see II-3.18). Figure 4 shows I as a function of ye. 
































P wb) 


Fic. 4. Here y: and ye are factors depending on the 
expected number of spurious pulses per second, on the 
period of the modulation, on the input carrier strength, 
noise level, etc. The precise relationship is given in Eq. 
(I-2.11) and the “improvement” (db) of circuit No. 1 over 
No. 2 for broad-band noise and short pulses follows from 
Eq. (I-2.10). 
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Fic. 5. Improvement (db) of circuit No. 1 over No. 2 
after half-wave linear rectification of narrow-band noise 
and short pulses d.c. component of the envelope removed. 


A more important case is that of narrow-band 
signal and noise triggering a circuit after rectifi- 
cation by a half-wave linear or quadratic de- 
tector. For the moment we assume that the d.c 
is passed by the video or audio filter. The total 
low frequency and d.c. outputs of the respective 
rectifiers are’ 


(linear): In=B8R/m2; Ia-c=68(bo/27)', 
(quadratic) : In=8R?/4; Ia-.=Bbo/2. 


(1-2.12a) 
(1-2.12b) 


Here £ is a tube constant. Again the signal is 
assumed to be on a negligible portion of the 
period. Our result (I-2.4) may thus be used in 
(1-2.8) to give us 


To2 
Lain.) = 20 logu(—) 
Tou 


y2exp (7 02”) 


= 20 logw( ) = $T quaa.(db), 


yrexp(ro1”) 


(rox > Tor; yo), (1-2.13) 


where now the transient term a(vo,@0; 7) in y 
(cf. I-2.11) has the value 1/7. When the d.c. 
component is filtered out, which is the more 


7 Reference 1, Section 4.1. J;; may also be obtained 
from (II-2.1). 
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common case in practice, we must use the 
definition (I-2.9), which becomes with aid of 
(I-2.12) 





roo — (w)!/ -) 
ro —(m)*/2 


L ain.) ( —d.c.) a 20 logw( 











ye2exp (1o2”) — 1/2 
= 20 logu( ), (1-2.14a) 
yrexp(ro:") — 1/2 
and 
Yo2" —1 
Tgeuna.)( —d.c.) = 20 logu( ) 
Yoi- —1 
y2exp( 2702”) — 1/2 
= 20 logu( ). (1-2.14b) 
yrexp(2ro?) —1/m 


The various values of 7»9 are obtained most 
easily from a plot of y=(A+a)/p as a function 
of ro, viz: 
ro 
y=— 
ee”? 





exp(—7o"). 


From this relation (Fig. 7) we note that y;,2 
<(27)-!, so that T must exceed zero. 

Curves for the improvement following linear 
or quadratic detection are shown in Figs. 5 and 
6 as a function of ro: =702/8, where B (>1) relates 

















_Fic. 6. Improvement (db) of circuit No. 1 over No. 2 
after half-wave quadratic detection of narrow-band noise 
and short pulses, d.c. component of the envelope removed. 
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the critical envelope levels ro; and ro2 at which 
triggering takes place in the two circuits. We 
observe that, for the cases when the d.c. is passed, 
if we set our triggering level below the point of 
most probable amplitude (~ 79), i.e., 701 < (70) max 
<ro2, the improvement decreases from its maxi- 
mum value at 7o:=(7o)max until 79; reaches a 
value for which y;= ye, at which [ =1; 79; below 
this level gives a response for circuit No. 1 less 
favorable than for circuit No. 2. When the d.c. 
(7) is filtered out, however, this situation cannot 
occur, since (7o)max <7o, and 71> 7. 

In general, the improvement is not great (a 
few db) if the trigger levels of the two circuits 
are not close to the most probable value of the 
amplitude and if the ratio of the number of 
spurious signals allowed/sec. (Ai, A:2) to the 
bandwidth of the disturbance is small. This may 
be observed from Fig. 4, or from Figs. 5 and 6 
used in conjunction with Fig. 7. Significant 
improvement (>10 db) follows only if circuit 
No. 1 can tolerate a much larger number of 
spurious pulses/sec. than circuit No. 2 and if 
this number is the same order of magnitude as 
the bandwidth factor p. 

The author wishes to express his gratitude to 
Miss Marilyn Lang, who performed the compu- 
tations. 



































Fic. 7. The parameter y is given by Eq. (I-2.11) and 
has the same significance as in Fig. 4; ro is defined as in 
Fig. 3. 
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II. MATHEMATICAL SECTION 


1. The Expected Number of Noise Peaks 
Exceeding a Given Magnitude 


The expected number of noise peaks that 
exceed a given level may be determined in the 
following fashion. We begin by calculating the 
probability that a function 


y=y(ay---ay;by:--by;x) = (II-1.1) 


passes through the value £& in the interval 
(xo, Xo+dx) with positive (or, negative) slope. 
Here the a;, 0; (t=1, 2, ---w) are random param- 
eters, possessing certain statistical properties (cf. 
Section 3). We choose dx so small that the parts 
of all but an ignorable (infinitesimal of higher 
order) fraction of the possible random curves 
lying in the strip (xo, x»+dx) may be regarded 
as straight lines. If y=yo at x=xo and passes 
through y=é (at x=x) for x9<x<xo+dx, the 
slope 7 of our line becomes 


n=(E—yo)/(x—Xo) ; (11-1.2) 


see Fig. 8 above. Since x» <x <x9+dx we obtain 
the inequality 


0> (yo— §) > —ndx. (II-1.3) 


Now let W:(yo, 7; x0) be the probability density 
function for the random variables yo and 
n=0y/dx|z-29 for a specified value of x=Xo. 
Then the probability that yo and 7 satisfy our 
inequality (II-1.3), going through y=£é with 
positive slope in the interval (xo, xo+dx), is 


~ yo=t 
dx f If Wi(yo, 7; sad in 
0 


vomt—ndz 
= p(t; xo)dx, (II-1.4) 
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which gives us 


d pf 
as - — -f dn f Wi(yo, n; wd] 
d(Ax)J t—ndz A200 


D d | 
ada | del -——(—vdx 
| f 7 dias) - »)| 


KX Wi(E— nAx, n; «)| 


427-0 


after differentiation in the limits; and thus we 
find® ** 


pies xo) = f nWilé, n; xXo)dn. (I1-1.5) 
0 


In our present case we have y= V(t), whose 
nature will be specified more precisely in the 
following sections. The expected number of 
crossings of y=t=V, by V(t), with positive 
slope in the time interval t—f, or what is the 
same thing, the expected number of ‘‘peaks”’ in 
excess of Vo in that interval, is given explicitly 
by (I-1.1) on a simple change of notation. The 
expected number of peaks per second is obtained 
from (1-1.2). 

Our problem has now been reduced to calcu- 
lating the probability density p(Vo;t) for a 
random noise wave, or one that contains noise 
and a signal, and evaluating the integrals (I-1.1) 
or (1-1.2). 


2. The Effect of Rectification 


We examine first the case of broad-band noise; 
see Fig. la. Let J(t) be the output (signal and 
noise) current after rectification by a device the 
Fourier transform of whose dynamic character- 
istic g(V) is f(iz). The output may be expressed 
as 


1 
1) =— f flisde" de, (11-2.1) 
2r vc 


* This is a generalization of Rice’s result, Eq. (3.3-5), 
Bell Sys. Tech. J. 24, 46 (1945), in which he computes the 
expected number of zeros, i.e., =0. In this instance 
~ (€=0) results equivalent to those of Rice (see also Am. J. 
Math. 61, 409 (1939)) have been obtained independently 
by M. Kac, Bull. Am. Math. Soc. 49, 314 (1943) for 
somewhat less stringent conditions on the random param- 
eters (a;,b;). See also M. Kac, Am. J. Math. 65, 609 (1943). 

** Since this was written, the result (II-1.5) has been 
stated by Rice, Bell Sys. Tech. J. 27, 109 (1948), Eq. (2.1); 
see also the remarks in footnote 10 of that paper. 
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where C is the familiar contour®* extending 
along the real axis from —« to + and is 
indented downward about a possible singularity 
at the origin. It is further understood that any 
video or audio filter following the detector is 
essentially uniform and infinitely wide for our 
purposes. 

From (11-1.5) the probability density for the 
expected number of peaks in excess of an arbi- 
trary output J) (>0) is 


pit)= f IP(Io, 1; t)dI, (11-2.2) 


where J and J are related to V and Vo by 
(I1I-2.1), and P is a density function analogous 
to W,. Differentiating (II-2.1) we observe that 
since V= Vo, dV =0, we have finally 


a 
pilest)= f (—) VP(1o, I; t)dV. 
0 OV/ 1=19 


Our assumption of no distortion because of 
frequency selection requires that the total prob- 
ability for the number of peaks before and after 
detection be unity. This permits us to write 


(11-2.3) 


’ a(I, 1) | 
P(Io, 1, t) =Wil Vo, v9 /— ad 


ol 2 : 
-(—) Wil Vo, ¥:h. 
OV/ 1=19 


When (II-2.4) is substituted into (II-2.3) the 
result is 


(11-2.4) 


p(1o; t) -{ VWi(Vo, V; dV 
0 


=p(Vo;t) (II-2.5) 


as expected. © 

For narrow-band noise (Fig. (1b)) and a signal 
in its immediate spectral vicinity the same 
general procedure may be used to show that for 


the low frequency output J1.r., i.e., the rectified 
envelope, 


p(L Tur. jo; t) -{ RW,(Ro, R: t)dR 


= p(Ro; t), 


ry See, for example, reference 1, Sec. 4.3, and Appendix 
4A, 


(11-2.6) 
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where R (>0) is the input envelope, Ro the level 
exceeded, and W,(Ro, R;t) is the probability 
density obtained from W,(Vo, V;2) on trans- 
forming from (V, V)-space to (R, R)-space. 


3. Peak-Number Distribution Density for Signal 
and Broad-Band Noise. 


Let our disturbance, consisting of a signal 


component V,(t) and noise Vy(t) be represented 
by 


V(t)=Va()+ V6.) =Y=Xi+ Vo 


=V,+)> (a, coswrt+b, sinw,t), 


n=l 
wn=2rf,r, (II-3.1) 
and 
V(t} =Vyt+V.= ¥2=Vet+X2 
= v.45 (—Aawn SiNW al 
n=1 
+b,w, cosw,t), (I1-3.2) 


in which a, and b, are the familiar random 
parameters® having the properties 


(aib;))=0; (a:)=(b) =0; 


(a,a ;) = (b:b;) = 6 ;'Afw(fi). (I1-3.3) 


Here 6;' is Kronecker delta, such that 6,;'=0, 
i~j; 6,/=1, w(f;) is the spectral intensity (mean 
power density) at the frequency fi, and Af is 
the separation between adjacent components in 
the Fourier development (see II-3.1, 3.2). The 
brackets ( ) indicate a statistical average. The 
quantities X, and X_2 are the gaussian random 
variables whose distribution W,(X,, X2) we seek. 
In terms of Y;= Vo, the level to be exceeded, and 
Y,=X.+V,, we observe that 
1 


Wi( Vo, Y3; t) => ation —exp[ — ( Vo- V.)? 2bo 
22 (bobe)! 
—_ ( Y.- V.)2, 2be |; (I 1-3.4a) 


where 


b= f w"w(f)df. (I1-3.4b) 


0 


W, is now a function of ¢, in as much as V,= V,(¢) 
and V,=V,(t) in the general case. Equation 
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(1I-1.5) or (11-2.5) then integrates to 
Pp r 

P( Vo; t) on (Vo— Vz)?/2bo | 
T 


X {exp(—V.?/2bop*) + Va(m/2bop?)! 
X(1+O[V./(2bop*)#])}  (II-3.5a) 


or 





(2 opt (Vo— Va)?/2bo ] 


V. 1 
x + 
2(bop”)* (2)! 





» XiFi(—4; 45 — Ve2/2bop?) , (II-3.5b) 


where p=(b2/bo)'. The quantity © is the error 
integral defined by 


2 zx 
OW) = f exp(—y")éy, 


(4 

and tabulated, for example, in B. O. Pierce’s, 
A Short Table of Integrals (Ginn and Company, 
New York, 1929). The expression ;F; is a con- 
fluent hypergeometric function, whose more im- 
portant properties are listed in references 1 
and 5. 

When there is no signal, we obtain Kac’s* ® or 
Rice’s result (reference 8, Eq. (3.3-13)), namely 


p 
p( Vo) = =n — Vo? /2bo) =exp( — Vo?/2bp) 
7 


<[% the expected no. of 


zeros/sec. |. 


(I1-3.6) 


Since now the process is stationary, p(Vo) is 
independent of ¢. 

Notice in the general case (II-3.5), that the 
probability density becomes infinite (~V,) for 
signals possessing an infinite time derivative. 
The physical significance of this point is dis- 
cussed more fully as we consider (a) and (6) 
below. 


a. D.C. Signal 


Here the signal appears simply asad.c. voltage. 
Then V, vanishes, except for the leading edge, 
marking the point (in time) at which the signal 
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is first applied. Let us assume for the moment 
that the change from no signal to a steady 
component of amplitude A» takes place linearly 
so that we have 


V(t) =Aot/t0; 
=Ay 


From (II-3.5) the probability density for the 
expected number of peaks/sec. in excess of Vo in 
the transition region (0O<t<7») is 


O0<t<70; 
t> rT; 


Ve =A o/To 


=0 (II-3.7) 


p( Vo; t) = —exp[ ~ (v9 — do ‘t9)? ] 
TT 


X { (projexp[ — ac?/ (p70)? ] 
+ado(r)*L1+0(ao/(p70)) ]}, (11-3.8a) 


where vo= Vo/(2bo)* and ag=Ao/(2bo)!, and in 
the steady-state interval (t>70) we have 


b( Vo) = (p/2m)exp[ — (vo — ao)? ] = (b2/2)! 
(probability density for 
noise distributed about 


the average value A»). (11-3.8b) 


The probability density is greater than it would 
be without signal, since the noise is ‘‘boosted”’ 
to a new average level (vp—dp). 

Observe in (II-3.8a) that as to—0 (t-0 also, 
since 0<t<70), p( Vo; t) becomes infinite as ro7!. 
The quantity (p70), which is proportional to the 
ratio of filter width to signal width (freq.), 
however, remains finite and possibly different 
from zero as to>—0. Usually (p70) vanishes as 
to—0, since the amount of noise power is assumed 
to be constant, with non-zero intensity Wo; this 
means that p=a finite constant. It may be 
objected that to achieve the linear transient of 
our example (II-3.7) above, even when ro>0, 
requires an infinite filter width. This is true, but 
the slope of the transient remains finite, and so 
p(Vo;t) presents no singularities. In practice, 
moreover, our filters never pass the highest 
harmonics, with the result that the slope of the 
transient, or signal, is always bounded. The 
quantity (p70), or an equivalent expression for 
_waves of other shapes, never quite vanishes, 
although it may become very small in the case 
of filters with high resolution. In these instances 
the limiting form as t>—0 is a good approxima- 
tion, i.e., the ideal vertical edge is a safe model 
for the real, not-quite vertical one. 
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In an interval 7(> 70) from the time the signal 
first goes on, the expected number of peaks per 
second exceeding Vo is, from (I-1.2), 

(V.) {1 =| [—(m»—a0)*} 
No+N =—| 1—— Jexp[ — (%—a»)? 
+N 0 oe T 0 


+ [4 Tao(x)! |" O(a) — ©O(v — ao) 1 
X (proexp(—ao?/(p70)*) 


+ado(r)'[1 +0(ao/ pro) }). (I1-3.9) 


For long times (T>7o) the transient’s effect 
becomes quite negligible, as we would expect, 
and ,,n(Vo) is now stationary and reduces to 
(1I-3.8b). The precise condition when p( Vo) may 
be used in place of ,,.(Vo) for signals of the 
above type may be obtained from (II-3.9) on 
requiring that p/2mrexp[ —(vo—ao)*] exceed the 
remaining terms by a sufficient amount. 

In the case of a vertical leading edge our result 
(I1I-3.9) assumes a simpler form since as to—0, 
po vanishes also, in consequence of the condition 
of finite noise power. We have accordingly 


p 
No+n( Vo) =—exp[ — (v9 — ao)? | 
2r 


1 
+—[O(v») = ©(v — dp) |. (1I-3.10) 
2T 


Again, for sufficiently long times the contribution 
of the second term proves negligible and n,,~( Vo) 
= p(Vo), Eq. (II-3.8b). 


b. Signal As a Trapezoidal, Triangular, or 
Rectangular Pulse 


Before considering a rectangular pulse let us 
evaluate p(Vo;t) and n,,n" for a trapezoidal 
pulse of width +, with the following analytic 
representation : 


V,=Aot/To, O<t<To, V.=Ao/To 
=Ap, tT9<t<1t—To, =0 . 
=A (r—t/r0), r—T0St<7, = —A>/t 

(I1-3.11) 


The rectangular pulse, of duration 7, may then 
be treated as the limiting form of the above as 
to—0. From (II-3.5) we see that the expressions 
for p(Vo;¢t) are precisely those given in (II-3.8a) 
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and (II-8b); the former applies in the interval 
(O<t<7o), and during the time (r—19<t<r) 
provided we replace t by r—?; the latter holds 
for the period (r><t<7—7 9). The expected num- 
ber of peaks/sec. exceeding Vo while the pulse ts 
on follows at once from (I-1.2) and is 


p 270 
Ne-+ v( Vo) = —{ 1— —Nexpt — (Up — dp)? | 


7 ee 

Qrao(x 

X L(pr0)exp( —ao?/ (p70)*) 
+o(1)!O(ao/ pro) |. 





[ O(v) — O(v9— do) | 


(II-3.12) 


When the transients fore and aft are short com- 
pared to the total duration of the pulse (17>>27») 
we obtain the familiar result (II-3.8b). 

An expression similar to (I1-3.12) for triangular 
pulses may be deduced immediately on setting 
t=27 9; we have then 


1 
——[O(v) — @(v9 — av) ] 
4 r9o(2)! 


Noin( Vo) = 


x [ (pro)exp( —Ay (pto)”) 


+ao(r)'O(ao/pro) | (I1-3.13) 


while the pulse is on. 

The most important special case is the rec- 
tangular pulse. Here rt» and top both vanish in 
the limit of finite noise power and vertical pulse 
front, giving us finally from (II-3.12) 


p 
Nein ( Vo) = 4 —exp( — (v9 —ao)?) 
2r 


[O(v) — O(v — ao) | 


2r 


(11-3.14) 


The expected number of peaks exceeding Vo 
while the pulse is on is rm,,n~(Vo), cf. (11-3.14) 
above. 

For a negligible transient contribution we 
require the pulse width to be related to the 
noise band width by 


T 
T>>—exp (09 — do)? O(v) — O(vo—ao) J. (I1-3.15) 
p 
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Now, the expected number ,, of peaks/sec. 
exceeding V» during a single pulse is given by 
(II-3.14), and when there is no pulse, mor is 
expressed by (II-3.6). Let T be the repetition 
period of the pulses. The expected number ms 
of spurious peaks per second, greater than Vp, is 


ns= {total no. peaks/sec.—no. signal peaks/sec. 











in absence of noise}, (11-3.16) 
=(r/T)mont+(1—17/T)morr—1/T ) 
~(: T 
oe f .(11-3.17) 
x Lexp(—ax-+20w) -1]) 
| — a(vpo, ao; T)3 
where 
a(Vo, ao; T) = (1 T)A —_ © (v0), 2 
+O(vp—ado)/2). (II-3.18) 


Now if the pulse is on only a short fraction of 
the period (J7>r), which is the customary case 
in most instances, (II-3.17) simplifies to 


pexp( — v0") 
Rex - ———a(Up, do; t), 


(11-3.19) 
2n 


as we would expect, since the number of added 
noise peaks while the pulse is on is negligible 
compared with the number exceeding Vo in the 
interval T—+r. Note that in the limit of large 
signals (ao?>>v,") the transient effect gives a the 
value 1/2T, while for small signals (vo?>a¢?)a 
becomes 1/7, as then the contribution of the 
transients can be ignored. 

Other pulse shapes, such as the gaussian or 
“optical,”” may in principle be handled in the 
same way, with the aid of (II-3.5), but the 
integration does not prove readily tractable, 
which is the pragmatic reason for our choice of 
the rather idealized pulse shapes treated in the 
present paper. 


4. Peak-Number Distribution Density for 
Narrow-Band Noise 
In this case (see Fig. 1b) our voltage wave 
V(t) may be represented by (II-3.1) with 
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Wn=Wn +wo, where wo is an arbitrary (angular) 
frequency in the band, usually the central one 
if the power spectrum is symmetrical, and w,’ is 
a “local” frequency, referred to wo. Then we have 


(11-4.1) 


V(t) = V c(t) COSw ol + V s(t) SINWof. 


Here Ve and Vs are, respectively, the slowly 
varying parts of the noise wave in phase with 


W(X 1X2X3X4) = 


—— exp} — 
42? (bob. — b,”) 
where, in analogy with (I1-3.4b), 


zx 


b= f (w—wo)"w(f)df. 


(11-4.3) 


Now, since the noise is narrow-band, we may 
introduce the concept of the envelope and write 
(11-4.1) as 


V(t)=R(t) coslwot—O(t)]. (1 1-4.4) 


The quantity R is the envelope and @ is a phase, 
both of which are slowly varying functions of 
the time compared to coswof. Comparison of the 
coefficients of coswof and sinwof in (11-4.4) and 
(11-4.1) shows at once that 

R cos@= Ve=Xi1; 


Rsind@=Vs=Xe. (I1-4.5) 


The Jacobian of the transformation from X-space 
to (R, R; 6, 6)-space is R*, so that 


W(X 1X 2X3X,) = R2wi(R, R: 8, 6) 


=W,(R, R: 0,6), (11-4.6) 


where w,(R, R; 6, 6) is the density obtained on 
substituting for X,, Xo, etc., their values in R, 86, 
etc., obtained from (I1-4.5) and its time deriva- 
tive. To obtain W,(R, R), the probability density 
of the envelope and its rate of change, we must 
integrate over 6 and 6, the ranges of which are 
respectively (0O—27) and (—* to +). From 
(11-4.2)—(11-4.6) we obtain 


R . 
—_—exp(— u*R?/2| M|) 
(27y*)! 
-exp( — R*(u!! uw — (Cul?) /2u| M |], 


where 


Wi(R, R) = 


(11-4.7) 


| M | = (bob. — b,”); 
us =bo| M|!; 


p'=be| M\}; 
pit =, | M\ i 
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1 — be Ve2+ Vs?) —bo( Vie? + V 32) —2b1( VeV s— VsVe) 





coswof and sinwof. To determine the probability 
density W, we require now the joint distribution 
of the four random variables X,;= Vc, X2= Vs, 
X;=Vc, Xs=Vs. The approach is exactly the 
same in the case previously considered, of broad- 
band noise (Section 3), but extended to four 
variables instead of two. The desired probability 
density is® 


. : , (IT-4.2) 
2 (bobs — 6,7)? 


When (I1-4.7) is substituted into (II-1.5) we 
obtain for the probability that the expected 
number of envelope peaks exceeding Ro 


| M| Ro 


Pi Ro) = —exp[ — Re ( ga™* 


(2e(u*})! 
—(u'}) 2u*| M pA 


This simplifies somewhat when we observe that 
b,=0 because of the narrow-band condition and 
assumed symmetry of the noise spectrum. Our 
result (11-4.8) reduces to 


(11-4.8) 


be 4 Ro 
p( Ro) = (~) ———— —exp( ae Re’ ‘2bo) 
bo (2mbo)! 


- ( “) X (Rayleigh probability 


2r density) 


(11-4.9) 


as the expected number of crossings per second, 
or equivalently, the expected number of envelope 
peaks per second greater than Ro. Again because 
of the stationary character of the process, p(Ro) 
is independent of the time. 


5. Peak-Number Distribution Density for 
Narrow-Band Signal and Noise 


The results of the preceding section may be 
extended without difficulty to include the case 


of a modulated signal. The 
question is now represented by 


V(t) = Vy(t)+ V(t), 


disturbance in 


where 


(11-5.1) 
V(t) =Ao(t) cosw-t 


and where Vy is specified by (11]-4.1). Here we 
let w.-=wo+twza, so that wy is the difference be- 
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tween the central noise (angular) frequency wo 
and the carrier frequency w,. Then it is evident 
that 


V(t) =(Veta) coswot+(Vs—B) sinwot 


=Rcos(wot— 6), (II-5.2) 
with 
a=Ao(t) coswat; B=Ao(t) sinwet; 
R=((Vce+a)?+(Vs—)*}; 
' : (I1-5.3) 





Vs—B 
é=tan™ , 
Veta 
the only additional requirements being that A o(¢) 
and A(t) do not vary faster than V¢ or Vs and 
that wa<wo, w-, as otherwise we would no longer 


be able to consider R as an envelope. As before 
(Section 3) we may set 


Y,=Veta=Xi+a; Y2=Vs—B=X2—8; etc., 
and use (I1-4.2)—(I1-4.3). The Jacobian of the 
transformation from X-space to (R, R, 0, 6)-space 
is again R*, where now R is given in (II-5.3). 
We obtain finally the following general expression 
for W,(R, R; 6, 6;2): 


DY 


W,(R, R; 0, 6; 1) =————exp(— [Afi (R, 9, 2) 
4n’| M|! 


+Bf2(R, R, 6, 6, t) 


+Cf;(R, R, 6,6, t)]), (11-5.4) 
where 
A=p"/2|M|; B=yp®/2|M|; 
C=2y'*/2| MI, (11-5.5) 


and 


fi=R’+0°+6?—2R(a?+6*)! cos(0+¢), | 
fo= R°+OR?+c++(2R6 sin(6+y) 
—2R cos(6+y) ](a2+6)!, 
fs= R°6—[R6 cos(0+¢) 
+R sin(6+¢) ](a?+6*)!+48 — Ba 


eee 























+R(c?+ 6)! sin(@+y), 
with the angles ¢ and y defined by 
a B 
cos¢ = , sing=— . 
(a? +")! (a? +)! 
(11-5.7) 
a 
cosy = ————_,_ siny=- —, 
(a?+ 6")! (a?+ 6")! 


The desired probability Wi(R, R) requires inte- 
gration over @ and 6; the latter is easily per- 
formed, but the former proves excessively com- 
plicated in the general case. Hence we introduce 
the approximation };=0, which makes C=0; 
see (II-5.5). Integrating over 6 and expanding 
the trigonometric exponentials in a triple series 
of Bessel functions, we may effect the integration 
over 6 with the aid of 


1 2r 
im=— [ cosj(@+ 7) 
2r 0 


Xcosl(6+ 5) cosm(0+7)d0 
1 
=— {6:-m/ cos[/(5— 1) —m(y—1n) J 
2€; 
+ 514m? cos[1(5—) +m(y—1) ]. 


If we note that \=¢—W= —weA o(t)/Aol(t), we 
find at last that the probability density W, is 


(11-5.8) 


R ;' 
Wi(R, } +. t) = eee a, = (R?+Ao(t)?), 2ho— ( R? + Ao(t)?/2 + w4?Ao(t)?/2)/2be | 


bo 2rbp 


wa oa 





em (E+E) (Rae +6")! 
ee =1( ):.(——— 


k=0 m=0 2 2b. be 


When the carrier is tuned to the center of the 
IF filter, i.e., wo =w,, then w,=0 and Eq. (I1-5.9) 
is altered accordingly, but the simplification is 
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RA, 
| cuendassn( ; ) cos(2k-+m)n 


0 


RAg 


+m milianni( 
b 





) cos(2k—m) ; (11-5.9) 


0 


slight. However, for an unmodulated carrier on 
center (w.=woe), considerable reduction is pos- 
sible, for then only the terms for m=k=0 
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contribute, 


(I11-5.9) 


and we may write at once from 


| R 
Wi(R, R) = ——— Io( RAo/bo)expl — (R?+Ae’) 


bo( WO2 


—R?/2b.], (11-5.10) 


which is independent of the time, since now the 
process is stationary (for the envelope and phase, 
but not for the entire wave). For no signal at 
all (Ao=0), Eq. (II-5.9) reduces further to 





Ro 
P(Ro; t) —— —[Ro?+Ao(t)* ] aimed J/4b2) 
émE(m +1, 2) sa? + 8? 
Sete ee 
k=) m=( 2m! 2b. 


RoAo(t) 
x €2k+4 m1 2% | _ cos(2k-+-m)d+ €2x— ml ex—m 


z Zz 


Xd 





0 


Again, when there is no modulation and wa=0, 
(I1I-5.11) assumes a particularly simple form, viz: 
Ro 


Ro) = p——— 


RoAo 
1o( ~ = )expl — (Ro? + Ao’) 2bo | 
(2xbo)! bo 


b. \! 
= (~) X (Rayleigh dist. for 


2r envelope V,+Vy). (II-5.12) 


R 
Wi R, R) carpe Ge R?/ 2bo— 


R?/2b-), 
bo(2b.)! /2bs] 


(I1-5.10a) 


a simplified version of (II-4.7). 

The probability density p(Ro;¢t) for the ex- 
pected number of times per second that our 
wave will exceed Ro follows at once from (I1-2.6) 
with the aid of Eq. (2). 13.3 of Watson, Theory 
of Bessel Functions (Cambridge University Press, 
Teddington, England, 1945). We have 


. 


a? + B? a+ B? 
): A(= m+1; —~) 
4b, 2b. 


RoAo(t) 
4a (11-5.11) 


0 





We may employ (II-5.12) when the carrier is 
modulated by a rectangular pulse,® for the time 
the pulse is on, and use (II-5.13) below, for the 
interval of noise alone: 


Ro 
P( Ro) = "Omb)i Pe a Re, 2bp |. (11-5 1 3) 


TOO 





Calendar 


September 


20-24 Illumination Engineering Society (National Technical Con- 
ference), Boston, Massachusetts 

International Rheological Congress, Scheveningen (The Hague), 
Netherlands 

27-29 Technical Association of the Pulp and Paper Industry, Poland 
Spring, Maine 

American Society for Professional Geographers, Madison, Wis- 
consin 


29-October 2 Institute of Radio Engineers, Los Angeles, California 


21-24 


27-30 


October 
5*+7 American Institute of Electrical Engineers, Washington, D. C. 
6-9 Society of Automotive Engineers (National Aeronautic Meeting 
and Aircraft Engineering Display), Los Angeles, California 
7-9 International Scientific Radio Union (American Section) and 
the Institute of Radio Engineers, Washington, D. C. 


12-16 Fifth National Chemical Exposition, Chicago, Illinois 

13-15 American Society of Civil Engineers (Fall Meeting), Boston, 
Massachusetts 

13-16 Electrochemical Society, New York, New York 

18-22 American Institute of Electrical Engineers (Midwest General 
Meeting), Milwaukee, Wisconsin 

20-22 American Standards Association (Annual Meeting), New York, 
New York 

21-23 Optical Society of America, Detroit, Michigan 
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of Meetings 


24-27 American Institute of Mining and Metallurgical Engineers 
El Paso, Texas 

Technical Association of the Pulp and Paper Industry, Buffalo, 
New York 

25-29 American Society for Metals, Philadelphia, Pennsylvania 

25-29 Society of Motion Picture Engineers, Washington, D. C. 

26 Association of Consulting Chemists and Chemical Engineers, 
New York, New York 

Gas Discharge Conference, Brookhaven National Laboratory, 
Upton, Long Island 

30 American Mathematica! Society, New York, New York 


25-28 


27-29 


November 

3-5 American Institute of Electrical Engineers (Southern District 
Meeting), Birmingham, Alabama 

4-5 Society of Automotive Engineers (National Fuels and Lubri- 
cants Meeting), Tulsa, Oklahoma 


4-6 Acoustical Society of America, Cleveland, Ohio 

4-6 Nationa! Electronics Conference, Inc. (Annual Technical 
Forum), Chicago, Illinois 

7-10 American Institute of Chemical Engineers, New York, New 
York 

8-10 Institute of Radio Engineers, Rochester, New York 

11-13 Geological Society of America (Annual Meeting), New York, 
York 

12 X-Ray and Electron Diffraction Conference (Southeast Section) 
Detroit, Michigan 

13-15 Conference on Applied Spectroscopy, Pittsburgh, Pennsylvania 
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Frequency Stabilization of Microwave Oscillators by Spectrum Lines. II* 


José L. Garcia DE QuEvEDO** AND WILLIAM V. SMITH 
Department of Physics, Duke University, Durham, North Carolina 
(Received February 2, 1948) 


The design of both narrow band high stabilization and broad band lower stabilization 
spectrum line discriminators is developed in a generalized form applicable to various spectrum 
lines and cavity designs. Comparison with experiment is made for two broad band discriminators 
utilizing the NH; 3,3-absorption line at 23,870 megacycles as the frequency determining 
element. These discriminators used with an amplifier of 2000 gain, reduce drifts by factors of 250 
and 1000, respectively, compared to an unstabilized tube. 





INTRODUCTION 


HIS paper is an extension of the work of 

the authors and two co-workers! on the use 
of the NH; 3,3-line at 23,870 megacycles to 
control the frequency of microwave oscillators. 
The theoretical analysis has been generalized to 
facilitate comparison with other spectrum lines 
and to include a greater variety of possible 
discriminator designs. A stabilizer design has 
been achieved that reduces frequency drifts by a 
factor of 1000 in comparison with an unstabilized 
tube—i.e., from long period drifts of a few 
megacycles to long period drifts of a few kilo- 
cycles—and means of obtaining still higher 
stability are discussed. The design of a broad 
band frequency standard based on microwave 
spectrum lines is developed. Saturation of the 
gas by the microwave radiation is found to be a 
major factor influencing the discriminator design. 


CIRCUIT THEORY 


The circuit theory developed in | for the 
admittance of a magic tee whose two symmetric 
arms are balanced by two equal lengths of wave 
guide, one length being filled with a low pressure 
gas with a resonance absorption at vo (I, Fig. 1b), 
may be reduced to the problem of a magic tee 
with symmetric arms balanced by two identical 


* The research described in this report was supported by 
Contract No. W-28-099-ac125 with the Army Air Forces, 
Watson Laboratories, Air Materiel Command. 

** This work is based on a thesis presented in partial 
fulfillment of the requirements for the degree of Doctor of 
Philosophy in the Graduate School of Arts and Sciences of 
Duke University. 

1W. V. Smith, José L. Garcia de Quevedo, R. L. Carter, 
and W. S. Bennett, J. App. Phys. 18, 1112 (1947). This 
paper will be referred to hereafter as I. A different approach 
to the problem is presented by W. D. Hershberger and L. 
E. Norton, RCA evles IX, No. 1 (1948). 
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cavities, one filled with gas, and both gas and 
cavities resonant at vp. A particularly simple 
cavity is a shorted wave guide of length /. The 
impedance of this cavity when filled with gas is 
found to be (J): 


Z=Z, tanhLao+j8+(y/2) ll, 


where ao represents the skin loss of the guide, 
B=2n/d,, A, being the wave-length in the guide, 
and y is the complex gas attenuation in cm 
(or twice the attenuation in nepers/cm) including 
both absorption and anomalous dispersion terms. 
From the theory of collision broadening,’ 


¥ = YoAwo/LAwo+j(w—wo) J. (1) 


At resonance y=‘yo, the maximum value of 
absorption, while at w—wo=Awo, the absorption 
is reduced by a factor of two and the dispersion 
is a Maximum. 

Expanding (1) around w—wo, and tanh2zl/dgo 
=0, for small values of a and y, where Ago is 
the value of \, at vo: 


mL lf dgoco 2 (w— wo) AwoyorA,go 
Z= +j + 


Ago L T Wo 








TWO 





1 
teccaciesals 
(2Awo ‘wwo) + j2(w = wo) /wo 
which is of the form 
1 2(w—wo) 
Z= 2] —+j— 
Q. 


wo 








1 
+ : | (3) 
ysl (1/Q,) + j(2(w—wo) /wo) ] 


2 J. H. Van Vleck and V. F. Weisskopf, Rev. Mod. Phys. 


17, 227 (1945). 
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where Z, = (#lZ)/Xgo; Oc = 4/ (Agomo); Vs = (24Q,)/ 
(YoAgo) ; Os =wo/(2Awo). Thus the shorted trans- 
mission line may be thought of as a series reso- 
nant cavity of parameters Z,. and Q, determined 
by the properties of the cavity without gas, in 
series with a parallel resonant cavity of param- 
eters Y.y, and Q, where y, and Q, are functions 
of the gas only. The fact that the gas is a parallel 
circuit while the cavity is a series circuit is a 
consequence of the anomalous dispersion prop- 
erties of the gas near resonance and means that 
the rate of change of impedance with frequency 
caused by the gas is in the opposite sense to that 
caused by the cavity alone. Thus filling a 
cavity of resonant frequency vo with a gas of the 
same resonant frequency will reduce rather than 
agument the effects of the cavity. 

If, however, the impedance of an empty cavity 
Z' is subtracted from (3), by means of a magic 
tee as shown in I, the input admittance to the 
tee becomes 


Y=1/(Z—2Z’') = Y-y.[.(1/Q.) +7(2(w— wo) /wo) | 

=[(1/a)+j(a/a)]Vo. (4) 
In the special case where the cavities are shorted 
transmission lines, with arms of lengths / and 
1+ ,/8, a discriminator curve may be obtained 
as in I, giving a slope of 


(dV/dv)|»=»0=DPLQ,/vo][20/(1+a)*] 1(7) 


volts per megacycle, where P» is the input power 
to the discriminator and D is the rectification 
efficiency of the crystal. More generally, the 
magic tee plus two cavities may be treated as a 





Fic. 1. Discriminator characteristic of 2.16 meter cell 
using the NH; 3,3-line. Vertical line establishes axis. The 
spacing between the center of the linear portion of the 
discriminator curve and the first satellite is 1.72 mega- 
cycles, establishing the frequency scale. 
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single cavity and then used as in I, circuit 1a, in 
which case Eq. I(7) is still valid. 

By adjusting the coupling to the cavity, 
maximum slope is obtained for a=1, in which 
case the gas is equivalent to a cavity of unloaded 
Q=Q.,. If cavity dimensions are greater than the 
mean free path in the gas and if the residual 
doppler line width is not approached too closely, 
Q.=wo/(2Awo) for the NH; 3,3-line may be made 
equal to 23,870/0.2 or about 100,000 with the 
attenuation still given by (1) and the discrimi- 
nator slope by 1(7). There is one important 
restriction, however. The r-f fields in the cavity 
must be low enough to avoid saturation, which 
implies either a large cavity or a low power 
input, as will be discussed in the following 
section. 


SATURATION EFFECT 


When high powers are passed through a gas at 
low pressure, it is found that the thermal equi- 
librium of energy states is altered, with a re- 
sultant change in the absorption at the resonant 
frequency and in the line width, so that® 


Yo 





[1+(AP/(Aw)2)] | 
2Av =2Avol1+(AP/(Avo)?) }}, 


(y)-0o= 
(5) 


where P is the incident power in ergs/cm? 
monochromatic at the frequency in question, 
. 


A =8x| p;;|2/3ch?. 


uig=matrix element of the dipole moment; 
vo= ((rh?v?)/(RTAvo))noA ; mo =number of mole- 
cules/cc in ground state. For the NH; 3,3-line, 
A is 1X107.*** While relations (5) were derived 
for plane traveling waves, they are very nearly 
true for traveling waves in wave guide if P is 
replaced by the v2 times the average power and 
both yo and A are multiplied by A,/A. For 
standing waves a second averaging over the 
guide length replaces P by twice the average 
power, or four times the incident power. 


3R. L. Carter and W. V. Smith, Phys. Rev. 73, 1053 
(1948); also C. H. Townes, Phys. Rev. 70, 665 (1946); B. 
Bleaney and R. P. Penrose, Proc. Phys. Soc. 60, 83 (1948); 
R. Karplus and J. Schwinger, Phys. Rev. 73, 1020 (1948); 
R. Karplus, Phys. Rev. 74, 223 (1948). 

*** Appropriate averaging (see reference 3) over the 
Zeeman components of the line replaces A by 1.710’. 
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If the coupling to the balanced cavities is 
adjusted to give a=1, all power incident at the 
resonant frequency is adsorbed. Assuming the 
gas absorption to be large compared to the skin 
losses, as is the case for the NH; 3,3-line and 
most silver cavities operating in high modes, 
then the Q, of the gas may be related to the 
cavity volume. V, the cross section ab of guide 
feeding the cavity, and the ratio k of energy 
density in cavity and guide as follows: 

2x Xstored energy 
= a 
energy dissipated per cycle 
€o0( cay”) Vw T Ago Vk 


2ab Peuide ‘ 2 a 








(6) 


in m.k.s. units. (In the special case where the 
cavities consist of sections of wave guide of the 
same cross section as the tee, k=} from the 
power division at the tee.) 

Substituting (5) in 1(7), for a=1: 








| DP, 
(dV/dy) | » =——— . (7) 
4APok X07! 
sanf 1+ | 
ab(Avo)? oN 


For fixed input powers, relation (7) increases as 
Avo decreases irrespective of the value of k, 
although it increases more rapidly as k ap- 
proaches zero, corresponding to a large cavity, 
low field in the cavity, and absence of saturation. 
If a small cavity is employed and the energy 
density within the cavity is large enough to 
cause saturation, the gas attenuation soon be- 
comes smaller than the copper losses and the 
above analysis no longer holds. 

While for a given volume the steepest dis- 
criminator curves are obtained when the coupling 
is adjusted to get a=1, for many purposes the 
broad band discriminator design employing 
balanced wave guide arms is desirable. By using 
tapered sections from the tee to arms employing 
larger sized guide, k may be made less than } 
while the broad band feature of the discriminator 
is preserved. For a given volume and input 
power however, a is no longer a_ variable. 
Assuming @ to be much smaller than 1, substi- 
tution in I(7) then gives 


(dV/dv)| »=DPoa/Av=KyoPo/Av, (8) 
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LINE WIDTH IN MC 


Fic. 2. Maximum discriminator slope and corresponding 
input power for 2.16-meter cell as function of line width. 
Curve A: observed maximum discriminator slope. Curve 
B: input power for maximum slope. 


where K=(D)goZ,)/(24Zo). Differentiating (8) 
with respect to Avo and equating the resulting 
expression to zero for a maximum gives 


4P okA Ago, ‘ab(Avo)*d = 3. (9) 


Converting units from c.g.s. to P in milliwatts 
and line width 2Av») in megacycles, with ab 
= 1.07 X0.43, and \,/A=1.2: 


Po/(2Avo)?=0.115/kA X 10-7, (10) 


Since the maximum value of discriminator 
slope occurs for AP/(Avo)*=}$, where saturation 
is just beginning, the ratio of signal to crystal 
noise is also near its optimum value. 
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Fic. 3. Schematic diagram of the laboratory model of the 
spectrum line stabilizer. The switch s is shown in the 
“stabilized” position. The cathode-ray oscilloscope was 
used to observe the discriminator curve, and the spectrum 
analyzer was used to measure the stabilization factor. 
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TABLE I. 











Frequency* Strength> 
Designation (megacycles) relative to 3,3 
2,2 23,722.63 0.4 
3,3 23,870.13 1.0 
4,4 24,139.41 0.6 
5,5 24,532.98 0.5 
6,6 25,056.02 0.9 





* W. E. Good and D. K. Cole, Phys. Rev. 71, 383 (1947). 
>» James W. Simmons and Walter Gordy, Phys. Rev. 73, 713 (1948); 
also C. H. Townes, Phys. Rev. 70, 433 (1947). 





In the experiments to be discussed, the cavities 
were shorted wave guides of varying dimensions 
connected to the arms of the tee by tapers, in 
which case substitution of (10) in (8) yields 


(dV/dv) | »=5 X10-*(D Vy0/ab) 


x[Pok/(AX10-7)]! (11) 


volts per megacycle, where the crystal rectifica- 
tion efficiency D is in volts per milliwatt and V 
is the volume of the cavity. The independent 
variable is chosen as Py as practical considera- 
tions, such as crystal saturation or tube pulling, 
limit its maximum value to about 10-' mw in the 
20,000 to 30,000 megacycle region. From the 
definitions of yo and A, the maximum discrimi- 
nator slope is seen to be proportional to the 
fraction of molecules in the ground state and to 
the effective dipole moment of the transition. 


COMPARISON OF EXPERIMENT WITH THEORY 


A discriminator was constructed using balance 
arms of 2.16 meters and 2.16+),/8 meters of 
1.07 0.43 cm silver wave guide. Photographs 
of the discriminator output as a function of 
frequency were taken by sweeping the oscillator 
through the NH; 3,3-resonant frequency. Several 
series of photographs were taken for different 
fixed input powers to the tee, and a continuously 
variable pressure (proportional to line width 
2Avo) as the cell was. evacuated. A typical 
photograph is shown in Fig. 1 in which the fine 
structure can be seen (establishing the frequency 
scale’) and a reference line shows the wcrtical 
-axis. The maximum and minimum discriminator 
outputs occur at the half-power points, so that 
their separation is the line width 2Ay. Also, the 

‘B. P. Dailey, R. L. Kyhl, M. W. P. Strandberg, J. H. 


Van Vieck, and E. B. Wilson, Jr., Phys. Rev. 70, 984 
(1946). 
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maximum discriminator slope at the resonant 
frequency can be shown to be twice the ratio of 
the difference of maximum and minimum dis- 
criminator outputs to the line width. The vertical 
scale was established by a separate calibrating 
signal not shown on the photograph. 

As the cell was evacuated, the maximum 
discriminator output continually decreased, at 
first very slowly, corresponding to the region of 
approximately constant yo, and then more 
rapidly as saturation was approached. During 
this time the line width decreased, so that the 
discriminator slope first increased, then went 
through a maximum as saturation was ap- 
proached, finally decreasing again at the lowest 
pressures. A plot of maximum discriminator 
slope as a function of line width is shown in 
Fig. 2, curve A being the observed points. Also 
plotted is the square root of the power at which 
the maximum slope occurred. Both relations are 
linear, as predicted theoretically. The plot of the 
square root of the power versus the line width 
yields an experimental value of P»/(2Avo)? =0.25, 
which from Eq. (10) gives a value of A of 
0.97 X 107. The best theoretical and experimental 
value of A from other investigations? is 1.7 X 10’. 
The discrepancy can at least in part be ascribed 
to the more complicated averaging of electric 
fields necessary in the discriminator investigation 
(an averaging that was done only semiquantita- 
tively). The maximum observed slope was 34 
millivolts per megacycle. The decrease of a factor 
of two from the figure of 70 millivolts per mega- 
cycle quoted in I is because the present figure is 
for one crystal only, whereas the previous higher 
value was for the output of two balanced 
crystals. 

A second cell of dimensions 1.07 X 2.40 X 386 
cm was also tested, and a maximum slope of 
152 millivolts per megacycle for a single crystal 
was observed. Theoretically, this large cell should 
give a slope steeper than the small one by the 
ratio ((386)/(216))(2.40/0.47)!=4.22. Different 
crystals were used in the two cases however, 
and the observed improvement with the large 
cell was somewhat greater than anticipated. 


FREQUENCY STANDARD DESIGN 


The small cell discriminator discussed in the 
previous section was built into a frequency 
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stabilizer according to the design of reference I, 
except that the correction signal was adjusted to 
be zero at the resonant frequency of the gas by 
balancing the output of the magic tee against 
the output of a directional coupler instead of a 
second magic tee (see Fig. 3). The method 
yielded a much smaller background of reflections 
than did the two-magic-tee set-up previously 
used, but resulted in a loss of a factor of two in 
sensitivity. 

The stabilization was measured by modulating 
the reflector of the 2K50 both when stabilized 
and when unstabilized, and comparing the fre- 
quency shift per unit voltage in the two cases. 
With the measured d.c. amplifier gain of 2000 
and tube tuning characteristic of 283 millivolts 
per megacycle, the observed maximum discrimi- 
nator slope of 34 millivolts per megacycle should 
give a stabilization factor of 240. The observed 
stabilization factor was 125, which is in satis- 
factory agreement since only half of the output 
signal of the amplifier is available for stabiliza- 
tion. Typical fluctuations in power supply volt- 
age of 0.1 volt applied to the reflector correspond 
to frequency fluctuations of about 1.5 kilocycles. 
A stabilizer based on the large cell would have a 
stability approximately fourfold greater while, 
by using appropriately large cavities and ad- 
justing the coupling for maximum results, much 
greater stability can be achieved in less flexible 
systems. It should be emphasized that higher 
gain amplifiers can also be used and that the 
modifications in the r-f circuit necessary to use 
an i-f amplifier® do not affect the basic discrimi- 
nator design. 

In applying the theoretical analysis of the 
preceding section to the design of practical fre- 
quency standards, it is evident that the require- 
ments of maximum stabilization and flexibility 
of operation are to some extent mutually incom- 
patible. Thus, a practical broad band system 
would have a cell of 10.7 4.3 mm cross section 
to match the tee cross section, and a length of 
about 6 meters (which can, of course, be coiled 
compactly), as beyond this length the sharpness 
of the discriminator curve is reduced by de- 
parture of tanh(ao+~7o0/2)/ from linearity with 
(aot+vo/2)l, at least for a silver cell and the 


5R. V. Pound, Rev. Sci. Inst. 17, 490 (1946). 
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NH; 3,3-line. The discriminator response is not 
critical to the eighth wave-length difference in 
balance.arm length, but some retuning is neces- 
sary for different lines. Some of the frequencies 
and discriminator characteristics available for 
ammonia are cited in Table I. To get the 
discriminator slope in millivolts per megacycle 
for the above cell, multiply column three by 100. 
The table is for the stronger lines only, and at 
room temperature. If the cell is immersed in 
dry ice, the density of molecules is greater for a 
given line width by (300/190)!=1.25, thus in- 
creasing yo and reducing saturation effects, 
while the relative intensities of the low J number 
lines is increased over that of the high J number 
lines. The ultimate limiting factor of the doppler 
width is also reduced at low temperatures. 

To date, ammonia is the only gas with an 
absorption in the centimeter region strong enough 
to be practical in a broad band frequency 
standard. It should be noted, however, that many 
transitions in the millimeter region are even 
stronger than the ammonia lines. 

The use of a cell of large cross section con- 
nected by a taper to the magic tee is an attempt 
to gain sensitivity while preserving the broad 
band feature of the design. This can indeed be 
done, but it should be noted that for a given cell 
length, the gain in sensitivity is proportional to 
the square root of the cross section when coupling 
is by means of a taper, as for a constant input 
power the field strength in the guide varies 
inversely with the cross section, and the optimum 
line width therefore inversely with the square 
root of the cross section. If size is no objection, 
however, large cross-section cells can be used in 
broad banded systems. Some trouble with ex- 
traneous modes of propagation was observed in 
the large cell discussed previously, but these 
could always be tuned out by squeezing the cell 
walls to slightly detune the resonance. 

For the most precise control at one fixed 
frequency, it is necessary, of course, to maintain 
the balance arms at the same temperature, 
although it is not necessary that the temperature 
remain constant. A more compact design, capable 
of better temperature control, might be accom- 
plished by using high Q cylindrical cavities of 
minimum surface to volume ratio, and coupling 
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to them rather heavily so that the reflection 
coefficient at resonance is 0.9 or more and a 
slight unbalance of the two cavities would be 
unimportant compared to the high Q gas ab- 
sorption. With this reflection coefficient and 
short sections between tee and cavity, a uniform 
but not necessarily constant temperature would 
again be necessary for high stability. 





In conclusion, the authors wish to thank Mr. 
Casimir Nawrocki for his invaluable assistance 
in pointing out the peculiar adaptability of the 
2K50 oscillator to this problem, and in helping 
familiarize us with the suitable operation of the 
tube. We also wish to thank Dr. W. O. Gordy 
and Dr. W. M. Nielsen for their continuing 
interest in the problem. 





On the Calculation of Radiation Patterns of Dielectric Rods* 


R. B. Watson anp C. W. Horton 
Defense Research Laboratory, The University of Texas, Austin, Texas 


(Received February 13, 1948) 


A previous paper described the radiation field of a dielectric rod used as a termination for a 
wave guide. The theory therein made use of a standing wave distribution of equivalent 
currents over the surface of the rod to describe its radiation effect. It is shown in this paper that 
the equivalent currents may be expressed in terms of traveling waves without appreciably 


changing the computed radiation patterns. 


HE theory for the radiation field of a 

dielectric rod as a termination for a wave 
guide has been developed and a summary of 
results presented in the paper, ‘“The radiation 
patterns of dielectric rods—experiment and 
theory.’ In this paper certain equivalent cur- 
rents distributed over the surface of the dielectric 
rod are used to describe its radiation effect. In 
the way in which the theory is there formulated, 
these distributions of currents take the form of 
standing waves over the surfaces of the rod. 
Since there is appreciable radiation occurring, 
it might be more appropriate to approximate 
these currents by expressions representative of 
traveling waves. The surface magnetic currents 
described in Eq. (9) of the above-mentioned 
paper should then take the following form: 


wy 
M=jM, —," (1) 


. The analysis proceeds as before except that the 


* The work described in this paper was done at the 
Defense Research Laboratory under the sponsorship of the 
ae of Ordnance, Navy Department, Contract NOrd- 

195. 

'R. B. Watson and C. W. Horton, J. App. Phys. 19, 661 

(1948). 
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integral J, [Eq. (19) ] has the form 


I 


I,= J eik’2e—ikz cosd] z (2) 
0 


which can be readily integrated to give 








singkl 1—cosgkl | 
q q 
or 
16 gkl 
| J;|2=— sin*—, (4) 
q° 2 
where 
g=n—cosé, n=k'/k. 


This is a much simpler expression than the prior 
value for J2, and the maxima and the minima of 
the absolute value are readily found by analysis. 
For the values of m=1.10, ki = 26.4, the function 
is found to have minima at @= 30.5°, 51.5°, 67.3°, 
etc. The function has maxima at @=0°, 40.6°, 
59°, etc. The corresponding values for the elec- 
tric field for the first three of these maxima are 
0 db, —12.5 db, and —18 db. Comparison with 
Fig. 4 of the reference shows that the differences 
of these values from the values there plotted are 
insignificant. The: terms other than J, in: the 
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calculation of the field remain the same, and 
therefore it is concluded that the traveling wave 
picture does not appreciably modify the radiation 
pattern as compared to the standing wave 
picture. The calculation for k/ = 64.2, because of 
the periodicity of the function sin?(gk//2), leads 
to a very small value for the maximum at 0° 
and a very large value, of the order of 20 db 
greater, for the next maximum which occurs at 
16.2°. This is in agreement with the calculation 
for the standing wave type of currents. 

A similar evaluation may be carried out for 
the equivalent electric currents on the sides of 
the dielectric rod. When this is done the integral 
I, takes the form of Eq. (2), and the evaluation 


proceeds in the same way. This process has 
already been shown to produce the same maxima 
and minima as the integral for the standing wave 
description, and hence the use of traveling waves 
for the description of the electric currents pro- 
duces no new results. 

It thus appears that neither a standing wave 
description nor a traveling wave description of 
the equivalent currents on the surface of the 
dielectric rod is satisfactory for describing the 
radiation patterns for long rods, that is, greater 
than about six wave-lengths, but that either 
description gives agreement with experimental 
patterns for rods in the range of length of three 
to six wave-lengths. 





Heat Transfer from a Rotating Disk to Ambient Air 


CaRL WAGNER 
Ordnance Research and Development Division, Sub-Office Rocket, Fort Bliss, Texas 


(Received February 19, 1948) 


The coefficient of heat transfer from a rotating disk to ambient air in case of a laminar 
boundary layer is theoretically calculated on the basis of von Karman’s calculation of flow 


conditions. 


INTRODUCTION 


EAR the surface of a rotating disk, air 

flows mainly in the tangential direction. of 
rotation. An additional component of airflow in 
radial direction is caused by the centrifugal 
forces due to rotation. Accordingly, ambient air 
moves from the surroundings towards the sur- 
face of the disk. Calculation of flow velocity as 
function of the spatial coordinates and the angu- 
lar velocity has been carried out by Theodore 
von K4rman! in order to obtain the torque mo- 
mentum due to friction. On the basis of these 
calculations, the coefficient of heat transfer from 
a rotating disk to ambient air in case of a 
laminar boundary layer can readily be obtained 
as shown below. 


MATHEMATICAL ANALYSIS 


According to von Karman, the velocity com- 
ponents, v, and 2, of the air flowing in the radial 
and the tangential direction inside the boundary 


1 Theodore von Karman, Zeits. f. angew. Math. u. Mech., 
1, 233 (1921); NACA Technical Mem. No. 1092 (1946). 
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layer are approximately* 


x x\? x 
s,=rw-| 1.026-—- (1--) (142°) 
6 6 6 
1 /x\? x\? 
-—(-)-(:--) | if x88 (1) 
2 \6 6 
1 x x\? 
=re-—-(24-)-(1--) if xsé. (2) 
2 6 6 


Here 7 represents the distance from the axis, 
x the distance from the disk, w the angular 
velocity, and 6 the thickness of the boundary 
layer which is independent of 7 and is given by 
the formula 


§=2.58-(v/w)! (3) 


where » is the kinerfatic viscosity of air. Equa- 
tions (1) and (2) hold true inside the boundary 
layer, whereas outside the boundary layer the 





* In Eq. (2) of the present paper an erratum in Karm4n's 
Eq. (27) is corrected; the term (1—x/5) must be em- 
ployed in the second power instead of the first power, as 
can be derived from von K4armén’s subsequent equations. 
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velocity components v, and v, vanish according 
to the above approximations. 

The velocity component v,, normal to the 
surface of the disk, is not explicitly given in von 
K4rman’'s paper, but can be obtained from a 
differential equation. Thus, inside the boundary 
layer 


Zz 


vy 
v-= -f 2-—-dx 
0 r 
z x sy" x 
= 20: f | 1.026--. (1-") (1422) 
s\ 5 é 
1 /x\? x\? 
---(-) -(:--) |-as if x6, (4) 
2 \6 6 


whereas outside the boundary layer the velocity 
component v, is constant, viz., 
C) 
v,= -f 2(v,/r) -dx = —0.708- (ww)! if x>6. (5) 
0 
At a steady state the general equation of heat 
balance for a fluid in motion is 


k-div grad’ —grad(pvc,3) =0, (6) 


where # is the difference between the local tem- 
perature and the bulk temperature of air, k the 
heat conductivity, c, the specific heat at con- 
stant pressure, p the density, and v the vector of 
flow velocity. The first term in Eq. (6) represents 
the increase of enthalpy due to heat conduction 
in the case of air at rest. The second term equals 
the gradient of the flux of enthalpy irrespective 
of heat conduction, given by the product of the 
vector pv of mass flow and the excess of enthalpy 
cp? per unit mass. In cylindrical coordinates, 
there is no dependency on the angular coordinate 
because of symmetry. Thus, 
fo lad ad 
div grad § =—_+-—-—+— (7) 


Or? ¢r or ax? 


Ata steady state the divergence of the product 
pv vanishes according to the condition of con- 

tinuity of mass flow. Consequently the second 
* term of Eq. (6) becomes 


grad(pvcpd?) =c,d0-divvp+vpc,: gradd 
ov ov 


=Vpc,* grad’ = pc,- (». —-+v,-— }. (8) 
or Ox 
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On introduction of expressions (7) and (8) 
into Eq. (6) and division by cpp, it follows that 


ed 
a-( + 
or 
where a represents the thermal diffusivity of 
air, given by 





1 ad dd dvd ov 
)-%- =v, =0 (9) 


r or ax? or Ox 


a=k/cyp. (10) 


The over-temperature #o of the disk shall be 
given so that the boundary conditions pertaining 
to differential Eq. (9) are 


v=, at x=0 (11) 


(12) 


According to von Karman, the thickness of 
the boundary layer is independent of the dis- 
tance r from the axis and the radial flow velocity 
at a given distance x from the disk is propor- 
tional to r. Consequently the total mass rate of 
air flowing in the radial direction at any distance 
r with the circumference 2xr is proportional to 
the area ar* of the circle of radius 7. Thus, a 
temperature distribution dependent on x but 
independent of r can be assumed, so the partial 
differential Eq. (9) becomes an ordinary differen- 
tial equation 


v=0 at x=. 


a(d?3/dx*) —v,(dd/dx) =0. (13) 


Upon first integration, 
(d3/dx) = ~exo| (1 a) f eds+A | (14) 
0 


and subsequently, upon second integration, 


o-f exp| (1 a) f saa |-ds. (15) 
n) 


z 


Here A is an integration constant to be 
calculated from boundary condition (11); the 
second boundary condition (12) is automatically 
satisfied by the upper limit of the first integral 
of Eq. (15). Introduction of Eq. (15) into (11) 


gives 
2 Zz (16) 
Io= (1/a) nde+A |-dr. 
Lolo 


On division of Eqs. (14) and (15) by (16), the 
integration constant A drops out. Thus 
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in exp| (1/0) f ods| 
== — 8 (17) 
x i 2) z 
exp} (1/a) ae|-as 
Jf exs[are > 
J exp] (1/a) f ede |-as 
9 = Io. ——— (18) 











J exp| (1/a) f vas! * 


On the basis of Eq. (17), the heat transfer 
coefficient h, defined as the quotient of heat flux 
by temperature difference #o, is 


1 od 
HG). 
do Ox z=0 


kb . 
= — ——_—. (19) 


f exo| (1/0) f vate |-dx 


0 








Dice comsienicenin 


r) a 1 zr 
fom Jf exo (=): f vae]-ax | -as 


The integral involved in Eq. (19). where v, 
is defined by Eqs. (4) and (5), is to be calculated 
by numerical integration. Thus 

h=0.339-k- (w/v)? (20) 
for v/a=0.74, valid for air. 

In analogy to von K4rmA4n’s calculation of the 
torque momentum, the heat transfer coefficient 
can also be calculated by means of the following 
consideration. The enthalpy flow from the disk 
of surface rr’? to the boundary layer must equal 
the enthalpy flow in the boundary layer in the 


radial direction at the distance r from the axis. 
Thus 


5 
rrhoa= f pu,’ 0Cp* 2ar-dx (21) 
0 


and consequently, on division by zr’d> and 
introduction of Eqs. (18) and (10), 





2 


(22) 


; f exp| (1/a)- f vas |-as 
0 0 


where the velocity components v, and v, are 
given by Eqs. (1) and (4). Numerical evaluation 
yields 


h=0.335-k- (w/v)! (23) 


for v/a=0.74, valid for air. 

The method employed for the result in Eq. 
(23) is to be regarded as a closer approximation 
than that employed for Eq. (20), for representa- 
tions of velocity components in Eqs. (1) and (2) 
are only approximations, and small deviations 
from the true velocity functions cause a greater 
error in differentiation of the function #(x) as 
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applied in (19) than in integration as applied in 
(21). However, the difference between the re- 
sulting Eqs. (20) and (23) is insignificant. 
Equation (23) is expected to hold true as long 
as the flow in the boundary layer is essentially 
laminar. According to friction measurements by 
W. Schmidt,? this condition is met if the Rey- 
nolds number, 
Re=(Arr’w/v), 


is less than 500,000. f 


~ 


(24) 


2W. Schmidt, Zeits. d. Vereins Deuts. Ing. 65, 441 
(1921). 
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An Interferometric Method for the Calibration of Electron 
Microscope Magnification 


J. L. FARRANT AND A. J. HODGE 
Division of Industrial Chemistry, Council for Scientific and Industrial Research, Melbourne, Australia 
(Received February 6, 1948) 


Interferometric methods have been used to measure uniform cylindrical silica and glass fibers 
as small as 5u in diameter with better than 1 percent accuracy. Such fibers yield sharp images 
when examined in the electron microscope and constitute a convenient and rugged specimen 
for accurate direct calibration of the instrument for magnifications up to 10,000 diameters. 
Higher magnification ranges are then readily determined by means of an internal calibration 
obtained by taking micrographs of glass spheres about 0.5 micron in diameter. Such spheres 
are also convenient specimens for determining the image distortions which must be known to 
correct the widths of the fiber images. This method appears to be capable of greater accuracy 


than any employed hitherto. 


INTRODUCTION 


HE increasing use being made of the electron 
microscope for the measurement of long 
spacings in fibrous materials'? and of the dimen- 
sions of macromolecules* * and the desirability of 
correlating these results with those obtained by 
x-ray diffraction and ultra-centrifuge methods, 
respectively, have created a demand for accurate 
means of electron microscope calibration. 

Four different calibration methods have been 
described in the literature. Borries and Ruska® 
used a “two-step” procedure in which the mag- 
nifications of the first and second stages are 
determined separately by imaging suitable aper- 
tures with each stage individually. The aperture 
diameters are large enough to be measured 











Fic. 1. Fiber No. 10 between a cover glass and an optical 
flat illuminated with the Hg \5461 line. 





1C, E. Hall, M. A. Jakus, and F. O. Schmitt, J. Am. 
Chem. Soc. 64, 1234 (1942); J. App. Phys. 16, 459 (1945). 
2 J. L. Farrant, E. H. Mercer, and A. L. G. Rees, Nature 
159, 535 (1947). 
3W. C. Price, R. C. Williams, and R. W. G. Wyckoff, 
Arch. Biochem. 9, 175 (1946). 
( *W. C. Price and R. W. G. Wyckoff, Nature 157, 764 
1946). 
5 B. v. Borries and E. Ruska, Naturwiss. 27, 577 (1939). 
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accurately with a light microscope. Another 
“‘two-step’” method developed by the 
Columbian Carbon Company® in which a speci- 
men is moved through known distances by means 
of a micrometer screw. A number of workers have 
used a direct comparison method’ in which 
micrographs of the same specimen are taken with 
both the electron microscope and a calibrated 
light microscope. Burton, Barnes, and Rochow’ 
have developed a method in which a thin 
replica of a diffraction grating of known perio- 
dicity is used as an electron microscope specimen. 

The use of the two-step methods mentioned 
has lapsed as the errors involved amount to 
about 20 percent. 

The direct comparison method is useful, par- 
ticularly at magnifications up to 1000 diameters, 
but is limited in accuracy for two reasons. The 
light microscope when used without oil immer- 
sion has a limit of resolution of 0.25 micron cor- 
responding to a useful magnification of 750 
diameters, so that a single particle must be at 
least 25 microns in diameter to enable its size to 
be measured with an accuracy better than 1 per- 
cent by means of a light microscope. Such a 
particle fills the 5-cm field of the RCA electron 
microscopes at a magnification of 2000 diam- 
eters; so the higher ranges up to 20,000 diameters 
must be calibrated by means of internally deter- 


was 


® Columbian Carbon Company: ‘The particle size and 
shape of colloidal carbon as revealed by the electron 
microscope,” Binney and Smith Company, Distributors, 
New York (1940). 

7C. J. Burton, R. B. Barnes, and T. G. Rochow, Ind. 
Eng. Chem. 34, 1429 (1942). 
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mined ratios of the electron microscope magni- 
fications. This difficulty can be avoided by using 
a specimen consisting of individual isometric 
particles about 0.2 micron in diameter, separated 
by distances of the order of 5-10 microns and 
supported by a thin membrane. The second ob- 
jection, that the supporting membrane may 
temporarily stretch or contract while under the 
influence of the electron beam, is then en- 
countered. 

The fourth method is that most generally used. 
It has been found’ that the periodicity of the 
replica is rarely the same as that of the grating 
from which it is prepared. It is usual to attempt 
to avoid this by using either silica replicas pre- 
pared by the method of Heidenreich and Peck® 
or metal replicas as suggested by Wyckoff and 
Williams.*® The periodicity of the replica may be 
checked by mounting it on a spectrometer and 
determining the angles through which it diffracts 
light of known wave-lengths. This procedure, 
however, cannot be applied to determine any 
permanent or temporary local irregularity of the 
spacings Over areas as small as those encountered 
in the electron microscope at high magnifica- 
tions. Such permanent irregularities may appear 
owing to local strains during preparation of the 
replica or as a result of variations or periodic 
errors in the rulings. We have noticed temporary 
variations in the spacings of the replicas while 
under the influence of the electron beam in the 
case of collodion and Formvar replicas even when 
“‘shadow cast’’ with up to 30A thickness of gold. 

An alternative to these methods would be to 
employ as a calibration specimen a uniform 
cylindrical fiber a few microns in diameter. Such 
fibers are easily drawn from glass or silica. Fine 
quality glass cloth is a convenient source of 5—6 
micron fibers. 

However, the accurate determination of the 
diameter of fibers less than 10 microns in diam- 
eter presents some difficulty. Burton, Barnes, and 
Rochow’ suggested measuring the diameter by 
weighing on a microbalance, but rejected. the 
proposal because fibers cannot be obtained uni- 


®§R. D. Heidenreich and V. G. Peck, J. App. Phys. 14, 
23 (1943). 

*R. C. Williams and R. W. G. Wyckoff, J. App. Phys. 
17, 23 (1946). 
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Fic. 2. Microphotometer trace across negative of Fig. 1. 


form over the considerable length required to 
give sufficient weight. 

We have found it possible to determine the 
diameters of 5 micron fibers as short as 3 mm 
to an accuracy better than 0.3 percent by the use 
of an interferometric method. Once such a 
measured fiber is mounted on a specimen screen 
it is available as a permanent direct calibration 
specimen suitable for magnifications up to 10,000 
diameters. Magnifications up to the maximum 
value of 20,000 diameters of commercial instru- 
ments can be determined by the use of a 1-2- 
micron fiber or glass sphere calibrated by direct 
comparison with the interferometrically meas- 
ured fiber. This small ratio of 2 to 3 may be 
readily and accurately determined provided 
allowance is made for image distortion. 

It has been stated’ that sufficient charge may 
accumulate on the silica fibers to cause variations 
in the image size when the objective current is 
varied during focusing. This has not been our 
experience, possibly because our micrographs 
have been taken with the condenser lens current 











Fic. 3. Fiber No. 10 between a convex lens and an optical 
flat illuminated with the Hg 5461 line. 
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Taste I. 











Image Image 
width width 
Diam- (uncor- (cor- 
Fiber Illumi- eter rected) rected) Magni- 
Method No. nation “ mm mm fication 
Cover glass 
and flat 1A* Na D 6.77 44.4 41.0 6060 
Cover glass 
and tlat 2A NaD 6.18 39.3 36.8 5960 
Cover glass 
and flat 3A Na D 5.89 38.0 35.7 6060 
Cover glass 
and flat 5A NaD 5.96 38.4 36.1 6060 
Cover glass 
and flat 9B NaD 4.54 28.1 27.0 5950 
Cover glass 
and flat 10B Na D 5.91 
Cover —. 37.4 35.2 5960 
and flat 10B Hg AS461 5.91 
Convex lens 
and flat 10B NaD 5.90) 
Convex lens > 7.4 35.2 5980 
and flat 10B Hg S461 5.88) 








* An interval of two months elapsed between the measurements of 
series A and series B. Visual estimates only of the band order were 
made in series A. 


adjusted so that a considerable area of the wire 
mesh supporting the fiber is bombarded by the 
electron beam. This is a well-known method! of 
avoiding charging of the specimen. 


I. THE MEASUREMENT OF FIBERS 


The procedure in the methods we have used 
to measure fiber diameters involves placing the 
fiber between two flat glass surfaces so that it 
determines their separation at the fiber, and then 
measuring this separation interferometrically. 
The system is illuminated by well collimated 
monochromatic light, and the interference fringes 
are viewed by reflection with a low power micro- 
scope. The diameter d of the fiber is given by the 
formula d=n(\/2), where n is the interference 
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Fic. 4. Microphotometer trace across negative of Fig. 3. 


10 J. Hillier and E. G. Ramberg, J. App. Phys. 18, 48 
(1947). 
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band order at the fiber and ) is the light wave- 
length. The band order at the fiber is obtained by 
counting the bands back to the region of zero 
order where the two glass surfaces are in contact. 
Fractions of a band order are determined by use 
of a micrometer eyepiece or by photography and 
subsequent measurement either directly or from 
a microphotometer trace. 

In practice it is found difficult to locate a 
region of zero band order when thick optical 
flats are used, as the areas in contact are ex- 
tremely small. This difficulty is easily overcome 
by replacing the upper optical flat with a thin 
microscope cover glass. Slight pressure applied a 
few millimeters to one side of the fiber brings the 
glass surfaces into contact over a considerable 
area. Figure 1 is a low power photomicrograph of 
fiber 10 (Table I) and the associated fringes ob- 
tained when such a system was illuminated by 
means of the Hg \5461 line. The negative was a 
Kodak Super XX plate; the print was made on 
normal contrast paper and normal contrast de- 
velopment conditions were maintained. A micro- 
photometer trace taken across the fiber shown 
in the negative of Fig. 1 is reproduced in Fig. 2. 
This curve yields 21.66 as the interference band 
order at the fiber which gives 5.91 microns as the 
fiber diameter. When the same system was 
similarly examined by means of the Na D line, 
the interference band order was 20.06, also cor- 
responding to a fiber diameter of 5.91 microns. 

An alternative method is to use a plano-convex 
lens in place of the cover glass mentioned above. 
The fiber is placed upon a glass optical flat with 
the lens adjusted so that its spherical surface 
touches both the flat and the fiber, the fringes are 
photographed, and a microphotometer trace of 
the negative is obtained. Figure 3 shows the 
Hg 5461 fringes obtained in this manner with 
fiber 10. From the corresponding microphotom- 
eter trace shown in Fig. 4, the interference band 
order at the fiber is given as 21.54 and the fiber 
diameter as 5.88 microns. Use of the Na D line 
gave a band order of 20.03 and 5.90 microns as 
the fiber diameter. These results show that it is 
possible to determine the fiber position to within 
0.1 of a fringe. When the Hg 5461 line is em- 
ployed the band order at a 5.5 micron fiber is 
about 20, so accuracy of 0.5 percent is obtainable. 
This is true only if the presence of dust on the 
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fiber and glass surfaces is scrupulously avoided. 
Because of the curvature of the cover glass or 
lens, the interference bands are not equidistant, 
but conditions may be readily arranged to ensure 
that the curvature at the fiber is insufficient to 
introduce an appreciable error if the fraction of 
the band order is determined by direct inter- 
polation. 

Tolansky" has found that the sharpness of 
interference fringes can be greatly increased by 
evaporating about 500A thickness of silver onto 
the glass surfaces. Multiple beam interference is 
then obtained owing to the high reflectivity of 
the films, and the resultant sharpening of the 
fringes is sufficient to allow the fiber to be 
located to within 0.02 of a band. Our attempts to 
make use of this fact were eventually abandoned 
for lack of sufficiently accurate means of deter- 
mining the phase change occurring when light is 
reflected from silver films of this nature. To 
avoid the necessity of an accurate knowledge of 
the phase change, another method was tried. 
The fiber was held in contact with an optical 
flat by means of a collodion film picked up from 
a water surface, a silver film was evaporated onto 
the collodion, and observations were then made 
on the interference fringe system formed on 
illuminating the wedge between this surface and 
another silvered flat. The maximum displace- 
ment of the fringes as they crossed the fiber 
should then have given the fiber diameter to a 
high degree of accuracy. Attempts to observe 
the fringes as they crossed the fiber were unsuc- 
cessful owing to the small size of the fiber and the 
low light intensities available from mono- 
chromatic sources. 

After measurement, the application of a little 
collodion cement to each end of the fiber serves 
to hold it permanently across a standard wire 
mesh specimen screen. 


II. CORRECTION OF FIBER IMAGES FOR 
DISTORTION 


As shown in Fig. 5, the images of the fibers are 
extremely sharp but are considerably distorted. 
This is largely due to the aberrations of the 
projection lens; the matter has been discussed at 
length by Hillier.” 


1S. Tolansky, Proc. Roy. Soc. A184, 41 (1945). 
2 J. Hillier, J. App. Phys. 17, 411 (1946). 
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Fic. 5. Electron micrograph of Fiber 9 showing sharp but 
distorted image. 


One effect of this distortion is, in the case of a 
small sphere measured radially, to increase the 
magnification by a factor f(r) as it is moved 
outwards from the image center by an image dis- 
tance r. Plots (Fig. 6) of f(r) against r are readily 
obtained from a series of exposures, taken on one 
plate, of a glass sphere as it is moved progres- 
sively from the field center to one corner. In 
order to use such a curve to correct the relatively 
wide fiber images some convenient method of 
carrying out the integration involved is required. 
This is apparent from the fact that an image 
radial distance r measured from the image center 
would be represented by the distance 


’ dr 
= 70 
in the absence of distortion. The function 1/f(r) 
is well represented by the expression 
1—Ar—Br’, 


and the constants A and B may be determined by 
curve fitting to the plot of 1/f(r). In this way r’ 
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Fic. 6. Radial distortion factor f(r) plotted against image 
distance r from center of field. 


843 








is readily calculated from the equation 


r’ =r(1—}Ar—4Br). 


The importance of this correction is shown by 
the fact that at a magnification of 6000 diameters 
the correction to a 5.9-micron fiber amounts to 
6.0 percent in the case of our RCA Model U 
electron microscope. This instrument has been 
fitted with the high intensity electron gun recom- 
mended by Hillier," and caution is necessary to 
avoid focusing too heavy a beam current onto 
the fibers and so fusing them. Glass spheres for 
the determination of distortion are readily made 
by feeding a coal gas flame with air or oxygen 
containing suspended powdered glass. 


Ill. ACCURACY OF THE METHOD 


In Table I are set out the final values for the 
magnification for one setting of the microscope 
obtained by means of six different fibers measured 
by the two methods described. Fibers 1, 2, 3, 5, 
9, and 10 were measured between a cover slip 


8]. Hillier, J. App. Phys. 17, 307 (1946). 


and an optical flat, while fiber 10 was also 
measured between a convex lens and a flat. The 
deviation from the mean magnification is less 
than 1 percent. 

The good agreement of the various results 
shows that the fibers are uniform and cylindrical 
to a high degree and that neither they nor the 
glass flats are seriously distorted by the pressure 
applied during measurement provided reasonable 
care is taken. It is important to realize that good 
reproducibility of the magnification is obtainable 
only if the battery potential, which controls the 
electron accelerating potential, is maintained 
constant, if hysteresis effects in the lenses are 
avoided by switching the lens currents off and on 
several times before making an exposure, and if 
the specimen is always supported at the same 
height. 
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Electronic Radiography and Microradiography* 
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Two techniques of radiography by means of secondary photoelectrons liberated by hard 
x-rays (180 kv), one by reflection and one by transmission, are described. Special fine-grained 
photographic emulsions of the Lippmann type are required which will not be affected by x-rays 
traversing the film, but are sensitive to photoelectrons liberated from the surface of a specimen 
in intimate contact with the photographic film or paper. In this way photographic negatives 
may be reproduced (from the silver of the image), microradiographs prepared by enlargement 
of any area of the radiograph, surface topography studied, ores analyzed, etc. For electronic 
radiography by transmission, instead of the specimen being its own electron generator, a thin 
metal foil (usually lead) pressed against the specimen generates the electrons which then 
traverse the thin specimens, and register on the photographic film or paper. Applications are 
papers, insect wings, varnish or rubber films, plant leaves, tissue sections, plastics, and alloys. 


HE two methods described herewith are 
totally different from the usual radio- 
graphic processes because they make use of 








* Edited, submitted for publication, and proofread by G. 
L. Clark, University of Illinois, Urbana, Illinois. 
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secondary photoelectrons instead of x-rays as a 
means of registration of images on the photo- 
graphic emulsion. We shall call these methods 
electronic radiography by reflection and elec- 
tronic radiography by transmission, respectively. 
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I. ELECTRONIC RADIOGRAPHY BY REFLECTION 
Principle 

In April 1940, in a paper published in the 
Revue Scientifique,! it was suggested that a 
method designed to obtain radiographs or micro- 
radiographs from finished articles without having 
recourse to a destructive method of sample 
preparation should prove useful and interesting. 

The method developed to achieve this purpose, 
designated ‘‘method by reflection,’”’ is based on 
the following principle. : 

Let us consider a plane material surface 
(having any dimension) of an article in which 
components have quite different atomic numbers. 
When irradiated by x-rays having adequate 
wave-lengths, the surface becomes the source of 
a secondary scattered fluorescent wave emission 
and simultaneously of a photoelectron emission 
whose kinetic energy, according to Einstein’s 
and M. de Broglie’s equation, depends upon the 
incident quantum and the extraction work of 
the electrons, that is to say, upon the atomic 
number: 

mv = hv(incident) — W,,, 


where W,, is the extraction work corresponding 
to an m level; it varies proportionately with the 
atomic number. 

The number of emitted photoelectrons depends 
upon the wave-length and the intensity of the 
incident x-rays. It increases when the atomic 
number increases, heavy atoms actually having 
more electrons than light ones. Thus, use is 
made of the photographic effect of the secondary 
electronic emission which is closely related to the 
fluorescent phenomenon. It should be noticed 
that this is merely a surface effect, electrons 
being much more absorbed than x-rays in a 
material substance.?~> 

In order to obtain the desired result, it is 
necessary that the emulsion which is in closest 


an J. Trillat, ‘‘La microradiographie,” Rev. Sci. 4, 212 
(1946). 

3 J. J. Trillat, ‘Sur l’action photographique des electrons 
secondaires resultant de Il’action des rayons x sur, les 
metaux,’’ Comptes Rendus 216, 179 (1943). 

4J. J. Trillat, “Sur une extension de la methode de 
ma par reflexion,”” Comptes Rendus 216, 683 
1943) 

* J. J. Trillat, “‘Radiographie et Microradiographie par 
transmission et par reflexion,” Bin Ste Fr. Electriciens 6, 
III, 1-20 (1943); and Rev. ‘Sci. 8, 364, 370 (Sept., Oct. 
1942). 
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Fic. 1. Diagram of the method of electronic radiography by 
reflection. 


possible contact with the surface to be investi- 
gated and through which the incident x-rays 
will pass first should be back impressed by the 
photoelectrons only and not by the incident or 
simply scattered x-rays. It is necessary that the 
absorption coefficient of the photographic emul- 
sion should be very low for this particular 
radiation and, on the other hand, that the 
photographic emulsion should be impressed by 
the photoelectrons in a sufficiently strong way. 

In order that the absorption coefficient of the 
emulsion for x-rays should be low, it is necessary 
that the photographic layer be very thin; it must 
consist in very fine and non-sensitive grains. 
Furthermore, the x-rays used must be hard, that 
is to say, they must be emitted by a tube 
operating at high voltage ; thus the absorption in 
the emulsion and therefore the photographic 
effect are lessened. 

The first condition is realized when one uses 
special Lippmann films which are used in micro- 
radiography, the thickness of the sensitive layer 
being approximately 1/300 mm. Excellent results 
have been obtained with Lippmann films pre- 
pared by Eastman Kodak Company. Thus, 
highly contrasting emulsions may be employed 
and even silver chloride photographic papers 
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Fic. 3. 


(Velox) which are only slightly sensitive to 
x-rays but very sensitive to electrons. 

The second condition is satisfied when using 
constant voltages higher than 100 kv and up to 
200 kv. 

It may be seen that this method is quite 
different from ordinary microradiography or 
radiography by transmission, as illustrated by 
Fig. 1, in the following respects: 

(a) The photographic film is first traversed by 
incident x-rays (which looks like a paradox). 

(b) X-rays are not active, but photoelectrons 
are. 

(c) Instead of using soft x-rays as is generally 
the case with ordinary microradiographs, hard 
x-rays are used here. 
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(d) Instead of investigating the whole of a 
radiated mass which results in the appearance of 
superposed layers, the structure of a purely 
superficial film is investigated. 

The technical procedure is very simple; it 
comprises applying the film very closely against 
the surface to be investigated and exposing the 
whole system to x-rays. It is essential that the 
film should be in very close contact with the 
surface because photoelectrons are emitted in all 
directions and are greatly scattered in air. 


Results 
(a) Investigation of metallic surfaces 


Let us place in an adequate holder (Fig. 2) a 
Lippmann film in contact with a plane surface 
consisting of metals having different atomic 
numbers (26 Fe—28 Ni—29 Cu—50 Sn—82 Pb) 
and let us expose the surface to x-rays which are 
emitted by a tube operating under 180 kv; the 
soft part of the continuous spectrum is elimi- 
nated by means of a 4-mm copper filter. Under 
these conditions, a picture is obtained (Fig. 3) 
where the photographic density measured with a 
Capstaff-Purdey densitometer is exactly propor- 
tional to the atomic number as shown in Table I. 

Thus, elements of consecutive atomic number 
such as Fe, Ni, and Cu, are separated easily and 
one may foresee interesting applications of this 
fact in the field of alloys and mineral studies. 

Provided a certain care is taken, it is even 
possible to use ordinary commercial fine-grained 


Fic. 4. Comparisor. of an ordi- 
nary radiograph of an alloy of 
aluminum, antimony, tin, and 
lead and an electronic radio- 
graph by reflection of the same 
alloy (Saulnier). The ordinary 
radiograph is made at 50 kv, the 
electronic at 180 kv; filtered 
2-mm complex and 2-mm alu- 
minum; exposure time 4 min.; 
distance 45 cm. 
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Fic. 5. Uranium mineral, enlargement 6 diameters. 


emulsions, such as ‘‘Panatomic’’ Kodak emul- 
sion. Exposure times are greatly reduced and 
the process requires a few seconds only; when 
such films are adopted, the field of application 
of the method is extended; true radiographs by 
reflection are made and then magnified 20 or 30 
times, later on. 

Systematic tests have shown that the optimum 
practical range of voltages for this method is 
150-200 kv; in this range the effects of fluorescent 
x-rays are largely negligible in comparison with 
the photoelectron effect. Ordinary microradi- 
ography by transmission would not have given 
the same results without using monochromatic 
x-rays and selecting the radiations so that their 
wave-lengths should fit inside the critical absorp- 
tion discontinuities of the elements to be 
differentiated. 


(b) Alloy inspection 


When using Lippmann films, one may convert 
the radiographs by reflection into as many 
microradiographs as required ; in fact, each point 
of the film might be magnified many thousands 
of times because of the excessive fineness of the 
grain; thus beside the whole image, any required 
detail photographs are obtained. Thus, a single 
process is sufficient to yield simultaneous infor- 











TABLE I. 
Fe N=26 Photographic density 0.83 
Ni 28 0.85 
Cu 29 0.86 
Sn 50 0.94 
Pb 82 1.10 
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Fic. 6. Pitchblende, black uranium, gray lead, white SiOz, 
enlargement 6 diameters. 


mation about the macro- and the microstructure 
of a specimen. 

Figure 4 illustrates a radiograph by reflection 
of an Al alloy containing Pb, Sb, and Sn. 

In case macrographic investigation only is 
aimed at (for instance, for segregation or hetero- 
geneity investigations), it is not necessary to use 
a special very fine-grained film, since the image 
obtained requires no magnification. Excellent 
results were obtained by simply placing a rela- 
tively insensitive photographic silver chloride 
paper, for instance a ‘Velox contrast’’ paper, in 
contact with the surface. All manipulations may 
then be made without resorting to complete 
darkness, and exposure time is greatly reduced. 


(c) Reproduction of photographic negative 


The silver image of a film or plate negative 
might be used as a secondary electron source 





Fic. 7. Galena containing nickel arsenide and SiQs, enlarge- 
ment 6 diameters. 
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which will register an image on a fresh photo- 
graphic paper of film.? A negative to be repro- 
duced is placed at the bottom of an ordinary 
photographic frame in which the glass has been 
replaced by a black celluloid sheet; next to the 
negative is placed a fresh photographic paper or 
film with its sensitive surface in contact with the 
silver image. The whole system is exposed to 
hard x-rays which first go through the celluloid 
sheet, then through the paper, and finally strike 
the negative plate. Photoelectrons emitted by 
the silver of the image will then impress the 
photographic paper which on development will 
yield the exact image of the negative. This 
reproduction process gives a direct copy of the 
original without any negative intermediate. 


(d) Surface topography investigation 


The reflection method enables us to get some 
information concerning the mechanical condition 
of a surface. When a surface has been mechani- 
cally worked and when superficial heterogeneities 
exist, these defects will appear prominently on 
the exposed negative. The process is reliable only 
when the path in air is not too long, on account 
of electron scattering. 


(e) Study of ores 


Lastly, we have recently applied this method 
to the study of polished surfaces of various ores 





Fic. 8. Electronic radiograph by transmission. Butterfly 
wing, enlargement 8 diameters. 
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containing elements with different atomic num- 
bers. The results were excellent, and it is to be 
supposed that this process will be of use to 
geologists and mineralogists.® 

For example, Figs. 5 and 6 represent a piece of 
uranium ore (pitchblend) in which the uranium 
present is clearly outlined in black. The bright 
parts correspond to the lighter elements such as 
silicon. 

Figure 7 represents ore containing plumbago 
(PbS), together with nickel arsenide (NiAs2). 
The whites correspond to the gangue of silica 
(SiOz), the blacks to the plumbago, the grays to 


*-NiAse. This method can even be used to find out 


the quality of the ore under examination, by 
measuring the blackening of the plate and 
comparing with a standard one. 


Limitations of the Electronic Radiography by the 
Reflection Method 


Limitations of this method depend on the 
electron diffusion (scattering) almost solely, 
since electrons are emitted in any direction by 
the surface layer; however, the inconvenience is 
not so great as when one uses the fluorescent 
secondary x-rays, as M. Guinier did.’ Neverthe- 
less, it is possible to minimize the scattering 
nuisance by applying the photographic film or 
paper as perfectly as possible against the investi- 
gated surface and by using a vacuum apparatus 
when necessary. In any case, it is absolutely 
essential to have a nearly perfect contact be- 
tween the object and the film, because of the 
great absorption and diffusion of electrons in air. 


Discussion of Results 


It was explained above that secondary elec- 
trons are active in our experiments and fluores- 
cent x-rays are not. Scattered x-rays, both 
unmodified and modified by the Compton effect, 
do not appreciably affect the film. On the con- 
trary, secondary fluorescent x-rays which are 
emitted upon impact with primary quanta—K, 
L, etc., rays—may possibly have fairly long 
wave-lengths, especially when light elements are 
involved, and thus might well play a part in the 

‘J. J. Trillat and Ch. Legrand, “Radiographie elec- 
tronique de minerals,’’ Comptes Rendus 224, 1000 (1947). 


7A. Guinier and J. Devaux, ‘‘Sur la radiographie par 
reflexion,’’ Comptes Rendus 214, 223 (1942). 
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investigated effect on account of their rather 
high absorption. 

Within the experimental conditions, character- 
ized by voltages higher than 100 kv, we may 
safely conclude that the film is affected by 
photoelectrons only. This may be checked by a 
very simple experiment. With a surface con- 
sisting of aluminum, iron, or lead, an exposure 
is made first by reflection of electrons generated 
below 180 kv; then a second exposure is made 
with conditions being exactly the same except 
that a very thin celluloid foil (1/20 mm) is 
placed in the path of the x-ray beam. In the 
last case, no image will appear, a proof that we 
are actually concerned only with an electronic 
effect. If characteristic x-rays, such as copper, 
iron, or lead K-rays which are strongly excited in 
these conditions, were involved in the effect, the 
low thickness of the celluloid foil (1/20 mm) 
would not stop them, and a blackening of the 
film almost as severe as the one observed in the 
first experiment would be observed. 

It may be interesting to compare these results 
with those obtained later by M. Guinier,* who 
made use of long wave-length incident x-rays for 
his researches instead of using short wave-length 
X-rays, as was done in this work. 

M. Guinier has shown that when illuminating 
an iron plate with rays having a wave-length 
below that of the K absorption discontinuity of 
iron (Ax = 1.74A), an intense blacking will appear 
on the photograph; it is due to fluorescent K-rays 
of iron. But when wave-lengths are shorter, for 
instance, about 0.25A (that is to say, the voltage 
is 50,000 v) there is a change in the phenomenon 
and the effect is mainly due to photoelectrons; 
this is checked by interposing a 2/100-mm thick 
Cellophane foil in the path of the rays. Our own 
experiments which were carried on partly at 
70,000 to 100,000 v demonstrated that photo- 
electrons only were decidedly active. 

Thus, the observed effect represents the com- 
bination of two other effects, the relative magni- 
tude of which depends on the incident radiation 
wave-length; when taking into account the 
number of photoelectrons, their energy, and 
decreasing sensitivity of the photographic emul- 
sion as wave-length decreases, it will be found 


8A Guinier, “Etude theorique de la radiographie par 
reflexion,’’ Revue Sci. 63, 52 (1943). 
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Fic. 9, Electronic radiograph by transmission of a sheet of 
cigarette paper, enlargement 8 diameters. 


easily that the combined effect is due solely to 
fluorescent x-radiation when \<Ax. This effect 
remains practically constant for medium wave- 
lengths, and then increases when wave-lengths 
become shorter. Then, electrons only act on the 
emulsion ; finally, the effect decreases when rays 
of high penetrative power (above 200 kv) are 
used. M. Guinier’s experiments constitute a 
favorable complement to our own, in the long 
wave-length field. Analytical applications may 
be foreseen for the method of using the fluores- 
cence radiation of illuminated elements instead 
of photoelectrons. 


II. ELECTRONIC RADIOGRAPHY BY 
TRANSMISSION 


Working with M. A. Saulnier,® the author has 
also attempted to use the photoelectric effect in 
order to obtain radiographs from a very thin 





Fic. 10. Electronic radiograph of a section of rabbit kidney. 
* A. Saulnier and J. J. Trillat, “‘Radiographie par elec- 


trons secondaires,’"’ Comptes Rendus 220, 772 (1945); and 
Revue Sci. 3244, 211 (1945). 
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Fic. 11. Electronic radiography by transmission of an 
onion peeling. 


object; instead of being its own electron gener- 
ator, as in the radiography by reflection, the 
article is traversed by photoelectrons, extracted 
by means of x-rays from a metallic surface 
against which the article is strongly pressed. 

To obtain the required result, it is necessary 
that the photoelectrons should be emitted in 
great quantity and that they should have suff- 
cient energy to go through a preparation which 
may be a few thousandths of a millimeter thick; 
these conditions are practically obtained when a 
lead radiator is used and when a 150- to 200-kv 
primary quantum is used. It is, of course, neces- 
sary for the thin sheet or film to have a constant 
thickness, any local thickness variation being 
made obvious by an additional absorption of 
the photoelectrons. 

Incident x-rays produced in a tube at 150-200 
kv first go through a blackened celluloid plate 
which acts to prevent light from striking the 
film, and for maintaining the investigated prepa- 
ration. Then, radiation is passed through a 
perfectly smooth and plane 0.2-mm thick lead 
foil which filters out the soft part of the rays 





Fic. 12. Bubbles of air in a film of paint of different 
thickness. The bubbles appear as the black dots, enlarge- 
ment 3.5 diameters. 
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and which emits from its back surface a large 
quantity of photoelectrons, thus playing the 
part of a radiation-transmitting unit. The 
preparation to be inspected is tightly held be- 
tween the radiating source (lead foil) and a 
photographic film or paper. Under these condi- 
tions, the photographic emulsion is affected by 
the photoelectrons which passed through the 
preparation and also by those which are emitted 
by the preparation proper. Incident or fluores- 
cent x-rays do not produce any noticeable image, 
as shown by our previous researches, provided a 
tube operating between 150 and 200 kv is used. 

In this way a true electron radiography is 
obtained, and it may be compared with the 
results given by an electron microscope, if magni- 
fying power is overlooked. When an extremely 
fine-grained emulsion is used, such as the one 
used in microradiography, the resulting image 
might be magnified and many details which did 
not appear on the original print are made clear. 
Exposure time varies from a few seconds to 


. several minutes. 


This method has been applied successfully on 
papers, insect wings, varnish or rubber films, 
plant leaves, plant or animal tissue sections, 
plastics, alloys, etc. Very clear pictures are often 
obtained and give added information to the 
microscopi¢ inspection; this result is interesting 
in particular when opaque specimens which can- 
not be inspected under a microscope are con- 
cerned. A few examples are given in Figs. 8-11. 
I have also applied this method to the study of 
thin oil films which are in contact with pieces of 
machinery. It is thus possible to disclose the 
structure of films about 2y thick. Further study 
of some lubrication phenomena and _ paints!® 
(Fig. 12) is in progress. It is quite obvious that 
the examples could be multiplied and the method 
applied to numerous cases. It is a very simple 
method which requires only ordinary medical 
radiographic apparatus, provided it may be 
operated above 150 kv. The x-ray apparatus 
need not even be kept at a constant voltage; a 
pulsating voltage is sufficient provided the 
radiation is filtered through a few millimeters 

© J. J. Trillat and Ch. Legrand, “Sur une application de 
la Radiographie electronique a la detection de couches 
minces organiques ou minerales,’’ Comptes Rendus 224, 


“a9 alee Rev. de |'Institut Francais du Petrole 2, 250 
1947). 
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copper so as to eliminate the soft part of the 
spectrum. 

We must point out that a quite similar method 
was lately utilized in Great Britain by H. S. 
Tasker and S. W. Towers;"! Our own researches 
had been initiated for a long time when we heard 
of the work which was performed quite inde- 
pendently. The English authors obtained results 
which quite agree with our own. 


Discussion of Results 


Electronic radiography by transmission is 
based on the absorption of highly energetic 
photoelectrons by matter; it is thus essentially 
different from soft x-ray radiography such as has 
been utilized by Dauvillier! and Lamarque.” 
The photoelectron absorption law is very differ- 
ent from that of x-ray absorption and is less well 
known; while x-ray absorption varies propor- 
tionally to the fourth power of the atomic num- 
ber, the electron permeability of matter which 
depends on the density has no simple relation to 
the atomic number. It varies little from one 
element to another when compared with the 
coefficient of absorption of x-rays; thus electron 
absorption in gold is nearly three times that in 
aluminum, while the x-ray absorption coefficient 
is twenty times greater in gold. ‘‘Electronic”’ 
negatives should then be explained in a different 
way than x-ray negatives; exactly the same 
problem arises in the explanation of images ob- 
tained with the electron microscope by means of 
electrons having essentially the same energy as 
that of the photoelectrons in our researches 
which pass through a very thin specimen. 

As might be expected, the main limitation to 
the method is the thickness of the specimen; 
it must be very thin (a few to several thousandths 
of a mm) and the thickness must be constant. 
A. Saulnier made measurements to show that 
the absorption of photoelectrons by aluminum 
or Cellophane screens does not depend on the 
quality of the absorbant matter, but in a rough 
approximation on its density only. 


4H. S. Tasker and S. W. Towers, Nature 156, 50 (1945). 

2A. Dauvillier, Comptes Rendus 191, 1287 (1930); J. 
Phys. 6, 229 (1932). 

%P, Lamarque, “L’Historadiographie, Union Intern. 
contre le Cancer,” 4, 205 (1939), 

4 A, Saulnier, “Adsorption des electrons secondaires par 
des ecrans minces,’”” Comptes Rendus 220, 772 (1945). 
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Moreover, photoelectrons are emitted in any 
direction and with varying energy from the lead 
surface; thus, the sharpness of the image might 
be altered. 

In order to determine the best conditions for 
applying the method, it is, of course, necessary 
to investigate theoretically the secondary elec- 
tron emission conditions. This research was 
executed by M. A. Saulnier and will be the 
subject of a publication elsewhere. It was found, 
for instance, that electrons ejected backwards 
(Compton effect) need not be taken into account, 
that photoelectron emission is about 1.2 times 
stronger along the direction of the incident 
radiation direction than in the opposite direction 
(Cooksey), and that the asymmetry increases as 
wave-length decreases and atomic weight de- 
creases. The computation shows that lead surface 
layers only (~8/100 mm for LZ electrons) are 
involved in the emission. 

Finally, it is interesting to note that electronic 
radiography is produced from hard x-rays which 
are entirely unsuited for the radiography of such 
thin articles directly, but that the lead screen, 
acting by converting the hard rays into cor- 
puscular rays which are strongly absorbable and 
possess a very high photographic effect, makes 
possible a highly successful technique. 


SUMMARY 


When x-rays hit an element of atomic number 
N, it is a known fact that the kinetic energy of 
the photoelectrons produced depends on the 
incident quantum (by application of Einstein’s 
law). 

These photoelectrons are strongly absorbed 
by the photographic emulsion which they deeply 
affect when their energy is sufficient. 

Under these conditions, in order to study any 
surface made up of elements having different 
atomic numbers (an alloy, for instance), it is 
sufficient to place a film or photographic paper 
in contact with this surface and to irradiate the 
whole with a beam of hard x-rays (150 to 200 
kilovolts) which affect only very slightly the 
silver bromide because of their weak absorption. 
An image of this surface is thus given by elec- 
trons, a picture in which the different elements 
come out in different shades of black; superficial 
irregularities can also be thus revealed. ‘By 
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making use of grainless films, of the Kodak- 
Lippmann type, for instance, each point of the 
photograph can be enlarged and details obtained. 

On the other hand, the photoelectrical emission 
from high atomic number metals, such as lead, 
can also be used for the purpose of examining 
the internal structure of thin films such as paper, 
vegetable or animal tissues, paintings, etc. The 
object under study will then be placed between 





a sheet of lead acting as radiator and a photo- 
graphic film, and the whole is submitted to the 
action of hard x-rays ; the photoelectrons radiated 
by the lead register on the film after passing 
through the object. This last process has been 
used jointly with Dr. Saulnier. 

Examples of applications are given by the 
author together with suggestions of further 
development of these methods. 





y A Temperature-Controlled Frevel Focusing X-Ray Powder Camera 


V. VAND 
Research Department, Lever Brothers & Unilever, Limited, Port Sunlight, Cheshire, England 
(Received February 20, 1948) 


A focusing x-ray powder camera of Frevel type is described, suitable for the study of large 
cells of organic crystals. It uses a convergent x-ray beam collimated by a multiple slit system 


and records transmitted reflections only. 


INTRODUCTION 


GREAT shortening of the time of exposure 
can be obtained by the use of ‘“‘focusing”’ 
methods in x-ray diffraction. The Seeman!-Bohlin? 
arrangement is well known for back-reflection 
diffraction patterns. Frevel* uses a convergent 








Fic. 1. Vertical cross section of the 12.5-cm diameter 
focusing x-ray powder camera. 


1H. Seeman, Ann. d. Physik 59, 455 (1919). 
2H. Bohlin, Ann. d. Physik 61, 421 (1920). 
3L. K. Frevel, Rev. Sci. Inst. 8, 475 (1937). 
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Soller slit for focusing the x-ray beam, the 
geometry of the arrangement being otherwise 
that of Guinier and Johansson,’ which uses a 
bent and ground quartz crystal for focusing. 
Other arrangements are those of Smith.® 

In the Frevel arrangement, the diffraction 
pattern is produced on transmission through a 
thin sheet of the specimen, and it is thus suitable 
for the study of organic compounds, having large 
unit cells. 

In order to fulfill the focusing condition, the 
film, the specimen sheet, and the focus must all 
lie on a common circle. 


GENERAL FEATURES OF THE 12.5-CM 
FOCUSING CAMERA 


In this paper a focusing camera of Frevel 
type is described which uses a convergent beam 
collimated by Soller slits, arranged to converge 
to a linear focus. The advantage of the Guinier’s 
camera of the monochromatic radiation is lost, 
but the camera is faster and can easily be made 
from available materials, which compensates this 
disadvantage. The camera and its components 
are shown in Figs. 1, 2, and 3. 


The camera consists of an adjustable base A, 
*A. Guinier, Comptes Rendus 204, 115 (1937); J. Sci. 
Inst. 22, 139 (1945). 


5 T. Johansson, Zeits. f. Physik 82, 507 (1933). 
6 C. S. Smith, Rev. Sci. Inst. 12, 312 (1941). 
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Fic. 2. The camera assembled. 


holding a bracket B, with a multiple slit system 
S and a temperature controlled curved sheet 
specimen holder H carrying the specimen Q, a 
bracket D with a hollow cylinder C holding the 
film F in position and a beam trap 7. The base 
is provided with three leveling screws; for this 
purpose, we have used the base of a Metrovick 
Universal Camera. 


THE COLLIMATING SYSTEM 


The collimating system S consists of twenty 
slits, made to converge to a focus X at an ap- 
proximate distance of 12.5 cm from their exit 
ends. The slits are constructed from a pack of 
twenty-one strips, spaced from each other by 
pairs of wedge-shaped spacers. The strips used 
were 4.3 cm long, 1 cm wide, and 0.1 mm thick. 
The spacers were each 4.3 cm long and 2 mm 
wide, placed at each side of the strip, which 
leaves 6 mm for the height of each slit. The thick- 
ness of the spacers was 0.075 mm at the exit end 
and 0.135 mm at the entrance end, which gives 
the correct angle of the wedge to produce a 
focus at 12.5 cm from the exit of the slit system. 
The strips were made from razor blades, which 
are sufficiently uniform and flat, and the spacers 
by folding a sheet of aluminum foil 0.025 mm 
thick in a steplike manner. The foils may be 
cemented together with a trace of lacquer and 
strips 2 mm wide may then be cut from the re- 
sulting wedge. During the assembly, the wedges 
and the steel strips may be held in position by 
applying a trace of lacquer to their surfaces. 
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Fic. 3. The camera opened. 


After the assembly, the whole pack is enclosed 
in a lead shield L. The half-width of the beam at 
the focus X obtained by the above arrangement 
is about 0.4 mm. 

As a precaution against fouling the ends of 
the slits S by the beams, the exit ends of the 
steel strips were made slightly thinner by grind- 
ing, serving thus as the side pieces of the ordi- 
nary camera slits. 


THE SPECIMEN MOUNTING 


The specimen Q is placed in a curved specimen 
holder H, which in turn is located by two pins 
on a temperature-controlled bracket E. The 
bracket is hollow and is maintained at any tem- 
perature between 0° and 100°C by water circu- 
lated from a thermostatically controlled bath by 
a pump. The temperature of the specimen is 
not exactly that of the water, so that a calibra- 
tion curve must be prepared; melting points of 
pure fatty acids can conveniently serve as 
calibration temperatures in the range used. 

The specimen holder H is made of a sheet of 
metal 1 mm thick, which is bent to a radius 
6.25 cm and provided with a hole of the same 
diameter as the incoming x-ray beam. If the 
specimen Q is self-supporting, it is fixed or pressed 
into the hole; otherwise the specimen may be 
mounted between two thin Cellophane sheets, 
fixed at each side of the holder. It is important 
to have a specimen of correct thickness. A thick- 
ness equal to or thinner than 1/y, where u is the 
linear absorption coefficient, should be used. For 
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heavier minerals, etc., specimens of a correct 
thickness may be conveniently prepared by 
sedimentation of a suspension of the specimen in 
methyl ether with a drop of collodion added, in 
a Petri dish, at the bottom of which was previ- 
ously placed a Cellophane strip. After the evapo- 
ration of the ether, the strip with the sediment 
attached to it may be peeled off the bottom of 
the dish and transferred to the specimen holder. 
By measuring the area of the dish and by weigh- 
ing the amount of the material added, the ab- 
sorption can be accurately calculated from the 
mass absorption coefficient of the material. 


THE FILM ARRANGEMENT 


The film F is pressed against the outside of 
the cylinder C by means of a strip of velvet V, 
pressed on the film by a rubber band R. The 
rubber band is hooked on a rivet M. 

A standard film 3.5 cm X28 cm is used. No 
paper casette is used, and the light is excluded 
by the velvet V and by a strip of black paper P 
fixed inside the brass cylinder. The paper can 
be, however, replaced by a nickel foil to serve 
as a filter for CuK, radiation. 

The cylinder is provided with windows which 
admit the convergent beam from the outside 
and which expose a suitable strip of film to 
diffracted beams. The center-line of the beam 
does not need to pass through the axis of the 
cylinder, so that the camera can be considerably 
offset and the film placed asymmetrically. This 
extends the range of the angles registered on 
one side of the film up to 20=90°. The pattern 
produced by the arrangement is independent of 
the amount of the offset, so that great accuracy 
is maintained without any need for fine adjust- 
ment of the beam, provided the film and the 
specimen lie accurately on the same circle. 

The cylindrical body of the camera is provided 


_with a lid (not shown on the photographs) and 


with openings O, and QO2, through which the 
camera can be filled with hydrogen. 


THE BEAM TRAP 


A lead beam trap T is placed about 3 cm in 
front of the focus X. At the bottom of the trap 
is fixed a small fluorescent screen, which serves 
as a useful check of the alignment of the camera. 
It is not practicable to place the beam trap too 
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close to the specimen, in order to reduce air 
scatter, as the broadness of the beam would 
require a large trap. The trap is on a hinge and 
it is convenient to swing it out of the beam, by 
means of an outside knob G, for a short time to 
register the central spot, which serves as a 
starting point for measurements of the spacings. 
Knife edges K, and Kz are provided at either 
end of the window covered by film to determine 
film shrinkage. 


PERFORMANCE OF THE CAMERA 


The camera is most economically used in con- 
junction with a broad x-ray source, so that the 
whole slit system is filled with the radiation. 
This can be obtained by de-focusing a demount- 
able x-ray tube. This lowers the specific loading 
of the target and makes it possible to run the 
tube at a slightly higher power. 

From the geometry of the arrangement, it 
follows that the camera gives a pattern on ex- 
actly twice the scale of an ordinary powder 
camera of the same diameter. 

As the reflections are integrated from a com- 
paratively large area of the specimen, no rota- 
tion or oscillation of the specimen is necessary. 
Typical exposure times are 1 hour with CuK, 
radiation filtered with Ni foil at 45 kv and 
20 ma. 

For smaller Bragg angles, the relative in- 
tensities of diffracted lines should be approxi- 
mately the same as in a conventional powder 
camera; for large Bragg angles, there is a 
strengthening of line intensities due to the x- 
rays penetrating the film at an angle, which is, 
however, compensated to a certain degree by 
weakening of the beams due to traveling a 
longer path through the specimen. The relative 
intensities are also affected by preferred orienta- 
tion in the specimen. This can, however, be 
detected by comparing line intensities on the 
right and left half of the film which are then 
different because the specimen is inclined to the 
beam. The errors in measurements of spacings 
arise from the following factors: 

(1) film shrinkage, which can be corrected by the use 
of the knife edges, 

(2) effective camera radius, which is slightly larger than 
the cylinder radius owing to the finite thickness of the film, 


(3) inaccuracy of setting of the specimen on the common 
circle and inaccuracy of reproducibility of this setting, 
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(4) for thick specimens, the effective depth from which 
X-rays are diffracted varies very slightly with the Bragg 
angle and specimen thickness, and 

(5) the diffracted lines are formed by a superposition 
of many circles, the centers of which all lie along the focal 
line. This tends to shift point of gravity of the diffracted 
lines towards smaller Bragg angles. 


The errors arising from points 2, 3, and 5 can 
be largely eliminated by calibration. The camera 
was calibrated by using standard substances, 
such as quartz, rocksalt, sugar, and gypsum; the 
calibration was extended to smaller Bragg angles 
by using the lower orders of the long spacings 
of a homologous series of silver soaps, which 
give exceptionally sharp reflections. 

The calibration proved that the errors of the 
camera are not serious, or can be eliminated by 


careful experimental technique; indeed, the 
camera seems to be more free of systematic 
errors than conventional cameras. 

After application of corrections, the probable 
error of a single line measurement was found to 
be of the order 1 percent at 20A and 0.1 percent 
at 2A, using CuK, radiation. 
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The Absorption of Incident Quanta by Atoms as Defined by the Mass Photoelectric 
Absorption Coefficient and the Mass Scattering Coefficient 
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Evidence presented indicates that the absorption of quanta as defined by the mass photo- 
electric absorption coefficient +/p and the mass scattering coefficient ¢/p of any element may 


be represented by a simple expression. 


Tabulated values of the calculated mass absorption coefficients u/p are given fog wave- 
lengths from 0.01 to 40 angstroms for hydrogen and carbon and from 0.01 to \x and Xx to Az 
for aluminum, copper, tin, and lead. Graphs are given showing agreement with observed 


values. 


INTRODUCTION 


HE author, in a previous article,’ presented 
an expression for the calculation of ab- 
sorption coefficients. The data were purely em- 
pirical in that every attempt was made to repre- 
sent observed values and, when necessary, to 
compromise differences between observers. It was 
believed, however, that a simple law could be 
derived from the empirical data. This appears to 
have been done satisfactorily, but as might be 
expected from such an approach to the problem 
facts may be indicated before they are fully ex- 
plained. Thus, the quantum states required to 
express the absorption of the elements have been 
found from agreement with observed data. 
In the present article an attempt has been 
made to preserve quantum conditions at the ex- 


"1. A. Victoreen, J. App. Phys. 14, 95 (1943). 
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pense of the observer. It appears that this has 
improved the general agreement with observed 
data and made the individual observer’s errors 
more evident. Examples of this may be seen in 
Fig. 3 for aluminum and Fig. 4 for copper. A 
change in the quantum numbers assigned will 
improve the agreement at wave-lengths ap- 
proaching Ax for copper and aluminum, but on 
close inspection one may possibly be justified in 
assuming the given calculated values to be 
correct. 

The generally excellent data of Allen? are 
shown to contain a large wave-length error on 
wave-lengths below 0.2A and have not been used 
below this in this work. Corrected, but unpub- 
lished data presented to the author by S. J. M. 

? A. H. Compton and S. K. Allison, X-Rays in Theory and 


Experiment (D. Van Nostrand Company, Inc., New York, 
1935), see Allen’s tables. 
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TABLE |, Calculated mass absorption coefficients 
wave-length in angstroms. 











» H 1 C6 Al 13 Cu 29 Sn 50 Pb 82 
0.010 0.1134 0.0571 0.0551 0.0523 0.0490 §=0.0492 
0.012 0.1240 0.0624 0.0602 0.0573 0.0542 0.0563 
0.015 0.1388 0.0694 0.0675 0.0644 0.0620 0.0688 
0.020 0.1577 0.0794 0.0767 0.0739 0.0743 0.0944 
0.025 0.1736 0.0874 0.0845 0.0826 0.0880 0.1298 
0.030 0.1870 0.0941 0.0910 0.0907 0.1040 0.1778 
0.040 0.2083 0.1052 0.1024 0.1073 0.1466 0.3210 
0.050 0.2252 0.1135 0.1112 0.1255 0.2078 0.5403 
0.060 0.2407 0.1215 0.1201 0.1486 0.2942 0.8507 
0.080 0.2628 0.1331 0.1353 0.2102 0.5581 = 1.780 
0.100 0.2794 0.1419 0.1502 0.3018 0.9733 3.163 
0.12 0.2926 0.1494 0.1672 0.4323 1.571 5.010 
0.15 0.3081 0.1592 0.1985 0.7180 2.871 2.040 
0.20 0.3263 0.1741 0.2740 1.495 6.324 4.535 
0.25 0.3366 0.1884 0.3882 2.741 11.49 8.459 
0.30 0.3498 0.2089 0.5573 4.561 18.61 14.05 
0.40 0.3599 0.2588 1.087 10.26 38 46 30.86 
0.50 0.3660 0.3357 1.950 19.25 12.21 55.74 
0.60 0.3736 0.4494 3.226 32.02 20.55 88.60 
0.80 0.3854 0.8119 7.318 70.38 46.37 
1.00 0.3939 1.402 13.95 126.9 86.26 
1.20 0.4009 2.280 23.68 201.4 141.7 
1.50 0.4318 4.256 45.24 45.83 255.6 
2.0 0.4927 9.765 103.6 105.0 $22.2 
2.5 0.5916 18.76 195.6 198.2 856.6 
3.0 0.7375 32.10 326.7 331.2 
4 1.208 74.92 720.7 730.7 
5 1.984 144.3 1304 1321 
6 3.142 246.2 2072 2100 
8 6.902 568.9 279.3 4112 

10 13.09 1082 536.9 6340 
12 22.32 1821 912.7 

15 43.17 3413 1739 

20 101.6 7522 3950 

25 197.6 1357 7377 

30 340.7 2154 

40 804.0 4198 





Allen are in general agreement with Cuykendall® 
and Jones‘ values. An example of this wave- 
length error may be noted in Fig. 6. A somewhat 
similar error may possibly account for the devia- 
tion in the K—L region of copper. These obvious 
discrepancies do not appear for most elements. 

It may appear to some that further elaboration 
of empirical expressions is unnecessary in view of 
developments in the classical theory of radiation, 
but as no simple method of calculating mass ab- 
sorption coefficients accurately from classical 
theory has been made available, there is no 
alternative. Apparently there has been little 
interest in the measurement of absorption coeff- 
cients, however, for the only available observed 
data are old and incomplete. Published values 
often differ by 100 percent. 


PREVIOUS EXPRESSIONS 
Early empirical expressions assumed that the 
mass absorption coefficient for a given element 
was given by 
(u/p)z=Rd’. (1) 
Further investigation indicated that the mass 
absorption coefficient for a given element Z was 


3? Trevor R. Cuykendall, Phys. Rev. 50, 105 (1936). 
4M. T. Jones, Phys. Rev. 50, 110 (1936). 
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roughly given by 
u/p=KZ4d’°. (2) 

All observers have noted that these equations 
given approximate values only under suitable 
conditions. For high values of Z and \ they do 
not give even a good approximation. 

Owen, in 19125 and 1918,° proposed Eq. (2) in 
its exact form, as did Bragg and Pierce.? The 
change in absorption with Z* was later verified by 
Richtmyer and Warburton.* 

The variation of the absorption coefficient with 
\* has been controversial, as may be expected, 
depending upon the wave-léngth under consider- 
ation. Richtmyer,’ in 1925, found that his data 
could be represented by 


(u/p)z=Kv\8 +b. (3) 


For the shortest wave-lengths for which K is 
valid this will be within the limits of experimental 
error, and } will represent the mass scattering 
coefficient. Deviations from the \* law were, how- 
ever, noted by Richtmyer,” Allen," and others." '* 

An expression of form similar to Eq. (2) was 
devised theoretically from J. J. Thompson’s old 
pulse theory" by Compton!" in 1924. 

Compton,!* Kramers,!? Oppenheimer,!* and de 
Broglie'® have suggested expressions based 


TABLE II. Values assigned for Eq. (10) for wave-lengths 
less than Ax. 











Z v1 v 3v3 
ny Ne ny neo ny ne Ss 
1 1 2 1 3 1 ea 0 
6 1 3 1 4 1 2 4 
13 1 3 1 4 1 *% 3 
29 1 4 1 2 1 20 1 
50 1 4 1 0 1 a 1 
82 1 3 1 4 1 2 1 








5 E. A. Owen, Proc. Roy. Soc. 86, 434 (1912). 

6 E. A. Owen, Proc. Roy. Soc. 94, 522 (1918). 

7 Bragg and Pierce, Phil. Mag. 28, 626 (1914). 

8 F. K. Richtmyer and F. W. Warburton, Phys. Rev. 22, 
539 (1923). 

*F. K. Richtmyer, Phys. Rev. 27, 1 (1925). 

10F, K. Richtmyer, Phys. Rev. 18, 13 (1921). 

1S. J. M. Allen, Phys. Rev. 27, 266 (1926). 

2K. W. Wingardh, Zeits. f. Physik 8, 363 (1922); ibid., 
20, 315 (1923). 

13 E. Dershem and M. Schein, Phys. Rev. 37, 1238 (1931). 

44 J. J. Thompson, Conduction of Electricity through Gases 
second edition, p. 326. 

16 A. H. Compton, Phys. Rev. 14, 249 (1919). 

16 A. H. Compton, Phys. Rev. 14, 249 (1919). 

17H. A. Kramers, Phil. Mag. 46, 836 (1923). 

18 J. R. Oppenheimer, Zeits. f. Physik 41, 268 (1927). 

191. de Broglie, J. Physique 3, 33 (1922). 
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Fic. 1. Calculated and observed values of the mass 
absorption coefficient of hydrogen. 


on theoretical considerations. Richtmyer,”*”! 
Jonsson,” and Walter? have discussed general 
absorption formulas ; Houston* and others*® have 
commented on the failure of previous empirical 
formulas to represent adequately experimental 
data. 

The author in a previous article! has shown 
that the photoelectric absorption coefficient 7, of 
any one element is given with excellent accuracy 
by the expression 

tTa=C/v®—D/v'=C'—D"™, (4) 
where C and D are constants which change value 
between absorption discontinuities, and X\ or vy is 
incident wave-length or frequency. 


TABLE ITI. Values assigned for Eq. (10) for wave-lengths 
between Ax and Xz. 


Z v1 v 3v3 
ny Ne ny nN» ny Ne 5 
13 2 3 1 3 2 0 4 
29 2 3 1 3 ? L 4 
50 2 4 1 2 2 2 4 
82 2 5 1 a 2 oo 4 














20 F. K. Richtmyer, Phys. Rev. 27, 1 (1925). 

*1F, K. Richtmyer, Phys. Rev. 30, 755 (1927). 

2 E. Jonsson, Thesis, Uppsala (1928). 

23 B. Walter, Forts. a. d. 35, 929, Geb. der Rontgen 1308 
(1927). 

24 R. A. Houston, Proc. Roy. Soc. 40, 35 (1920). 

2A. H. Compton, X-Rays and Electrons (D. Van 
Nostrand Company, Inc., New York, 1926). 
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Fic. 2. Calculated and observed values of the mass 
absorption coefficient of carbon. 


The constants C and D have been found to be 
given by 


C = vy, v2(3me?/mc), (5) 
D= V1V2V3(37e" ‘mc), (6) 


where v1, v2, and v3 are critical frequencies ob- 
tained without justification from the accepted 
spectral emission formula 


1 1 
v= Re(Z—3)*( -—) 
n\> MN? 
1 1 eos 
+R a2Z—-3)'(—-—)( --)+ | (7) 
n, no* k 4 


Where R is Rydberg’s number, c the velocity 
of light, Z is atomic number, a is Sommerfeld’s 
fine structure constant, s is a screening constant, 
and m, and mz are quantum numbers. 

Satisfactory agreement is obtained by as- 
signing the following values: S=0, $, 1, etc., 
n=1, 2, 3,4, etc., a=(2me?/ch), R=n, so n/k =1. 

It is essential that the bracketed term of 
Eq. (7), which represents the Sommerfeld rela- 
tivity correction, be included to sufficiently high 
powers of a? that the last additive term may be 
neglected. The values of C and D for any element 
can then be obtained by inserting the proper 
atomic number, Z, and assigning suitable values 
for s and n. The constants s and are consistent, 
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TaBce LV. For wave-lengths in angstroms shorter than Ax. 











Zz Cc D Z/A AK 
1 1.274 107? 4.661X10-* 0.99206 ...... 
6 1.2245 1.421107 0.49958 43.5 
13 14.60 0.8362 0.48201 7.955 
29 =—«:178.5 51.73 0.45619 =1.377 
50 940.2 853.8 0.42122 0.4239 
82 4158 11040 0.39573 0.1407 











so that it is possible to assign values to an 
element of atomic number Z which lies between 
two elements Z, and Z2 for which values of s and 
n are satisfactorily known. When C and D are 
satisfactorily known for a given element for any 
one wave-length between discontinuities, all 
other wave-lengths in that region can be cal- 
culated. 

Substituting from Eqs. (5) and (6) for Cand D 
in Eq. (4) and multiplying by No/A, we have 


T Vi¥2 ¥1Vev3\ 3re” No 
a 
p 





og 3 8 re* No 
» 4 3m*c* A 
lt+ta/s2[i+a] 1 
a ( —— log.(1+2a) 
a? i+2a a 


1 1+3a 
+— log.(1+2a) -——| (9) 
2a (1+2a)? 
adding the right hand terms from Eqs. (8) and (9) 
and dividing and multiplying by ?[(8/3)(2e*/ 
m*c*) |, we obtain 








b Vivo =V\VeV3 3 mc* 
n|( a(t 
p y3 y4 4 re’ 
i+ta/2[it+a] 1 
| ( —— log.(1+2a) 
a? 1+2a a 





. ¥§ 1+3a 3/8 wet \ No 
+—loge(1+2a)-———_||-(- \- (10) 

2a (1+2a)*JI4N3 m*%ct7 A 
where a=hy/mc’. 

The term 22[(3/4)(mc*/me?) |] has the dimen- 
sion of a limiting angular frequency, w. The term 
[ (8/3) (we*/m?c*) ]=o, the classical scattering per 
electron. 

It may be noted that the large bracketed part 
of Eq. (10) is now dimensionless, and the mass 
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Tase V. For wave-lengths in angstroms between 











AK and Ax. 
Z Cc D Z/A AL 
13 0.5796 4.294 10-3 0.48201 vandal 
29 14.79 0.8456 0.45619 11.8 
50 106.9 20.87 0.42122 2.75 
82 621.2 352.8 0.39573 0.78 








TABLE VI. Values of ¢.-No—wave-length in angstroms 
maximum value =0.4011. 











nN aeNo r oeNo r aeNo r aeNo 
0.0010 0.0262 0.010 0.114 0.10 0.282 1.0 0.384 
0.0012 0.0300 0.012 0.125 0.12 0.295 1.2 0.388 
0.0015 0.0350 0.015 0.140 0.15 0.310 1.5 0.392 
0.0020 0.0431 0.020 0.159 °0.20 0.329 2.0 0.394 
0.0025 0.0502 ‘0.025 0.175 90.25 0.339 2.5 0.396 
0.0030 0.0566 0.030 0.188 0.30 0.352 3.0 0.397 
0.0040 0.0680 0.040 0.210 040 0.362 4.0 0.397 
0.0050 0.0778 0.050 0.227 0.50 0.367 5.0 0.398 
0.0060 0.0866 0.060 0.243 0.60 0.374 6.0 0.399 
0.0080 0.101 0.080 0.265 0.80 0.382 8.0 0.400 








photoelectric coefficient may be represented by 


T V\V2 ViVoV3 3 No 

— = | es as Jala. 

p y3 y4 4A 
Values of u/p calculated by means of Eq. (10) for 
frequencies beyond a critical absorption fre- 


quency are indicated by the dotted line in Fig. 4 
for copper. 





(11) 


SUMMARY 


The data of Chao***’ and Lauritsen** indicate 
that Eq. (10) gives satisfactory agreement with 
experiment on wave-lengths as short as 0.005A. 
But at very short wave-lengths additional types 
of absorption occur as a result of pair produc- 
tion, etc. 

Equation (10) gives values of u/p which are, on 
the average, within several percent of the values 
given in the original paper ; and where differences 
exist, Eq. (10) appears to favor the observer's 
values. 

Constants for the determination of the mass 
photoelectric absorption coefficient 7/p by Eq. (8) 
have been determined at long wave-length where 
other factors are apparently negligible. The mass 


scattering coefficient, ¢/p, is determined by the 


Klien-Nishina formula, Eq. (9), which becomes 
the determining factor at short wave-lengths. 
Usually the total mass scattering coefficient, 


26 C. Y. Chao, Phys. Rev. 36, 1519 (1930). 

27C. Y. Chao, Science Reports, National Tsing Hua 
University 1, 59 (1932). 

28 J. Reed and C. C. Lauritsen, Phys. Rev. 45, 433 (1934). 


JOURNAL OF APPLIED PHYSICS 

















=| Trin TTT TUTinit Anil 
a I s / 
= ALUMINUM [3 f: / 
— — CALCAAATED J: 
© ALLEN 4: / 
@ CUYKENDALL / H 
1000 f- i} J 
. 1 
ay } / 
100 4 


MASS ABSORPTION COEFFICIENT 


we 
| 


PTT 


Litt 


on ao” 








g 
Litt 





tek.) 0.10 10 
WAVELENGTH IN ANGSTROMS 


Fic. 3. Calculated and observed values of the mass 
absorption coefficient of aluminum. 


a/p, is considered as being the sum of o,/p and 
o./p, where the latter represents the true absorp- 
tion associated with scattering. It appears that 
when o, the Klien-Nishina values of o,, are used, 
the addition of a term a, leads to values of u/p 
which are too high, when compared with experi- 
mental values over a number of elements. We 
may therefore consider that for practical pur- 
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Fic. 4. Calculated and observed values of the mass 
absorption coefficient of copper. 
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Fic. 5. Calculated and observed values of the mass 
absorption coefficient of tin. 


poses from Eq. (10): 
T € te Ge Ge 
> - 2} —F—, 
pp p p p 
The values of »/p given for several elements 
have not been selected as good examples of 
agreement with observed values. No obvious 
deviations are observed on most elements. Values 
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Fic. 6. Calculated and observed values of the mass 
absorption coefficient of lead. 
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of u/p for other elements are being computed and 
will be presented later in tabular form. 

The values of mass absorption coefficients 
previously published by the author are evidently 
a good representation of the probable values for 
practical purposes, but it is believed that the 
values obtained by the expressions given here are 
to be preferred. 

All calculations were carried to five significant 
figures, although only four are shown in the data. 

The author wishes to express his gratitude to 
Mr. William W. Managan for his generous effort 
in expanding Sommerfeld’s original relativity 
correction”® to the higher powers of a. 

CALCULATIONS 

Table | gives the values of u/p for six elements 
calculated from Eq. (10), using the following 
values of the physical constants, taken from 


279A. Sommerfeld, Atomic Structure and Spectral Lines 
(Methuen and Company, Ltd., London), third edition, 
Vol. 1, p. 257. 





Birge :* No, 6.0228 X 108, mole, h, 6.624 X 107’, 
erg sec., e, 4.8025X10-", abs e.s.u., c, 2.9977 
X10", cm _ sec.-!, m, 9.1066X10-*8, g, R, 
1.09737 X 105, cm, a, 7.2976 10-%. 

Table II shows the quantum numbers assigned 
for the several elements for wave-lengths shorter 
than the K critical absorption wave-length. 

Table II] shows the values assigned for wave- 
lengths between Ax and X,;. 

The values of the quantum numbers given in 
Tables II and III have been assigned arbitrarily. 
k is taken as always equal to n. 

The data from Table 1, computed from the 
values in Table 1V and V, are compared with 
observed values in Figs. 1-6. 

Tables IV and V give values to be used in 
Eq. (3). In this convenient form, u/p may be 
calculated for wave-lengths other than given in 
Table |. Table VI gives the values of ¢,No. 


8° R. T. Birge, Rev. Mod. Phys. 13, 233 (1941). 





A Solution of the Problem of Rapid Scanning for Radar Antennae* 


R. F. RINEHART 
Case Institute of Technology, Cleveland, Ohio 


(Received April 30, 1948) 


The problem of devising a rapid scanning focusing antenna based on wave-guide principles 
has been solved only approximately and not too satisfactorily. In the present paper, a known 
property of a plane optical system with variable index of refraction is shown to satisfy the 
requirements of such an antenna system. By a suitable transformation, the plane variable 
refractive index system can be transformed into a surface of revolution of constant refractive 
index which reproduces the optical properties of the plane system. This surface provides an 
apparently practicable solution to the scanning problem. It is shown that the diameter of the 
feed circle can be made arbitrarily small by judicious use of an ellipsoidal reflector. 


INTRODUCTION 


HE problem of designing a rapid scanning 

radar antenna employing the theory of 
parallel plate wave guides has been discussed by 
S. B. Myers.'! The fundamental property of such 
a curved wave guide is a Fermat least time 
principle: Electromagnetic energy is propagated 
along geodesics (i.e., shortest distance curves) of 
the hypothetical mean surface between the parallel 
conducting plates. The rapid scanning problem 





* This paper is a result of work performed on Contract 
No. W28-099-ac-141 for the United States Air Forces under 
the sponsorship of the Electronic Research Laboratory, 
Air Materiel Command. The particular results treated 
here grew out of a suggestion of Dr. Wade Ellis of the 
Watson Laboratories. 

1S. B. Myers, ‘Parallel plate optics for rapid scanning,” 
J. App. Phys. 18, 221 (1947). 
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for microwaves, as educed in Myers’ paper, is to 
find a surface = (to play the role of mean surface 
of a parallel plate wave guide) such that 

(a) = contains a circle (the feed circle), 

(b) = contains a straight line segment (the 
aperture), 

(c) (Optical requirement) all geodesics joining 
a fixed point of feed circle to points of the 
aperture meet the aperture at a fixed angle, 

(d) (illumination requirement) all geodesics 
normal to the feed circle pass through a fixed 
point (the center of illumination) of the aperture. 

Myers points out that it is not known whether 
a surface with such properties exists. If, however, 
requirement (d) is relinquished, a solution to the 
modified problem is the so-called RCA solution, 
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constructed from a portion of a frustum of an 
appropriate cone.' This solution furnishes un- 
satisfactory illumination characteristics at the 
aperture. 

The purpose of this paper is to describe a 
solution surface satisfying requirements (a), (b), 
and (c) and which yields what seems to be 
satisfactory illumination. 


A PERTINENT RESULT FROM OPTICS 


The solution stems from an appropriate adap- 
tation, or transformation, of the following known 
optical result 2? 

In the optical system composed of a unit circle 
with variable index of refraction »=(2—r°)}, 
where ¢ is the distance from the center, imbedded 
in a plane medium of refractive index 1, all light 
rays entering the circle from any fixed point P 
on its periphery energe from the circle parallel 
to the diameter OP * (see Fig. 1). Imagine for 
the moment that this unit circle is a mean 
surface for a variable refractive index wave 
guide. The unit circle would serve as feed circle. 
Any straight segment in the plane, not inter- 
secting the circle, could serve as aperture with 
respect to the optical requirement. As P rotates 
through a semicircle, the emergent beam scans 
180°. It may be remarked that this system also 
satisfies the illumination requirement with re- 
spect to a diameter of the circle as virtual 
aperture. 

It can be shown that the angular width of a 
beam from P is equal to the angle subtended at 

2 See, for instance, R. K. Luneberg, Mathematical Theory 
of Meory (Brown University Lecture Notes, 1944), pp. 


3 Note that » is continuous at r=1, and hence that 
ncident and excident angles are equal at the unit circle. 
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Fic. 3. Use of elliptical re- 
flector to reduce feed circle 
Fic. 2. diameter. 


0 by the emergent beam. From this it follows 
that the illumination intensity along a cross 
section of the emergent beam varies as seca. 
Thus a 90° beam width at P, for instance, 
would provide an emergent beam of linear width 
v2, and the illumination intensity would vary 
only 41 percent from its value at Q. This would 
seem to be an acceptable illumination pattern. 
Such a system has the obvious disadvantage 
that it can utilize only a portion of each revolu- 
tion (unless two opposite reflectors are used). 
On the other hand, this permits the feed to be 
active for an arc of only 90° in order to get a 
complete 90° scan by the emergent beam. 


THE TRANSFORMATION TO A SYSTEM OF 
CONSTANT REFRACTIVE INDEX 

We seek now a surface which duplicates the 
optical characteristics of the variable refractive 
index (v.r.i.) system just described. The straight 
wave front RS in the v.r.i. system results from 
the selective time delays induced in the various 
rays by the variable refractive index, which 
slows down the inner rays traveling the shorter 
distances. We seek to accomplish this delay by 
varying the path lengths, rather than the ve- 
locities, of the rays. We propose to achieve this 
effect by ‘‘bellying up’’ the unit circle into an 
appropriate surface of revolution. 

To this end we seek a surface = generated by 
revolving about the z axis a curve p= p(s), p being 
the arc length measured from the axis, which 
will reproduce the optical properties of the v.r.i. 
system (see Fig. 2). In the v.r.i. system, the 
differential of time (in appropriate units) along 
a curve is given by 


dt? = wd? = w*(dr?+rd@), (1) 
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where df is the differential of arc length in the 
plane of Fig. 1, and (7,6) are polar coordinates. 

Along a curve of the surface 2, the differential 
of time, in the same units, is precisely the differ- 
ential of arc length dé, since the refractive index 
is 1, and is given by 


df? =d?=ds*+p'd@, (2) 


where 2, p, 6 are cylindrical coordinates. 

We desire to so choose the function p(s) that 
dt as given by (2) will be identical with dé as 
given by (1). Let us introduce a new (to 2) 
variable r by means of the relation 


ds =[p(s)/r dr. (3) 
(2) then takes the form 
de = p?/r?(dr?+1r°d6*) = (p?/r*) dg? (4) 


which will be identical with (1) provided the 
relation 


p/r=p(r) =(2—r°*)! (S) 
is imposed. From (5) follows, since 0=r=1, 
r=(1+(1—p?)!}!. (6) 


Relations (3) and (6) provide a differential 
equation in s and p, which, with the boundary 
condition, s=0 when p=0, defines uniquely the 
functional relation between s and p which makes 
dt =dé. 

From Eqs. (5) and (6) we note that when 
p=0, r=0 and when p=1, r=1. Hence if we 
cut off the surface at p=1 to obtain our surface 
z, Eq. (5) defines a map of = on the unit circle 
in the (7,6) plane of the v.r.i., in which 


and 6-8. 


Now any curve arc which minimizes /fdt in 
the v.r.i. system will have an image on = which 
minimizes fdt, and conversely. Hence light 
rays in the v.r.i. system have images which are 
’ geodesics of 2, and conversely. 

Furthermore, Eq. (4) shows that this map is 
conformal (with magnification ratio p/r) and 
hence angles between curves are preserved under 
the mapping. Thus if a ray in the v.r.i. system 
_ leaves P(1,0;) at an angle ¢ with the unit circle, 
and recontacts the unit circle at M(1,@2) at an 
angle y with the unit circle, then the image 
geodesic of = behaves in precisely this same way 
with respect to the image circle p=1. 

Now let the surface = be continued planewise 


p—r 
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beyond p=1. It is well known that a geodesic 
makes equal angles, approaching and leaving, 
with an edge of a surface. It follows that the 
surface with the planar continuation for p=1 
duplicates completely the optical properties of 
the v.r.i. system, and constitutes a solution of 
the scanning problem. 


PROPERTIES OF THE SURFACE 


The equation of the surface = is readily 
computed from (3) and (6) to be 


s=1/2(p+sin~p). (7) 


Its equation in cylindrical coordinates can be 
computed in terms of elementary functions. It is 


2=1/2[4—3p?+4(1—p*)?}} 
—1/v3 log((3/2{1+(1—p*)#} } 
+(1/2{3(1—p*)'+1} }) 
—1/2+1/2 v3 log(2+v3). (8) 


The surface has a horizontal tangent plane at 
p=0 and a vertical tangent plane at p=1. The 
surface is concave down everywhere between 
p=0 and p=1. 


REDUCTION OF SIZE OF FEED CIRCLE 


The parallel plate system described here seems 
practicable from every standpoint except the 
large size of the feed circle, relative to the emer- 
gent beam width, or effective aperture. It can 
be proved that this difficulty is geometrically 
insurmountable in any surface of revolution. 
Specifically, the linear width of the emergent 
beam can never exceed the diameter of the feed 
circle (or of any intervening parallel circle of the 
surface). Consequently, in the solution given here 
the feed circle is the smallest attainable. It is 
possible, however, to reduce the size of the feed 
circle by optical means. A reflector in the shape 
of a section ‘of an ellipsoid can be utilized to 
produce a virtual point source on the unit circle 
from an actual point source situated much closer, 
or even on the axis of the surface. As the actual 
feed F, (Fig. 3) rotates on a small circle, the 
virtual, or image, feed F, describes the unit 
circle. Rough computations indicate that a 
reflector yielding a suitable angular beam width 
could be designed to be housed beneath the 
dome-shaped wave guide. With this refinement 
the antenna system described seems to furnish a 
workable solution.to the scanning problem. 
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Letters to the Editor 








Electron Optical Observation of 
Magnetic Fields 


L. MARTON 
National Bureau of Standards, Washington 25, D.C. 
June 21, 1948 


N a previous communication! it was pointed out that 
the electron optical ‘‘Schlieren’”’ effect can be used for 
the observation of minute magnetic and electrostatic 
fields. Reference was also made to the use of that method 
for the purpose of quantitative measurement of the fringe 
field, extending from magnetized regions of a steel wire, as 
used for magnetic recording purposes. In pursuing observa- 
tions on the same subject we used a center stop catching 
all direct radiation and allowing only the deflected elec- 
trons to reach the fluorescent screen. This center stop con- 
sisted of a 1-mm diameter thin copper disk supported by 
two silver cross-wires of 0.033-mm diameter. In one of our 
experiments it so happened that the image of the center 
stop (produced by projection from the virtual source 
between the lens and the stop) was displaced laterally by 
means of an external deflecting magnetic field. By adjusting 





Fic, 1. 





Fic. 2 








Fic, 4. 


the image of one of the supporting silver wires to be 
parallel with the image of the magnetized wire, we found 
that the image of the silver wire is not a straight line when 
in the neighborhood of the image of the magnetized wire. 
Figures 1-3 show a sequence in which the two images are 
first completely separated, then moved closer, and finally 
made to coincide. It is found that in the case of the two 
images, the image of the thin silver wire bulges out every- 
where where magnetized regions of the steel wire deflect 
the incoming electron beam. 

A schematic explanation of this observation is given by 
Fig. 4. The magnetic lens L used for forming a somewhat 
magnified image of the magnetized wire M forms a reduced 
image of the distant source S at the point P. This reduced 
image of the source serves as a virtual source forming a 
shadow image of the center stop C and of its supporting 
wires. If the wire M has magnetized regions, however, the 
electron beam coming from the source will be deflected at 
each such magnetized region with the result that the virtual 
source P will be displaced sidewise. Such a displacement 
produces the images as observed, and this displacement 
is obviously a function of the strength of the magnetic 
field. 

In close analogy to the electron optical ‘Schlieren’ 
effect, the displacement of that image can be used for a 
quantitative measurement of the magnetic field strength. 
For this purpose it is assumed, as before, that the magnetic 
field of the wire is produced by an ideal dipole and thus the 
field obeys the inverse cube law. The absolute value of the 
field can be calculated at least to a first approximation at 
the equatorial plane of the dipole from the measured 
deflection of the image and from the known geometrical 
constants of the apparatus. Calculations by Mr. S. H. 
Lachenbruch (to be published later) indicate that, for a 
commercially available magnetic recording wire, a record- 
ing, produced by 40-ma exciting current of the recording 
head, produces a magnetic field of 100 gauss at 0.1-mm 
distance from the wire axis. These data are based on 
measurement carried out on Fig. 5. 





Fic, 3. 
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It is to be pointed out that this method offers, in some 
cases at least, the distinct advantage over the “Schlieren” 
observation that the displacement of the wire image can 
be measured with greater accuracy than the intensity dis- 
tribution across the “Schlieren’”’ pattern. 

An extension of this method consists of placing a net- 
work of cross-wires beyond the back focal plane of the 
imaging lens (at the position C of Fig. 4). Such a network 
enables us to record a complete mapping of extended fields 
in a single operation. Experiments are being undertaken 
here to apply this method to numerous problems. 


'L. Marton, J. App. Phys. 19, 687 (1948). 





An Observation of Corrosion Cracking 
without Stress 


W. D. ROBERTSON AND H. H. UHLIG 


Department of Metallurgy, Massachusetts Institute of Technology, 
Cambridge, Massachusetts 


June 21, 1948 


TRESS corrosion cracking of metals and alloys is one 
of the most important unsolved problems in the field 
of corrosion. The reasons are largely unknown why metals 
subjected to residual or applied stress suddenly crack after 





a period of time in specific media. Various suggestions of 
mechanism have been made,'~ all pointing to the presence 
of metal or non-metal acting as anode at the base of a 
growing crack with the wall functioning as cathode. 
However, the detailed mechanism is not entirely clear, 
particularly with respect to transcrystalline cracking and 
to the extreme selectivity of the environment in which 
cracking occurs. The necessity of a specific environment 
for cracking indicates that the explanation based alone on 
galvanic effects of stressed metal may be inadequate and 
too general since, presumably, many conducting environ- 
ments should thereby cause stress corrosion cracking. Such 
is not the case. Furthermore, non-metals such as plastics 
exhibit stress cracking phenomena in various organic 
solvents where an electrochemical mechanism is obviously 
irrelevant. , 

Refined experiments on stress corrosion cracking or 
“season cracking”’ of brass, although not pointing to a 
specific mechanism, have indicated the relative importance 
of some factors.‘ Thus, impurities do not play a large role, 
and the addition of other elements does not greatly reduce 
susceptibility; furthermore, the critical stress necessary 
for cracking in an ammonia atmosphere has been reduced 
to lower and lower limits. Single crystals of brass in 
ammonia also crack by stress corrosion, so the grain 
boundary in polycrystalline alloys is relieved of total 
responsibility for failure. One is tempted to conclude from 
these facts that corrosion cracking in some alloys may be 


a2 


Fic. 1. Crystallographic disintegration of Mg2Sn single crystal. 
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associated with metallurgical factors inherent to the alloy 
itself. 

An observation we recently made, in connection with a 
study of the chemical properties of the intermetallic com- 
pound Mg)Sn, lends weight to this possibility. The com- 
pound was prepared in an argon atmosphere from sublimed 
magnesium (courtesy of the Dow Chemical Company) and 
from electrolytically purified tin. Large crystals were 
produced from the liquid by cooling slowly through the 
melting point to room temperature at a rate of 20°C per 
minute. The resulting cast rods were 3” in diameter and 
could be fractured (cleaved) easily to produce a bright 
bluish colored, essentially plane cleavage face across the 
entire diameter of the rod. 

These crystals react readily with oxygen-free sodium 
chloride solutions to form magnesium hydroxide and finely 
divided tin or, in the presence of oxygen, stannous oxide. 
In water, however, the compound is comparatively stable, 
and in a moist atmosphere (30 percent relative humidity) 
the cleavage planes remain bright for a period of weeks. 

Immersed in ordinary or triple distilled water for 
several hours, single crystals of the compound were ob- 
served to crack along definite planes, forming small geo- 
metrical fragments. Under the low power microscope, the 
fragments resolve themselves into three-sided pyramids, 
rhombohedra, and equilateral triangles, all of which can 
be derived from the octahedral cleavage planes of the 
crystal system involved. A typical group of crystal frag- 
ments, obtained by the action of distilled water on a single 
crystal of MgeSn, is shown in Fig. 1. Similar crystallo- 
graphic breakdown occurs in dilute sodium hydroxide, 
at a rate increasing with the hydroxyl ion concentration, 
in hydrogen peroxide (35 percent), and in water from 
which dissolved oxygen has been removed by boiling. 
Furthermore, the fragments continue to subdivide if 
allowed to remain in contact with water. 

The crystal structure of Mg2Sn is that of CaF» with the 
magnesium ions occupying the fluorine positions, and tin 
ions the positions of the calcium ions. The {111} planes are 
the cleavage planes across which cohesion presumably is 
least and are the planes of dense packed magnesium and 
tin ions, with the magnesium ions slightly elevated above 
the plane of tin ions. The action of water in splitting the 
erystal, although possibly related to traces of impurities 
in the compound, is more likely one of reaction or hydration 
of ions in the {111} planes, thereby destroying the bonds 
between these planes. Finite fragments are obtained, pos- 
sibly owing to the fact that the crystal is not perfect, and 
water, therefore, is drawn in between certain preferred 
planes rather than into all planes equally. Allied with this 
mechanism, the formation of a corrosion product between 
cleavage planes may result in a pressure that mechanically 
splits the crystal. 

The presence of hydrogen as a corrosion product does 
not appear to be involved, insofar as a crystal mounted as 
cathode in 0.02N NaOH (approximately 0.05 amp./cm*) 
fractured much less rapidly than another crystal, mounted 
in the same cell, as anode. 

The slow cooling of the Mg2Sn precludes large residual 
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stresses, and the brittle character of the compound is such 
that large stresses could not be supported. Therefore, the 
fracturing appears to be an instance of transcrystalline 
corrosion cracking in the absence of appreciable stress. 

The analogy with cracking in certain stressed mag- 
nesium or aluminum alloys in presence of water is obvious. 
The mechanism may have its counterpart in other metallic 
systems, with hydration and solvation supplementing 
galvanic action within the crack. Once a crack starts, 
application of tensile stress accelerates propagation of the 
crack and eventual failure of the alloy. 

The authors are indebted to the Office of Naval Research 
for support of this investigation. 

1 E. H. Dix, Jr., Trans. A.I.M.E. 137, 11 (1940). 


2J. T. Waber, H. J. McDonald, and B. Longtin, Trans. Am. Elec- 
trochem. Soc. 87, 209 (1945). 

3R. B. Mears, R. H. Brown, and E. H. Dix, Jr., “Symposium on 
stress corrosion cracking of metals,"” A.S.T.M. and A.I.M.E. ,Phila- 
delphia, 323 (1945). 

i ulow and G. Edmunds, “Symposium on stress corrosion 
to of metals,’’ A.I.M.E. and A.S.T.M., Philadelphia, 19 and 67 





New Booklets 








Tracerlab, Inc., 55 Oliver Street, Boston 10, Massa- 
chusetts, has issued its Catalog A, having 40 pages. 
Available on request. 


James G. Biddle Company, 1316 Arch Street, Philadel- 
phia 7, Pennsylvania, has announced a new bulletin, No. 
24-25, on its ‘“‘Ducter’’ low resistance ohmmeter, a hand- 
book on low resistance measurements in the shop and in 
the field. 12 pages. Available on request. 


Eberbach and Son Company of Ann Arbor, Michigan, 
offers free copies of its latest Announcer of Scientific Equip- 
ment, No. 48-6-31. The issue features an article on ‘‘De- 
velopment of chemical porcelain.” 


Burrell Technical Supply Company, 1942 Fifth Avenue, 
Pittsburgh 19, Pennsylvania, has published bulletin No. 
214 describing combustion tubes and boats and featuring 
a new tube, the “Zircum” combustion tube. 4 pages. 


T. C. Wheaton Company, Millville, New Jersey, has 
published a bulletin illustrating and describing its line of 
electronic glassware. 4 pages. 


Allen B. Du Mont Laboratories, Inc., 1000 Main Ave- 
nue, Clifton, New Jersey, announces the fifth edition of its 
Catalog of Equipment for Oscillography. 176 pages. 


Ward Leonard Electric Company, Mount Vernon, New 
York, has issued its Catalog D-130 on resistors, rheostats, 
and radio amateur relays. 24 pages. A copy may be pro- 
cured by writing to Radio and Electronic Distributor 
Division, Ward Leonard Electric Company, 53 West 
Jackson Boulevard, Chicago 4, Illinois. 
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National Bureau of Standards has prepared and made 
available as NBS Applied Mathematics Series 2 a 20-page 
Table of Coefficients for Obtaining the First Derivative without 
Differences. May be obtained only from the Superintendent 
of Documents, U. S. Government Printing Office, Washing- 
ton 25, D. C., at 15 cents per copy. 


The Scientific Instrument Manufacturers’ Association of 
Great Britain, Ltd., 26 Russell Square, London W. C. 1, 
has published the 1948 edition of its handbook of electronic 
instruments, which gives a brief review of the products of 
members of the Electronics Section of the Association. 
Copies available upon request. 28 pages. 
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New Appointments 


Haldon A. Leedy, Acting Director of Armour Research 
Foundation of Illinois Institute of Technology, has been 
named Director. 

William J. Youden, mathematical statistician and 
chemist, has been appointed to the staff of the National 
Bureau of Standards as Assistant Chief of the Statistical 
Engineering Section. 


L. B. Spinney Honored 


Professor L. B. Spinney, physics textbook author and 
head of the department of physics at Iowa State College, 
was honored recently by the Iowa State College Alumni 
Association when he was awarded a faculty citation in 
recognition of long, outstanding, and inspiring service to 
the college. Professor Spinney began his teaching duties at 
Iowa State in 1893. 


Necrology 


E. F. Burton, Chairman of the Department of Physics 
of the University of Toronto, passed away on July 6. He 
‘was outstanding for his pioneer work in electron microscopy 
on this continent. 

Professor Emeritus Ernest George Merritt, first dean of 
the Graduate School at Cornell University and a member 
of the Cornell Physics Department for 45 years until his 
retirement in 1935, died June 5 at the age of 83. 


Bulletin of Mathematical Biophysics 


Following is the table of contents for the September 1948 
issue of the Bulletin of Mathematical Biophysics: 


“Physical aspects of organic evolution’’—Alfred J. 
Lotka 

“Theory of the measurement of blood flow by the 
dilution of an indicator’’—John L. Stephenson 

“Interaction of cycles’’—J. B. Roberts 

“An analysis of theoretical systems of differentiating 
nervous tissue’’—Alfonso Shimbel 

“Cycle distributions in random nets’’—Anatol Rapo- 
port 

“On _ periodicities in 
Rashevsky 

“The kinetics of blood coagulation’’—John Z. Hearon 

“A note on the units of membrane permeability to 
water’’—H. D. Landahl 

“On the theory of blood-tissue exchange of inert gases: 
VI. Validity of approximate uptake expressions’’— 
Manual F. Morales and Robert E Smith 


metabolizing systems’’—N. 


The University of Chicago Press, Chicago 37, Illinois, 
Volume 10, Number 3. 


Conference at Oak Ridge 


The physics division of Oak Ridge National Laboratory 
sponsored a conference on ‘‘Nuclear Physics and Low 
Temperatures” the weekend of August 7 and 8. The con- 
ference, which was the first of a projected series of physics 
meetings was organized in response to the increasing 
interest in nuclear problems involving low temperature 
phenomena. Speakers and their subjects were as follows: 


Alignment of nuclear spins—M. E. Rose, Oak Ridge National 
Laboratory. 

Nuclear relaxation phenomena at low temperatures—E. M. Purcell, 
Harvard University. 

Thermal equilibrium at low temperatures—H. B. G. Casimir, 
Philips’ Works, Eindhoven, Netherlands. 

The superfluid state—F. London, Duke University. 


High Polymer Meetings at Brooklyn 


The following are among future meetings on high 
polymers to be sponsored by Polytechnic Institute of 
Brooklyn, 99 Livingston Street, Brooklyn 2, New York: 
Symposium on Ultrasonics and Macromolecules, October 
23, 1948, Brooklyn; Society of Rheology meeting, probably 
November 5-7, 1948, New York. 


U. S. Testing Company Expands 


In entering the field of supersonic research, the United 
States Testing Company, Hoboken, New Jersey, recently 
announced the completion of arrangements to take over 
in specific fields the research and development activities 
of the Crystal Research Laboratories of Hartford, Con- 
necticut. 
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